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Abstract 



Decoherence of many-spin systems in NMR: 
From molecular characterization to an environmentally induced 
quantum dynamical phase transition 

The control of open quantum systems has a fundamental relevance for fields rang- 
ing from quantum information processing to nanotechnology. Typically, the system 
whose coherent dynamics one wants to manipulate, interacts with an environment that 
smoothly degrades its quantum dynamics. Thus, a precise understanding of the inner 
mechanisms of this process, called "decoherence", is critical to develop strategies to 
control the quantum dynamics. 

In this thesis we solved the generalized Liouville-von Neumann quantum master equa- 
tion to obtain the dynamics of many-spin systems interacting with a spin bath. We also 
solve the spin dynamics within the Keldysh formalism. Both methods lead to identical 
solutions and together gave us the possibility to obtain numerous physical predictions 
that contrast well with Nuclear Magnetic Resonance experiments. We applied these 
tools for molecular characterizations, development of new numerical methodologies and 
the control of quantum dynamics in experimental implementations. But, more impor- 
tant, these results contributed to fundamental physical interpretations of how quantum 
dynamics behaves in open systems. In particular, we found a manifestation of an envi- 
ronmentally induced quantum dynamical phase transition. 
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Resumen 



Decoherencia en sistemas de espines interactuantes en RMN: 
De la caracterizacion molecular a una transition de fase en la 
dinamica cuantica inducida por el ambiente 

El control de sistemas cuanticos abiertos tiene una relevancia fundamental en campos 
que van desde el procesamiento de la information cuantica hasta la nanotecnologia. 
Tipicamente, el sistema cuya dinamica coherente se desea manipular, interactiia con un 
ambiente que suavemente degrada su dinamica cuantica. Es asi que el entendimiento 
preciso de los mecanismos internos de este proceso, llamado decoherencia, es crftico para 
el desarrollo de estrategias para el control de la dinamica cuantica. 

En esta tesis usamos la ecuacion maestra cuantica generalizada de Liouville-von 
Neumann para resolver la dinamica de sistemas de muchos espines interactuando con 
un bano de espines. Tambien obtuvimos la dinamica de espines dentro del formalismo 
de Keldysh. Ambos metodos nos llevaron a identicas soluciones y juntos nos dieron la 
posibilidad de realizar numerosas predicciones que concuerdan con las observaciones de 
experimentos de Resonancia Magnetica Nuclear. Estos resultados son usados para la 
caracterizacion molecular, el desarrollo de nuevas metodologias numericas y el control 
de la dinamica cuantica en implementaciones experimentales. Pero aun mas importante 
es el surgimiento de interpretaciones fisicas fundamentales de la dinamica cuantica de 
sistemas cuanticos abiertos, tales coma la manifestation de una transition de fase en la 
dinamica cuantica inducida por el ambiente. 
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Chapter 1 
Introduction 



Quantum Mechanics was developed to describe the behavior of matter at very small 
scales, around the size of single atoms. Today, it is applied to almost every device that 
improves our quality of life, from medical to communication technology. Since it involves 
laws and concepts that challenge our intuition, it keeps having a revolutionary impact 
on the formulation of new philosophical and scientific concepts not totally solved today 
|Omn92l ISch04| . While the foundations of quantum mechanics were established in the 
early 20 th century, many fundamental aspects of the theory are still actively studied and 
this thesis intends to contribute to this knowledge. 

1.1 What is quantum physics? 

One of the main characteristics of quantum mechanics is that it involves many coun- 
terintuitive concepts such as the superposition states. They were illustrated by the 
Austrian physicist Erwin Schrodinger in 1935 by his famous Schrodinger's cat thought 
experiment. In his words |Sch35] : 

"One can even set up quite ridiculous cases. A cat is penned up in a steel 
chamber, along with the following device ( which must be secured against direct 
interference by the cat): in a Geiger counter there is a tiny bit of radioac- 
tive substance, so small, that perhaps in the course of the hour one of the 
atoms decays, but also, with equal probability, perhaps none; if it happens, the 
counter tube discharges and through a relay releases a hammer which shatters 
a small flask of hydrocyanic acid. If one has left this entire system to itself 
for an hour, one would say that the cat still lives if meanwhile no atom has 
decayed. The psi-function of the entire system would express this by having 
in it the living and dead cat (pardon the expression) mixed or smeared out in 
equal parts. 

It is typical of these cases that an indeterminacy originally restricted to the 
atomic domain becomes transformed into macroscopic indeterminacy, which 
can then be resolved by direct observation. That prevents us from so naively 
accepting as valid a "blurred model" for representing reality. In itself it would 
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not embody anything unclear or contradictory. There is a difference between 
a shaky or out-of-focus photograph and a snapshot of clouds and fog banks. " 

Erwin Schrodinger 




Figure 1.1: Cartoon description of the Schrodinger 's cat paradox. After an hour 
the cat is in a quantum superposition of coexisting alive and dead states. Only af- 
ter opening the box we found the cat in a defined state. Figure extracted from 
http:/ /en.wikipedia.org/wiki/Image:Katze.jpg. 

Essentially, he states that if we put an alive cat in a box where, isolated from external 
interference, is in a situation where death has an appreciable probability, the cat's state 
can only be described as a superposition of the possible state results (dead and alive), 
i.e. the two states at the same time. This situation is sometimes called quantum 
indeterminacy or the observer's paradox: the observation or measurement itself affects 
an outcome, so that it can never be known what the outcome would have been, if 
it were not observed. The Schrodinger paper |Sch35j was part of a discussion of the 
Einstein, Podolsky and Rosen's paradox [EPR35J that attempted to demonstrate the 
incompleteness of quantum mechanics. They said that quantum mechanics has a non- 
local effect on the physical reality. However, recent experiments refuted the principle of 
locality, invalidating the EPR's paradox. The property that disturbed the authors was 
called entanglement (a superposition phenomenon) that could be described briefly as a 
"spooky action at a distance" as expressed in ref. [EPR35J. This was a very famous 
counterintuitive effect of quantum mechanics which leads very important physicists to 
mistrust of quantum theory. The entanglement property could be schematized by adding 
some condiments to the Schrodinger's cat thought experiment. First of all, we may 
consider that the indeterminacy on the cat's state is correlated with the state of the 
flask of hydrocyanic acid, i.e. if the cat is alive the flask is intact but if the cat is dead 
the flask is broken. We have here two elements or systems (the cat and the flask) in 
a superposition state |cat alive, flask intact) and |cat dead,flask broken) existing at the 
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same time. Assuming that after an hour we can divide the box with a slide as shown 
in figure LI and deactivate the trigger, we can separate as we want the two boxes. 
Then, if someone opens the cat's box and sees the cat's state, the state of the flask will 
be determined instantaneously without concerning the distance between them. This is 
only a cartoon description of what quantum entanglement is about, but for a further 
description we refer to Nielsen and Chuang (2000) |NC00] or chapter [6j 

One of the most interesting effects of quantum superposition is the interference phe- 
nomenon consequence of the information indeterminacy of the quantum state (dead or 
alive). The famous double slit ideal experiment, as Richard Feynman said, contains 



everything you need to know about quantum mechanics. As shown in fig. 1.2 a), the 
experiment consists of a double slit where a particle (photon, electron, etc.) can pass and 
a screen where it is detected. Behind it, there is a screen where we can register where 
the particle arrives. If only one of the slits is open, we have certainty that the particle 
only can pass through this slit. The probability to arrive to different places of the screen 
is shown in figure 1.2 b). There, we see that the most probable place for the particle 
arrival is obtained projecting the center of the slit to the register screen. Moving away 
from it, the probability decreases monotonically. The reciprocal situation occurs if only 
the other slit is open. However, if we leave the two slits open an interference pattern 



appears as in figure |L2 c). Figures 1.2 b) and c) represent mathematical probabilities 
(mathematical reality) describing the physical reality shown in figure 1.2 d) |TEM + 89j . 



Paul Kwiat, Harald Weintfurter and Anton Zeilinger making reference to quantum 
interference, in ref. |KWZ96j . express: 

"According to the rules of quantum mechanics, interference occurs whenever 
there is more than one possible way for a given outcome to happen, and the 
ways are not distinguishable by any means (this is a more general definition 
of interference than is often given in textbooks). In the double-slit experiment, 
light can reach the screen in two possible ways (from the upper or the lower 
slit), and no effort is made to determine which photons pass through which 
slit. If we somehow could determine which slit a photon passed through, 
there would be no interference, and the photon could end up anywhere on the 
screen. As a result, no fringe pattern would emerge. Simply put, without two 
indistinguishable paths, interference cannot occur. " 

Paul Kwiat, Harald Weinfurter and Anton Zeilinger 

Thus, the quantum mechanics is the physics of potentialities. When we have deter- 
minacy of some event, the classical physics appears. Nowadays, this appearance of the 
classical physics and state determinacy is considered a consequence of a phenomenon 
called decoherence |Zur03] which is the central topic of this thesis. 



1.2 Decoherence: the degradation of quantum su- 
perpositions 

The gedanken experiments introduced above must involve a perfect shielding from ex- 
ternal influences allowing the existence of quantum superposition. Realistic quantum 
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Figure 1.2: The double slit experiment, a) Schematic representation of the double slit 
device, b) The solid line is the classical probability prediction which is the sum of the 
individual one-slit probabilities, c) Interference pattern predicted by quantum proba- 
bilities accounting superposition. Panels b) and c) describe mathematical probabilities, 
in panel d) the physical reality is shown. The experiments were performed by A. Tono- 
mura, et al. |TEM + 89] where they showed that single electron events build up to form 
an interference pattern in the double-slit experiments. The number of electrons detected 
are 100, 3000, 20000 and 70000 respectively. Reprinted with permission from A. Tono- 
mura, et al. Amer. J. Phys. Vol. 57, Issue 2, Page 117, 1989. Copyright 1989, American 
Association of Physics Teachers. 
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systems are never isolated, because they are immersed in an environment that continu- 
ously interacts with them. A typical environment consists of a system with many degrees 
of freedom that are hardly fully controlled or are not relevant for the observation. The 
system-environment (SE) interaction degrades the quantum superposition leading to the 
phenomenon called decoherence |Zur03l ISch04j . Actually, the measurement process to 
observe if the cat is dead or alive involves an interaction between the system (cat, acid, 
box, hammer, etc.) and the environment (observer, apparatus to observe, etc.). When 
the observation is performed, the cat is found either dead or alive, but not in the state 
dead and alive. The last one is a pure-state and the first one is a mixed-state. The 
decoherence process leads the system from a pure-state to a mixed-state. It is impor- 
tant to emphasize that, although quantum mechanics is open to many interpretations, 
decoherence by itself is neither an interpretation nor a modification of the theory. Thus, 
their existence can be taken as a well-confirmed fact. However, the implications that 
derive from decoherence could need some interpretations and this is one of the reasons 
why nowadays many researchers are devoted to its study |Zur03t ISch04] . 

Decoherence does not exist if we consider the entire system. It arises when we are 
interested in a particular part of the system leading to the consideration of a system 
plus an environment which is called an open system. Looking at the properties of 
the system, the environment modifies them leading to decoherence. It is at this point 
when the concept of the reduced density operator appears as a tool to mathematically 
describe the quantum world. A system is described by an entity called density operator, 
but the density operator of the Universe is impossible to obtain, thus one decides to 
reduce it to describe a relevant subsystem. The concept of the reduced density operator 
appeared together with quantum mechanics introduced by Lev Landau 1927 |Lan27] and 
further developed by John von Neumann 1932 |Neu32j and W.H. Furry 1936 [Fur36j. 
To illustrate the idea of how the reduced density matrix works, and why by observing 
at a subsystem we can not distinguish between a pure and a mixed-state, we consider a 
system with two entangled elements in a pure-stata^ 

l*> = ^(i+>xi->2-i-)li+> 2 )- (1-1) 

For an observable O that belongs only to the system 1, i.e. O — 0\ ® I2, the expectation 
value is given by 

(d) # = (1.2) 
where the density operator of the pure-state is defined by 

P=l*)(*l- ( L3 ) 

This statistical expectation value is defined as the sum of the values of the possible 
outcomes, multiplied by the probability of that outcome. The same statistics is applied 
to the reduced density operator that is obtained by tracing over the degrees of freedom 
of the system 2. Thus, we obtain 

(d)^ = Tr |p6} = Tn {axOx} , (1.4) 



1 This entanglement is consequence of a previous interaction between the two elements. 
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where the reduced density operator is 

ft = Tr 2 {|*> (*|} = a (*|+> 2 + 2 (-|*> (*|-> 2 . (1.5) 

Therefore, when the observer has access to a particular part of the system (system 1), all 
the information obtainable through the subsystem is contained in the reduced density 
matrix (this assumes a statistical expectation value). 

Noting that the states of the system 2 are orthogonal, 2 (+|~ ) 2 = 0, the reduced 
density matrix becomes diagonal 

a x = Tr 2 {|¥> (*|} = \ (|+> (+\\ + (-\\ . (1.6) 

This result corresponds to the density matrix of a mixed-state of the system 1, i.e. in 
either one of the two states and |— ) 1 with equal probabilities as opposed to the 
superposition state \^f) . A suitable interference experiment could confirm if it is a pure 
or a mixed-state, but if the observable belongs only to system 1, the previous calculation 
demonstrates that it is impossible to distinguish between a pure or a mixed-state. We 
should not forget that this would not happen if the two elements, the system (1) and 
the environment (2) were not entangled. This demonstration could be extended to an 



arbitrary system of N elements as discussed in ref. |Sch04] . While eq. (1.6) could be 
misinterpreted as it means that the state of the system is in both states at the same time, 
it is important to remark that the density matrix is a mathematical tool to calculate the 
probability distribution of a set of outcomes of a measurement of the physical reality 
but it does not represent a specific state of the system. 

Thus, the interaction of a quantum system with an environment destroys the quan- 
tum superposition leading the system to a statistical mixture of states. This process 
called decoherence has many implications in the foundations of quantum mechanics 
like the problem of quantum measurements, the quantum to classical transition and 
irreversibility |Zur03t ISch04j . But questions arise not only at a basic level. As a real 
quantum system can never be isolated, when the technology gives us the possibility to 
work with systems where quantum phenomena appear, the understanding of decoherence 
becomes relevant to exploit the potentialities of quantum superpositions . 

In all the examples treated up to this point, the states constituting the quantum 
superposition have the same probability to exist. However, what happens when the 
probabilities are different? Moreover, what happens if the probabilities are time de- 
pendent? This leads to temporal interferences that appear in numerous experiments. 
For example, if we open each of the cat's boxes in an ensemble, at one minute or after 
one hour the probability distribution of the cat's state found is different. The same 
happens in the double slit experiment if there is an obstacle that blocks the slit oscillat- 
ing between them. The interference pattern will be different depending on the time of 
observation. What happens now if we consider the environment effects? Including the 
SE interaction the quantum evolution is more complicated. There is no simple expla- 
nation for the appearance of decoherence because as we said previously one deals with 
an environment that has many degrees of freedom. More importantly, the decoherence 
affects the phases of the quantum superposition states, whose consequences are difficult 
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Figure 1.3: Schematic representation of decoherence showed in ref. |Sch00j . Decoherence 
slides the quantum world into the classical world. The oscillations from a quantum 
superposition of a single atom, being in two places at once, gradually die out to lead the 
atom to one definite place (either position). Reprinted by permission from Macmillan 
Publishers Ltd: Nature (W. P. Schleich, Nature 403, (2000) 256), copyright (2000). 



to observe and understand. The first quantitative evaluation was given by Feynman 
and Veron (1963) [FV63J where they calculated dissipation through an environment of 
harmonic oscillators. Then, there were contributions from other people like K. Hepp 
and E.H. Lieb (1974) [HL73] and Wojciech Zurek (1981,1982) |Zur81l IZu782] who, while 
using less realistic models, suggest the universality of the effect and the relation with 
the measurement theory. However, the most complete work, in my opinion, was done 
by Caldeira and Legget (1983) |CL83c| [CL83al ICL83b] . 

One of the first techniques, if not the first, in allowing the experimental control of 
the temporal evolution of quantum states was the nuclear magnetic resonance (NMR). 
In this thesis, we consider NMR experiments in connection with the physical reality of 
the theoretical interpretations. 



1.3 NMR: The workhorse of quantum mechanics 

The origins of Nuclear Magnetic Resonance dates from 1930 when Isidor Isaac Rabi dis- 
covered a technique for measuring the magnetic characteristics of atomic nuclei. Rabi's 
technique was based on the resonance principle first described by Irish physicist Joseph 
Larmor, and it enabled more precise measurements of nuclear magnetic moments than 
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had ever been previously possible. Rabi's method was later independently improved 
upon by physicists Edward Purcell and Felix Bloch in 1945 |BHP46al IBlo46L IBHP46bl 
PTP46, PPB46]. Later on, the technique was improved by the advent of fast comput- 
ers and the development of pulse techniques that, through the Fourier transform, used 
the temporal evolution of the signal to notably optimize the acquisition time. The first 
experimental observations of the temporal evolution of a two-state system were done by 
H.C. Torrey (1949) |Tor49j and Erwin Hahn (1950) |Hah50a] where essentially a 1/2- 
spin system (two-state system) is under the presence of a static field H , which splits 
the energy levels of the states |+) and |— ) of each spin [see fig. 1.4 a)]. Then, through a 
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Figure 1.4: Oscillation between two states and the free induction decay, a) Energy 
splitting, HQ, of the states of a spin 1/2. b) Scheme of the spin precession around the 
static field, c) Experimental free induction decay obtained by E. Hahn (1950) [Hah50aJ. 
Reprinted figure with permission from E. L. Hahn, Phys. Rev. 77, 297 (1950). Copyright 
(1950) by the American Physical Society. 



transversal field with a radio-frequency (RF) pulse, one can build a superposition state 
a |+) + b | — ) whose dynamics can be interpreted as a classical precession around the 
static field direction with the Larmor frequency f2 [see fig. 1.4 b)]. Fig. 1.4| c) shows 
the original experimental data taken by Hahn [Hah50aJ, where one can observe, after 
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detection, a manifestation of the oscillation between the two states in an ensemble of 
spins. The attenuation of the oscillations is a consequence of the interaction with the 
environment, the other degrees of freedom that are not controlled and not observed. 
The simplest description of the experiment is to consider one spin and the other spins 
representing a spin-bath (the environment) whose interaction with the system (the se- 
lected spin) leads to decohere at a characteristic time T2 called the spin-spin relaxation 
time. 

From its fundamental beginnings, the NMR technique turned out soon into a pre- 
cise spectroscopy of complex molecules which triggered impressive instrumental develop- 
ments. However, nuclear spins and NMR keep providing wonderful models and continued 
inspiration for the advance of coherent control over other coupled quantum systems. It 
has gained the role of the workhorse of quantum dynamics. NMR was involved in the 
beginning of the experimental quantum information processing (QIP) applications, al- 
though nowadays, it is not considered feasible because its difficult scalability |QCR04|. 
However, in Vandersypen and Chuang words |VC04j . NMR 

"being one of the oldest areas of quantum physics[, give us the possibility 
to play with quantum mechanics because it] made possible the application 
of a menagerie of new and previously existing control techniques, such as 
simultaneous and shaped pulses, composite pulses, refocusing schemes, and 
effective Hamiltonians. These techniques allow control and compensation for 
a variety of imperfections and experimental artifacts invariably present in 
real physical systems, such as pulse imperfections, Bloch-Siegert shifts, un- 
desired multiple-spin couplings, field inhomogeneities, and imprecise system 
Hamiltonians. 

The problem of control of multiple coupled quantum systems is a signature 
topic for NMR and can be summarized as follows: given a system with Hamil- 
tonian 7i = 7i sys + 7^ contro i ; where 7i sys is the Hamiltonian in the absence 
of any active control, and H contro i describes terms that are under external 
control, how can a desired unitary transformation U be implemented, in the 
presence of imperfections, and using minimal resources? Similar to other 
scenarios in which quantum control is a welldeveloped idea, such as in laser 
excitation of chemical reactions [Walmsley and Rabitz, 2003], 7i C ontroi arises 
from precisely timed sequences of multiple pulses of electromagnetic radiation, 
applied phas e- coherently , with different pulse widths, frequencies, phases, and 
amplitudes. However, importantly, in contrast to other areas of quantum con- 
trol, in NMR 7i sys is composed from multiple distinct physical pieces, i.e., the 
individual nuclear spins, providing the tensor product Hilbert space structure 
vital to quantum computation. Furthermore, the NMR systems employed in 
quantum computation are better approximated as being closed, as opposed to 
open quantum systems. " 

Vandersypen and Chuang. 

Thus NMR inspired other techniques in the methodology of quantum control |PJT + 05j . 
In fact, the first realization of a SWAP operation in solids, an essential building block for 
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QIP, could be traced back to a pioneer NMR experiment by Miiller, Kumar, Baumann 
and Ernst (1974)0 |MKBE74j . While they did not intended it as a QIP operation, they 
described theoretically and experimentally the swapping dynamics (cross polarization) 
of two strong interacting spin systems and had to deal with the coupling to a spin-bath. 
Until that moment, all the experiments considering two interacting spins were treated 
through hydrodynamical equations |For90j using the spin-temperature hypothesis that 
leads to a simple exponential dynamics. Miiller, et al. (MKBE) showed that, in a case 
where the coupling between two spins is stronger than with the rest, one has to con- 
sider quantum coherences in the quantum calculations. They modeled the experiment 
treating quantum mechanically the two-spin system and considering the coupling with 
the spin-bath in a phenomenological way as a relaxation process. The original figure 
published in the paper is shown in fig. 1.5 where typical cross-polarization (swapping) 
dynamics for three different internal interactions (coupling between the two-spins) in 
ferrocene are displayed. One can clearly observe the frequency change of the quantum 
oscillation. More recent experiments, spanning the internal interaction strength were 
done by P. R. Levstein, G. Usaj and H. M. Pastawski |LUP98j . By using the model of 
MKBE [MKBE74], they obtained the oscillation frequency and the relaxation for dif- 
ferent interaction strengths. These results are shown in fig. L6 where one can observe 
striking changes in the relaxation time and frequency as a function of the control param- 
eter. Since this discontinuous change is not predicted by the standard model of MKBE, 
it remained unexplained. The description and interpretation of this striking behavior 
are among the main results of this thesis. 

Thus, in view of possible applications to fields like quantum information processing 
|Kan98l IBDOO] , the experimental manifestation of these dynamical interference phenom- 
ena in qubit clusters of intermediate size has a great interest. However, experimental 
realizations and control of a pure-state dynamics is still one of the challenges in nowadays 
physics |QCR04|. Therefore, one generally has to deal with an ensemble evolution, which 



is the case of the states involved in NMR, i.e. the dynamics of an initial mixed-state. 
One can generate mixed-states that are called pseudo-pure because they are constituted 
by a pure-state plus a mixed-state density operator. Numerous spin dynamics NMR 
experiments have shown surprising quantum phenomena |PLU95l lMBSH+971 IRSB+05] . 
The difficulty to produce pure-states in a high temperature sample leads to the develop- 
ment of the ensemble quantum computation [VSUQlJ ISSB05] . However, as we mention 
previously if the system is too complex, it is hard to mathematically describe its tem- 
poral evolution. This is a consequence of the exponential growing of the Hilbert space 
dimension as a function of the number of elements in the system. In order to overcome 
this limitation, we take profit of the quantum parallelism |SKL02j and the fragility of 
the quantum superpositions to develop a method that describes ensemble dynamics. 

As the dimension of the system increases, the sensitivity of the quantum super- 
position might lead to the inference that quantum phenomena will not manifest at 
macroscopic scales [M KT + 0Q[ ISchOOj . In contrast, an experimental demonstration of 
macroscopic quantum states done by Y. Nakamura, et al. j.Xl'T!)!). Aye99] shows the 



2 A similar work where transient oscillation where observed was presented the next year by D. E. 
Denco, J. Tegenfeldt and J. S. Waugh |DTW75j . 
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Figure 1.5: Transient oscillations in a cross-polarization experiment by Miiller, Kumar, 
Baumann and Ernst (1974) [MKBE74] . The two -spin dynamics coupled to a spin-bath 
is shown for three different internal couplings. Reprinted figure with permission from L. 
Miiller, A. Kumar, T. Baumann and R. R. Ernst, Phys. Rev. Lett. 32, 1402 (1974). 
Copyright (1974) by the American Physical Society. 



opposite. Indeed, there is no doubt about the high sensitivity of the quantum superposi- 
tion states in large systems which paves the way for an effective decoherence when there 
are interactions with the environment. As any environment usually has many degrees 
of freedom, it is very difficult to reverse the SE interaction constituting the dominant 
source of irreversibility in nature |Zur03l ISch04] . Numerous works are related to this 
topic, but we should begin discussing the pioneer work that made a temporal reversion 
of a quantum dynamics: the Hahn's echo experiment. It is based on the reversion of 
the dephasing among rotating spins due to inhomogenities of the static field [Hah50bJ. 
He observed an echo in the NMR polarization signal (see fig. 1.7) manifesting the de- 
terministic nature of quantum mechanics, but with an attenuation rate proportional to 
the spin-spin coupling. The forward dynamics is a consequence of the interaction of the 
spins with the static field and the spin-spin interactions, but only the interactions with 
the static field are reverted. Thus, the dipolar interaction remains working. Within 
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Figure 1.6: Fitting parameters of a two-spin cross-polarization experiment performed 
by P. R. Levstein, G. Usaj and H. M. Pastawski (1998) |LUP98] . Striking changes in 
the relaxation and in the oscillation frequency behaviour are observed. These effect are 
described in chapter |4} Reprinted with permission from P. Levstein, G. Usaj and H. 
M. Pastawski, J. Chem. Phys. Vol. 108, Issue 7, Page 2718, 1998. Copyright 1998, 
American Institute of Physics. 



the NMR field, there were many experiments using the deterministic nature of quantum 
mechanics to take out some interactions that disturb the relevant system evolution. But, 
the first work that emphasizes the deterministic nature of quantum mechanics, invali- 
dating the spin temperature hypothesis (thermodynamical approaches), was done by W. 
-K. Rhim and A. Pines and J. S. Waugh [RPW70] . They called a "Loschmidt daemon" 
to the process of reversion of the dipolar interaction in the "magic echoes" experiment. 
There, they observed an echo signal after undoing (reversion control) the evolution under 
spin-spin interactions that remain untouched in the Hahn's echo experiment. The pre- 
vious experiments evolve from multi-spin initial excitations. The local initial excitation 
version of the "magic echoes" was done by S. Zhang, B. H. Meier and R. R. Ernst (1992) 
[ZME92bJ. They called this experiment as "the polarization echo" where they used a 
very ingenious idea to observe a local magnetization X.\lK92bl IZME92a] . They used 
a rare nucleus, 13 C, bonded to a 1 H nucleus (abundant) as a local probe to create and 
observe the local polarization. However, we have to remark that while one increases the 
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Figure 1.7: First experimental time reversion in NMR: The Hahn's echo [Hah50bJ. 
Reprinted figure with permission from E. L. Hahn, Phys. Rev. 80, 580 (1950). Copy- 
right (1950) by the American Physical Society. 



quantum control on the Hamiltonians, a minimal decay of the echoes can not be avoided. 
Experiments performed in Cordoba suggest that the quantum states are so sensitive to 
perturbations that even a very small uncontrolled perturbation generates an intrinsic 
irreversibility characterized by the own system dynamics |LUP98l IUPL981 iPLU+OOj . 
By considering an analogy with the behavior of a simpler one body chaotic system, 
this was interpreted |JP01t IJSBOlt ICPJ04] as the onset of a Lyapunov phase, where 
I/t^ = min [1/tse, A] is controlled by the system's own complexity A. However, a theo- 
retical answer for many-body systems that do not have a classical analogue characterized 
by Lyapunov exponent remains open. This is also a topic that enters in this thesis' mo- 
tivation: the improvement of our comprehension and control of decoherence processes 
and irreversibility. The complexity of many-body systems leads us to study the forward 
dynamics of open systems to characterize the decoherence process before studying the 
time reversal. 
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1.4 Our contribution 

In this thesis, we solve the dynamics of many-spin systems interacting with a spin-bath 
through the generalized Liouville-von Neumann quantum master equation beyond the 
standard approximation. Further consideration of the explicit dynamics of the bath 
helps us to solve the spin dynamics within the Keldysh formalism, where the interac- 
tion with the bath is taken into account in a precisely perturbative method based on 
Feynman diagrams. Both methods lead to identical solutions and together gave us the 
possibility to obtain numerous physical interpretations contrasting with NMR experi- 
ments. We used these solutions in conjunction with experimental data to design new 
protocols for molecular characterization, develop new numerical methodologies and con- 
trol the quantum dynamics for experimental implementations. But, most important, 
these developments contributed to improve the fundamental physical interpretations of 
the dynamics in a quantum open system under the presence of an environment. In 
particular, we show a manifestation of an environmentally induced quantum dynamical 
phase transition. 

1.4.1 Organization of this thesis 

In Chapter 2 we use the standard formalism of density matrix to solve the spin dy- 
namics using the generalized Liouville-von Neumann quantum master equation. In the 
first part of the chapter, the spin dynamics of a two-spin system coupled with a fast 
fluctuating spin-bath is solved. This system describes the cross-polarization experiment 
of MKBE [MKBE74J. We start using the standard approximations and then we extend 
the solution without these restrictions. We compare the solutions and remark the main 
differences. We analyze the spin dynamics for different anisotropics of the SE interac- 
tions given by the different contributions of the Ising and the XY interaction. We show 
how the rates of decoherence and dissipation change depending on the anisotropy ratio 
between the Ising and XY coupling. In the second part of the chapter, we extend the 
solution to a three-spin system coupled with a spin-bath. The solutions obtained are 
applied to experimental data to get more detailed information for molecular character- 
ization. In particular, we use the three-spin solution to characterize the liquid crystal 
8CB and incorporating some memory effects, we conclude that the spin-bath has a slow 
dynamics. 

In Chapter 3 we solve the spin dynamics within the Keldysh formalism |Kel64] . The 
Keldysh formalism is well established in the electron transport description. Through the 
Jordan- Wigner transformation |JW28] . we map the two-spin system of chapter 2 into 
a fermion system. We find how to describe the SE interaction within the wide band 
approximation (fast fluctuation inside the bath) and we obtain a solution for the spin 
dynamics that improves the standard solution of the generalized Liouville-von Neumann 
quantum master equation. Here, we use a microscopic model to obtain the spin dynam- 
ics that avoids using a phenomenological description of the SE interaction. However, we 
obtain the same solution going beyond the standard approximation within the density 
matrix formalism. Then, we solve the spin dynamics of a linear chain including all the 
degrees of freedom of the environment in the calculations and we show how the memory 
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effects induce a time dependence in the oscillation frequency as is observed experimen- 
tally. We develop a stroboscopic model to describe decoherence which is optimized for 
numerical applications. This model converges to the continuous expression. 

In Chapter 4 based on the solutions obtained in previous chapters we describe a 
manifestation of an environmentally induced quantum dynamical phase transition. We 
show the experimental evidence and interpret the phenomenon in detail. In particular, 
we show how the anisotropy of the SE interaction has an important role in the critical 
point of the phase transition. An extension of this phenomenon to a three-spin system 
shows how to vary the control parameter to "isolate" two of them from the environment. 

In Chapter 5, inspired in the stroboscopic model developed in chapter 3, we propose 
a new NMR pulse sequence to improve the transfer of polarization through a specific 
pathway in a system of many interacting spins. The sequence effectively prunes branches 
of spins, where no polarization is required, during the polarization transfer procedure. 
Simulations of the spin dynamics in the 13 C backbone of leucine are performed. Possi- 
ble applications and potential fundamental contributions to engineered decoherence are 
discussed. 

In Chapter 6 we develop a novel numerical method to obtain the spin dynamics of 
an ensemble. It overcomes the limitations of standard numerical calculations for large 
number of spins because it does not involve ensemble averaging. We exploit quantum 
parallelism [SKL02J and the fragility of a randomly correlated entangled state to repro- 
duce an ensemble dynamics. 

In the final part of each chapter a brief summary of the main original contributions 
including references to publications is included. 

In Chapter 7 we summarize the whole work emphasizing the main conclusions and 
perspectives. 
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Chapter 2 



Many-spin quantum dynamics 
within the density matrix formalism 



The exact quantum dynamics of small quantum systems has regained interest during 
the last years [ALW92J, due to the technological advances that give us the opportu- 
nity to observe quantum phenomena. Spin systems are good candidates in this respect 
and provide beautiful playgrounds for fundamental studies. Besides, several challenging 
applications require a very fine knowledge of the spin interactions, such as molecular 
characterization, spin control in nanodevices |SKE + 01[ IKLG0 2] and quantum compu- 



tation |GC97} ICPH 98, BDOOj. In the introduction became evident the limitations of 
simple thermodynamical arguments |For90j based on the spin temperature hypothesis. 
The experiment of MKBE |MKBE74] showed the need to consider the system quan- 
tum mechanically keeping the quantum coherences to describe the transient oscillations. 
However, the first work that showed the weakness of the "spin temperature" hypothe- 
sis was done in 1970 |RPW70] . In it, a time reversal of the spin-spin interactions was 
performed. It was followed by numerous nuclear magnetic resonance (NMR) experi- 
ments that have demonstrated the time reversibility of the dipolar (many-spin) evolu- 
tion jZME92bl IEMTP98a[ IEMTP98bl ILUP981 IUPL98] leading to revise the concept of 



"spin diffusion" |PLU95l IPUL961 IMBSH+971 IWau98j . More importantly, by selecting 



appropriate systems and pulse sequences, one can investigate the sources of quantum 
decoherence |Zur03l ISch04] , ergodicity |PLU95t IPUL96t IWau98] , and quasi-equilibrium 
|SHE98j . 

From a practical point of view, spin dynamics observed by NMR has proved very 
powerful in order to characterize molecular structures and dynamics jSRS96j . Exper- 
imental observations together with simple analytical solutions for few-spin dynamics 
can provide detailed information on the intra and intermolecular interactions [MKBE74, 
ILUP98t IUPL98] . This is particularly important for the characterization of complex fluids 
in their native state, where one uses cross-polarization (CP) dynamics [HH62, Sli92] to 
evaluate order parameters [PR96] . However, the reliability of these and other structural 
and dynamical parameters depends on the accuracy of the spin dynamics description to 
which the experimental data are fitted. 

In this chapter, we use the standard formalism of density matrix to solve the spin dy- 
namics using the generalized Liouville-von Neumann quantum master equation [Abr6I4 
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EBW91J. In the first part of the chapter, we solve the spin dynamics of a two-spin sys- 
tem coupled to a fast fluctuating spin-bath. This system describes the cross-polarization 
experiment of MKBE |MKBE74] . As a first step, we use the standard approximations 
and then we extend the solution releasing these restrictions. We compare the solutions 
and remark the main differences. We analyze the spin dynamics for different SE interac- 
tions consisting of different Ising and XY contributions. We show how the decoherence 
and dissipation rates change depending on the anisotropy ratio between the Ising and 
XY couplings. In the second part of the chapter, we extend the solutions to a three-spin 
system coupled to a spin-bath. The solutions are applied to get more detailed informa- 
tion from our NMR experimental data. This leads to new methodologies for molecular 
characterization. In particular, we use the three-spin solution to characterize the liquid 
crystal 8CB. The slow dynamics of the smectic phase, experimentally observed, lead us 
to include some spin-bath memory effects. 

2.1 Quantum dynamics of a two-spin system 

For didactical reasons, we start solving the spin dynamics of an isolated two-spin system. 
Then, we will include the interactions with the spin-bath. 

2.1.1 Quantum evolution of an isolated two-spin system 

We solve the evolution of an isolated two-spin system during cross-polarization (CP). 

In this procedure, two different species of spins, S-I, which here will correspond to a 
13 C- 1 H system are coupled in such a way that they "believe" that they are of the same 
species |Abr61l fSli92t IEBW91J . In that situation, the most efficient polarization transfer 
can occur. The system Hamiltonian, in presence of a static field H and the radio 
frequency fields of amplitudes H\ t s and H\j with frequencies uj^s and u> T f y i respectively, 
is given by |Abr61l ISE92] 

n s = -hn , s s z - nn 0:I i z 

— hQi jS js™ cos (u^s t) + S y sen (uj t{)S t)\ — HQij (l x cos (k>rf,s t) + P sen (c<j rfj s £) j 

+ 2bFS z , (2.1) 

where 

n ,i = jiH ,i,i = S,I (2.2) 
are the precession Larmor frequencies in the static field and 

£l lii = 'y i H 1 , i ,i = S,I (2.3) 

are the Zeeman (nutation) frequencies of the RF fields. The last term is the truncated 
dipolar interaction assuming that 

|MV - f^oM > |6| • (2-4) 
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The amplitude of the interaction is |Sli92] 

6 = -i(^)(3^-l), (2.5) 

where r is the internuclear distance and 9 the angle between the static field and the 
internuclear vector. In the double rotating frame [Sli92], at the frequencies of the RF 
fields, the system Hamiltonian becomes 

H s =hAn s s z + hA^i* - hn hS s x - htttjp + 2bs z i z (2.6) 

where 

AQi = Q 0ti - u a>i (2.7) 
are the respective off-resonance shifts with i = I, S . We assume the conditions 

Aft/ = AQ S = (2.8) 

and 

H\n ljS + ^i,i\ > |6| (2.9) 

which are obtained when the RF fields are applied on-resonance and when the RF power 
is much bigger than the dipolar interaction. Thus, the doubly truncated Hamiltonian 
becomes 

H B =\ (S - A) S x + |(E + A) r + b [s z i z + S y I y ^j (2.10) 
= |s (s x + r^j + \a (i x -s^j+b (s z i z + S y F^ , (2.11) 

where the non-secular elements of the dipolar interaction with respect to the S (S x + I x ) 
term have been neglected. Here, S and A are defined by 

E = —h(fli,s + ^1,/) and A = h (Qi, s - ^1,1) . (2.12) 

Within the Hartmann-Hahn condition |HH62j ISli92] . A = 0, the two spins act as they 
belong to the same species improving the polarization transfer between them. To obtain 
the quantum evolution, we solve the Liouville-von Neumann equation |Abr61[ IEBW9lj 

^PW = 4[^s,p(t)], (2.13) 
where p is the density matrix operator of the two-spin system. Its solution is given by 

p(t) = U{t)mU-\t), (2.14) 

with 

= exp f-i# a ^ (2.15) 

the evolution operator and p(0) the initial condition. For the last, we consider the 
: H totally polarized and the 13 C depolarized. This can be experimentally achieved by 
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rotating, with a 7r/2 pulse, the equilibrium polarization of the 1 H with the Zeeman field 
H Q in the z direction to the XY plane assisted by a cycling pulse sequence |LUP98j . The 
initial condition (immediately after the ir/2 pulse) is expressed as 

where 

Z = tr<fexp I -^4^ I } (2.17) 



is the partition function. 

In the high temperature approximation 

, . i + 8 B m ii x 

p o = , ^ ? , 2.18 

Tr{l} 

where /?b = l/ksT. 

Thus, calculating the evolution operator, we obtain the magnetization of the 13 C as 
a function of the contact time t of the cross-polarization 

M s ,( ( ) = T i {p W S-}=M !l^^l, (2.19) 

and for the 1 H we obtain 

M„ (t) = Tr {p (i) f} = M ' 1 + C ° 2 S '"° t " . (2.20) 

Here, 

oo = b/h (2.21) 

is the Rabi frequency and 

M = ^-/feMV (2-22) 



is the initial magnetization at the 1 H. Figure 2.1 shows these two curves as a function 



of t. The Hamiltonian (2.10) has only Zeeman fields along the x direction. Thus, by 



changing the axis names: x — *> z, y — *> x and z — > y the Hamiltonian becomes 

H s =l (E - A) S* + \ (E + A) I* + b (s x I x + S y lA (2.23) 
= |E (§ z + I*\ + |A - 5") + §6 + S~I + ] . (2.24) 

Now, it is evident that the dipolar interaction is an XY (flip-flop) term that splits the 
energy level of the states |T, I) and ||, t) an d induces an oscillation between them. This 



is manifested in the oscillation of figure 2.1 where the magnetization is totally transferred 



forth and back from the 1 H to the 13 C with the Rabi (swapping) frequency u>o. Within 
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4-* 




Figure 2.1: Polarization evolution of the S spin (black line) and the / spin (red line). 
The oscillation frequency is Uq = b/h and the maximum amplitude is given by the initial 
polarization of the spin /. 



this representation^ in the new basis of I z , S z given by {|+, +), 
where |±, ±) = |±) <g> |±) = \M I>Z ) <g> \M S>X ) with P\M I)Z ) = M I: 
magnetization of the 13 C spin can be expressed as 



1+,-), I- +>, I-,")}, 
,\M ItZ ), M I)Z = ±±,the 



M s * (t) = Tr p(t) S 



(t) ~ P—,— (t) + P-+.-+ (t)) - P+-.+- (*) , (2.25) 



with p±±,±± (t) = (±, ±|p(t)|±, ±). Within this basis, the diagonal terms (populations) 
of the density matrix contribute positively to the magnetization when the carbon is in 
the state up and negatively when it is down. This means that the carbon magnetization 
is given by the difference between the populations of the state up and down. We can see 

that the elements of the density matrix p ++i++ (t) and p i (t) are constants of motion 

and give the first term of eq. (2.19). The difference between p + _ >+ _ (t) and p_ + _ + (t) 
(coherences) gives the oscillatory term, which describes the transition between |f, j) and 



Remember that now the z direction is the originally x direction. 
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t) . The Tr |p (t) S* 2 j is given by 



M 5 , (t) = l - [pn it) - p 44 (t) + p 32 (t) + P23 (*)] = \ [pn it) - P44 (t) + 2 Re {p 32 (*)}] 



in the basis of the eigenstates of the Hamiltonian 

|1> = !+,+>, |2> 
The states |1) and |4), previously defined as |+, +) and 



(2.26) 



^(l+,-) + |-,+)),|3) = i=(|+,-) 



- +»,|4> = |- -)|- 
(2.27) 

-), give again the first term 
of eq. (2.19 ). However, the oscillatory term is given by the real part of the P32 coherence. 



The magnetization on the I spin has the sign of the oscillatory term changed. Thus, the 
total magnetization of the system is a constant of motion described by 



M tot . (t) = Tr \p(t) [S z + i 



P++.++ (*) - P— ,— (*) = Pn (t) ~ P44 (*) = M , 

(2.28) 



where Mo is the initial magnetization. The blue line in fig. 2.1 represents Mo/2 which 
constitutes the mean magnetization of each species. 



2.1.2 A two-spin system interacting with a spin-bath 

We use the model proposed by Miiller et al. |MKBE74] to describe the experimental 
spin dynamics of a molecule of ferrocene and to characterize the quantum dynamics of a 
two-spin system interacting with a spin-bath. The model assumes that only one spin, Ji, 
interacts with the spin-bath which is described in a phenomenological way. The modeled 
Hamiltonian is 



Ti-SE 

He 



m z (s z + A 2 ) + \b [s + i{ + s~it 

al{F z + \(3 (i+F- + /f/F 



j>i 



2I Z T Z — i I T + T~ 4- T~T + 



(2.29) 
(2.30) 

(2.31) 

(2.32) 



where Mli is the Zeeman interaction and we are assuming the Hartmann-Hahn condition 
[HH621 ISlr92] 

Q 1>s = = n z . (2.33) 

7is is the system Hamiltonian of the two coupled spins, He is the spin-bath Hamiltonian 
with a truncated dipolar interaction and Use is the system-environment (SE) interaction 
with 

F u = J2bJr, u = x } y,z (2.34) 

i 
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and 



F x ± iF y 



(2.35) 



Hse is an Ising interaction if /3/ot = and an XY, isotropic (Heisenberg) or the truncated 
dipolar interaction if a/ (3 = 0, 1, — 2 respectively. This last case is typical in solid-state 
NMR experiments |Abr61t ISli92l IEBW9lj . In a quantum mechanical relaxation theory 
the terms F u are bath operators while in the semi-classical theory |Abr61l [EBW91J 
the F u (t) represent classical stochastic forces. The experimental conditions justify a 
high temperature approximation, and hence the semiclassical theory coincides with a 
quantum treatment |Abr61j . By tracing on the bath variables, the random SE interaction 
Hamiltonian is written as 



Use (t) = aF z (t) I» + \(5 \F~ (t) 7+ + F + (t) If 
The time average of these random processes satisfies 

F u (t) = 0, 

where their correlation functions are 

g {u ' v) (r) = F u (t) F v * (t + r). 



(2.36) 



(2.37) 



(2.38) 



Following the usual treatment to second order approximation, the dynamics of the 
reduced density operator is given by the generalized Liouville-von Neumann differential 
equation |Abr61| IEBW91] 



dt 



a it) = ~[H a , cr (t)] - if {& (t) - a^} 



where a (t) is the reduced density operator 

a(t) = Tr E {p(t)} 



(2.39) 



(2.40) 



with Trg denoting a partial trace over the environment variables. The relaxation super- 
operator T is given by the SE interaction. It accounts for the dissipative interactions 
between the reduced spin system and the spin-bath and it imposes the relaxation of the 
density operator towards its equilibrium value a^. 

We assume that the correlation times of the fluctuations are extremely short com- 
pared with all the relevant transition rates between eigenstates of the Hamiltonian, i.e. 
frequencies of the order of Qz and u = b/h. In this extreme narrowing regime or fast 
fluctuation approximation we obtain 



fH = ^^,J M (o) fes [if,* 



where 



J 



dr 
~h 



g M (r)exp{-iwr} 



(2.41) 
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is the spectral density and 

= [aS UjZ + (3 (8 U>X + S UtV )] [a5 VyZ + (3 (5 V)X + 8 VtV )\ . (2.42) 

Assuming that the spatial directions are statistically independent, i.e. g( u ^ (r) = if 
u 7^ v, the superoperator T can be written as 



T{a} = a 2 J z I{, I*, a 



(3 2 J X 



TX TX ~ 



+ (3 2 J V 



T v T y a 



where 



1 



(2.43) 



(2.44) 



Notice that the axial symmetry of Us around the z axis leads to the impossibility to 
evaluate separately J x and J y , so they will appear only as the averaged value 



jxy = ^(J X + J y ) 



(2.45) 



Thus, we obtain 



T{a} 



TIL 



jz jz ±. 



r 



XY 



I x , If, a + If, If, a 



where 



Tzz = oi 2 J z and T XY = [3 2 J xy . 



(2.46) 



(2.47) 



Note that T zz and T X y contain the different sources of anisotropy. The usual ap- 
proximation considers J x = J y = J z (identical correlations in all the spatial di- 
rections) and a — f3 — 1 (isotropic interaction Hamiltonian) [MKBE74J. A bet- 
ter appro x imation considers a dipolar interaction Hamiltonian, i.e. a = —2(3 = —2 
|CAL + 03 , ADLP06 , ALP07] . This is in excellent agreement with previous experi- 
ments in polycrystalline samples where fittings to phenomenological equations have 
been performed |LUP98| IRHGG97] . In particular, in the case of isotactic polypropy- 
lene [RHGG97J, a fitting where Rd p corresponds to a 2 J z /h and i?df = [3 2 J xy /h, gives 
Rdp/Rdf ~ 4. 

We consider the experimental initial local polarization (2.18) on the spin Ii, 



&(0) 



i + fennel? 

Tr{i} 



(2.48) 



and the spin-bath polarized, where /3b = 1/ (k B T). By neglecting other relaxation 
processes T\ p , etc.), the final state reaches the temperature of the spin-bath yielding 



Or. 



i + /3 B hn 0!l ^s z + it 
t7{i} 



(2.49) 



Here, commutes with Hs, not containing coherences with AM > 1. 
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2.1.2.1 Neglecting non-secular terms in the relaxation superoperator 

Following the standard formalism |Abr61l TEBW91j . we write the superoperator T using 



the basis of eigenstates of the system Hamiltonian (2.30). As the SE Hamiltonian is 
secular with respect to the RF fields a block structure results. If M is the total spin 
projection in the z direction, the first block couples the populations and off-diagonal 
elements with AM = 0, zero quantum transitions (ZQT), of the density matrix. Each 
of the following blocks couples one order \AM\ > 1 of the off-diagonal elements of the 
density matrix among themselves. As the initial and final conditions do not contain 
coherences with \AM\ > 1, the evolution of the density operator is reduced to a Li- 
ouville space restricted to populations and ZQT. Thus, in the Hamiltonian eigenstate 



basis (2.27), the generalized Liouville-von Neumann quantum master equation (2.39) 



restricted to this subspace becomes 

- =1 

dt a | h'"° h 



A <t={~JU-\t\{<t-g 00 ), (2.50) 
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(2.51) 



Here, the superoperator T is given by 
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and the superoperator of the system Hamiltonian, that is diagonal in its basis, results 
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-b 



(2.53) 



where its elements are transition frequencies (energy differences). After neglecting the 
rapidly oscillating non-secular terms with respect to the Hamiltonian, i.e. 



zz, J- XY 



-C 6, 



(2.54) 



a kite structure results |EBW91] . All the non-diagonal terms coupling the population 
block with the ZQT block are non-secular and can be neglected because the Hamiltonian 
(2.30) does not have degenerate eigenenergies. Although in this case the ZQT block 
is diagonal, in a general case only the diagonal terms of this block contribute to the 
evolution if there are not degenerate transitions. The differential eq. (2.51) is now 



d 
dt 
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(2.55) 
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where its solution, in the Hamiltonian eigenstate basis (2.27), is 



a(t) 



| + Mo(l-^f^ 







-Mpe-^e 



(Zzz+?£xy), 








( r zz+ 2r xy) , 



2 










M (l 



e -r XY t/fi 



At time t = 0, we obtain the initial state of the reduced density matrix (2.48) and when 
t — > oo, it goes to the final state ff m of eq. (2.49). The elements an (t) and cr 44 (t) , 



which are the populations of the states |+, +) and | — , — ) respectively, go to the equilib- 
rium state (Too with a rate Fxy/h- This accounts for the net transfer of magnetization 
from the spin-bath to the system and it contains the information of the net magnetiza- 
tion inside the system. The coherences 023 and 032, that take into account the swapping 
between the states |j, f) and |f, |) with the natural frequency u , decay to zero with a 
decoherence rate 

(r zz + 2r XY ) 



2h 



(2.57) 



Note that the coherences decay faster than the time that the system takes to arrive 
to the eq uilibrium state. We calculate the magnetization of the spin 5* obtaining an 
extension |ALP07] of the result given in ref. |MKBE7"4] 



M s * (t) = Tr { a (t)S 



1 _ ± e -T X rt/h _ 



1 



.i(r zz +2r XY )t/h 



cos (coot) 



(2.56 



where our essential contribution is that we specifically account for the anisotropy arising 
from the nature of the SE interaction reflected in Tzz = o?J z and Txy — (3 2 J xy . The 



first two terms of eq. (2.58) are given by [an (t) — cr 44 (£)] /2 and the oscillatory one by 
Re {(T32 (t)} [see eq. (2.26)]. The sum of the first two terms are the mean magnetization 
at each site or, multiplied by two, it represents the total magnetization of the 2-spin 
system which is given by 



Aftot. (t) =Tr\&{t)[S K + in\ = a n (t) - a M (t) = 2M 1 



rt/H 



(2.59) 



The time dependence of this quantity is due to a "diffusion process" from the spin-bath 
that injects magnetization at a rate Txy/^ through the XY SE interaction term. We 
define the SE interaction rate as 



Tse — Tzz + Txy 
and the weight of the XY interaction as 



(2.60) 



(2- 



Pxy — Txy/Tse- 



(2.61) 
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An Ising, dipolar, isotropic, and XY SE interactions are obtained when pxv = 0, |, |, 1 
respectively. Equation (2.58) becomes 



M s * (t) = M 



I _ I e -Px Y r SE t/?i _ l e -f (i+PxY)r SE iA 



COS ((dot) 



(2.62) 



Figure 2.2 shows typical curves for different SE interactions (black lines). The blue 




Figure 2.2: Temporal evolution of the polarization in the two-spin system interacting 
with a spin-bath for different system-environment interactions. Panels a), b) and c) 
show the polarization curves for the isotropic, XY and Ising SE Interaction, where pxv 
takes the values 1/2, 1 and respectively. The black line is the magnetization of the 
spin S and the red line is the magnetization of I\ for a ratio Tse/ (<^o^) = 0.1. The blue 
line is the temporal evolution of the mean magnetization at each site. The green lines 
show the coherence decays with respect to the mean magnetization evolution. Panel c) 
shows the evolution of the S magnetization with a dipolar SE interaction for different 
values of the ratio Tse/ (^o^)- The orange, black and magenta lines show the temporal 
evolution for the ratios 1, 0.1 and 0.01 respectively. 



line is the temporal evolution of the total magnetization divided by two or equivalently 
the mean magnetization at each site. We see that the curves go to Mq for long times 
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manifesting that the system arrives to the equilibrium state with the exception of 
the Ising SE interaction. In this case, the system has no injection from the spin-bath and 
goes to the quasi-equilibrium state of the 2-spin system, i.e. the initial magnetization is 
spread over both sites. This quasi-equilibrium is described by 

~ qe = exp{-7V(fc B T qc )} 

Tr{exp{-ft s /(£;BT qc )}} 

with 

T qc = 2{n 1J /n ,i)T. (2.64) 

In general, there is a competition between the Ising and the XY SE interaction terms that 
leads the system t o a 2-spin quasi-eq uilibrium state or to the total system equilibrium 
state respectively CAL+03 ALP07j . This quasi-equilibrium is time dependent and is 



given by the mean magnetization, M tot (t) /2, represented by the blue line in fig. 2.2 The 
green lines in the figure show the coherence decay relative to the mean magnetization at 
each site. We see that the XY interaction is the most coherent because its decoherence 
rate is equal to the magnetization transfer rate, while in the other cases, decoherence is 
faster than magnetization transfer. The red lines show the magnetization on the spin Ii 
described by the expression 



M /z (t) = M 



1 - I e -RxvW» + i e -|( 1 +PxY)rs Ei /^ cos ^ ot ) 



(2.65) 



where only the sign of the oscillatory term Re {032 (t)} changes. Figure 2.2 c) shows 
curves with dipolar SE interaction for different values of the ratio T^e/ (ujqK) . It shows 
how the decoherence and magnetization transfer are stronger as Tse/ (^0^) becomes 
higher. Here, we observe the decoherence's role described in the Introduction. The 
temporal interference pattern is described by the oscillatory term which contains the 
entangled two-spin superposition. A strong SE interaction leads to an efficient degra- 
dation of the two-spin quantum entanglement. This drives the system to a mixed-state, 
described by the diagonal elements of the density matrix, which constitutes the quasi- 
equilibrium state represented by the blue line. When decoherence is not too strong, 
we observe that it is not necessary to wait long times to obtain the maximum magne- 
tization at the spin 5* (totally polarized). It is enough to wait for a maximum of the 
oscillation at time k/ujq where the magnetization reaches a value close to the maximum 
obtainable (M ). But a more important result is that for an XY SE interaction, one can 
achieve the biggest gain of polarization at the first maximum of the oscillation. This is 
a consequence of the different behavior of the decoherent processes arising on the Ising 



or XY interactions. Moreover, for an XY SE interaction, expression (2.58) yields all the 
maxima of the oscillation equal to Mo, regardless of the magnitude of the SE interaction. 
However, we should not forget that this expression is valid only for r zz , T X y *C b. 

2.1.2.2 Non-secular solution 

If we release the condition Tzz, Txy <C b, i.e. we do not neglect the non-secular terms for 
the superoperator T, the dynamics still occurs in the Liouville space of the populations 
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and ZQT. The solution of the generalized quantum master equation is now, 

M 5 , (t) = M (l - a e' Rot - a x cos [(u + ir])t + O ] e~ Rlt ) , (2.66) 

where the real functions u>, Ro, R\ and r\ as well as a , a\ and 0o depend exclusively 
on b, 1/tse = 2 (r zz + r XY ) /h and p XY = r XY / (Tzz + r XY ). This expression will be 
discussed in chapter [3j where it is obtained from a microscopic derivation. However, it 
is important to remark that the short time evolution, t <C tse, of the secular expression 
does not satisfy the correct quadratic quantum behavior while the non-secular expres- 
sion does. The relevance of this inertial property reflected in the quadratic short time 
evolution will become evident in chapter |4| We will see how it leads to the manifesta- 
tion of so mething that we calle d an environmentally induced quantum dynamical phase 



transition ADLP06 DALP07 . 



2.2 Three-spin quantum dynamics 

In this section, we analyze theoretically and experimentally the quantum dynamics 
of a three-spin system coupled to a spin-bath during cross-polarization (CP) [HH62, 
MKBE74J. Our analysis takes into account a pure Hamiltonian behavior for a carbon 
13 C coupled to two protons 1 H, while the coupling to a spin-bath is treated in the fast 
fluctuation approximation. This model is inspired and then applied to the methylene 
and biphenyl groups of the smectic and nematic phases of the liquid crystal 4-n-octyl-4'- 
cyanobiphenyl (8CB). We make use of the Hartmann-Hahn CP technique as a function 
of contact time to measure 1 H- 13 C and 1 H- 1 H effective dipolar interactions. This tech- 
nique has proved very useful in order to evaluate order parameters in liquid crystals 
[PR96j . Most of the previous works where transient oscillations were observed during 
CP were analyzed in terms of a single 1 H- 13 C interaction incorporating the interaction 
with other protons as a thermal bath or reservoir in a phenomenological way. However, 
many liquid crystals have alkyl chains and aromatic groups in their structures, where 
the carbon is coupled to more than one proton and the carbon-proton and proton-proton 
dipolar interactions are of the same order of magnitude. This led us to consider a set 
of three strongly dipolar coupled spins 1/2 as the main system, which in turns interacts 
with the protons of the bath. Combining detailed calculations of the three-spin dynam- 
ics with structural information which provide the relative sign of the 1 H- 13 C couplings, 
we are able to obtain separately the 1 H- 13 C and 1 H- 1 H effective interactions. In order 
to test the suitability of the formula obtained, we compare the values of the 1 H- 13 C 
couplings obtained by two procedures. One involves fitting of the data from a stan- 
dard CP experiment to the calculated dynamics while in the other the 1 H- 13 C couplings 
are obtained directly from a CP under Lee-Goldburg conditions, i.e. when the dipolar 
proton-proton interactions have been cancelled out. The advantages and disadvantages 
of each procedure are discussed. 

An interesting aspect we could observe during the CP dynamics in 8CB is that the 
rate of attenuation of the oscillations (representing the coherences) is much faster than 
that of the polarization transfer from the bath in a factor several times larger than the 
one calculated assuming isotropic interaction with the bath (MKBE74]. We analyze here, 
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the origin of this highly anisotropic behavior, not observed in solid molecular crystals 
|LUP98UHH 94|. A well differentiated relaxation behavior among the two phases seems to 
indicate that while the extreme narrowing approximation is appropriate for the nematic 
phase, the description of the smectic phase requires the consideration of the slow motion 
limit. 



2.2.1 An isolated three-spin system 

We will consider the quantum evolution of a system of three spins 1/2 coupled through 
the magnetic dipolar interaction during the contact time in a cross-polarization exper- 
iment |HH62i [MKBE74J. The system is constituted by one spin S and two spins I 
representing a carbon-13 and two protons, respectively, under the presence of a static 
magnetic field Hq in the z direction and radio-frequency (RF) magnetic fields H\ i and 
Hi t s in the x direction. The Hamiltonian including the dipolar interactions truncated 
with respect to the Zeeman field H and in a double rotating frame [Sli92] can be written 
as 

H s = hAVt^r + hAQ s s z - hQ hI i x - hn hS s x + h kS z it + d(2ifr 2 - ifi* - if if), 



fc=l,2 



(2.67) 



where, as in the previous section, 

AQj = Q j — u r fj and AQ S = Q 0tS — (2.68) 
are the resonance offsets, 



/« = 1% + II (2.69) 

with u = x,y, z, 

n 1J = ll H 1: j and n liS = ls H 1:S (2.70) 
where 7/, 75 are the gyromagnetic factors of the / and S spins. The constant, as defined 



in eq. (2.5 ) 



hk = - l - ^ sh2 / (3cos 2 (M-l) ) fc = lj2 {m 

2 47T \ f c 



Sk 

and 

1^07? ft 2 / (3 cos 2 (#i2 



d = --^-^-( v V ^) (2.72) 

are the heteronuclear and homonuclear effective dipolar couplings respectively. However, 
here the angular brackets in the equations indicate that the dipolar couplings in liquid 
crystals are averaged over both molecular tumbling and any internal bond rotations. 
Thus, the molecular variation of the spin-spin distance, r^, and the angle between the 
internuclear vector and the external field, 9ij, are taken into account. Because of a 
special geometry of the oriented nCB liquid crystals, we will consider two different cases 
where the dipolar constants are related by b\ = 62 = b and b\ = —b 2 = b. 

As in the previous section, for a standard CP experiment, one can neglect the reso- 
nance offsets, and considering that 

|fii,/ + fii,s| > \b k \,\d\, (2.73) 
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the truncated Hamiltonian can be written as 



H s = ^[S X + I X ) + ±A [I* -S x )+}_^ h(S z F k + S y I k 



fc=l,2 



(2.74) 



with 



-h(£li s + £li,i) and A = H (Q 1>s - Oi,/) ■ 



As the Hamiltonian (2.74) has only Zeeman fields along the x direction, we change the 

x and z 



names of the axis as we did in section 
Hamiltonian becomes 



2.1.1 



x 



z, y 



(2.75) 

ge the 
y. Hence, the 



h s = |s (s a + r) + |a (i* -s z *) + Y, b k (s x i x k + s y if) - \d{2i z x r 2 - r x r 2 

fc=l,2 



nil). 

(2.76) 

In eq. (2.76) the non-secular elements of the dipolar interaction with respect to the 



I z + S z j term have been neglected. Similar as in the previous section, this allows us 
to write the matrix representation of the Hamiltonian in a simple block structure using 
the basis {|Mj, Ms)}, with Mi = Mi + M 2 and Ms denoting the spin projections of 
the / and S systems in the direction of their respective RF fields. Now, each block is 
characterized by the total spin projection M = Mj + Ms, i-e. nonzero matrix elements 
exist only between states with the same magnetic quantum numbers M. Thus, the 
heteronuclear dipolar Hamiltonian has non-diagonal terms different from zero generating 
transitions between spin states {|Mj, Ms)} and {|Mj ± 1, Ms =F 1)} • The eigenstates of 
this Hamiltonian can be denoted in the form \M,um), with um = 1, ••,5'M) where gu is 
the degeneracy of M (n± 3 / 2 = 1 and n ±1 / 2 = 1, 2, 3). It is very interesting to note that 
in each space of M — ±1/2 there are only two of the three eigenstates that are involved 
in the dipolar transitions that give rise to the oscillations. This is a consequence of the 
symmetry of the system, i.e. the flip-flop can occur only between the carbon and one 
(the symmetric or the antisymmetric) combination of the proton states depending on 
the relative signs of the heteronuclear couplings (pi = 62 or bi = —62)- The symmetric 
and antisymmetric combination of the proton are 

(2.77) 
(2.78) 

where the vectors are denoted by \M\, M 2 ) . Hence, in the ordered basis 

{\S) ® |+) ; |+, +) <g> I-) ; \A) ® |+)} (2.79) 
with \Mi, M 2 ) (gi \Ms), the M — 1/2 block of the system Hamiltonian is given by 



\s) = i/V2 (|+, -> + |-+» 

L4) = l/V2(|+, ->-|-+» 



n 



S,M= 



\ (S - A) + |d f fa + b 2 ) 
f (bi + b 2 ) (i£ + §A) - \d f(b 2 -h) 
4{h-h) \(Z-A) 



(2.80) 
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and similarly for the M = — 1/2 block. It is evident from the previous equation that un- 
der the conditions b± = b 2 or b\ = —b 2 one of the states, the symmetric or antisymmetric, 
is involved in the dipolar transition. 

The Liouville-von Neumann equation |Abr61t IEBW91] for the density matrix of the 



system is (2.13) 



dt 



P(t) 



h 



(2.81) 



where, similarly as in the 2-spin case (2.18), the initial density operator p(0) considering 
the situation after the it/ 2 pulse in the / system is given 



P(0) 



1 + f3 B htt jP 

Tr{i} 



(2.82) 



The solution of this equation is given by 

p(t) = U(t)p(0)U-\t), 



(2.83) 



where U{t) = exp y—^TCst 

In the simplest case, where the Hartmann-Hahn condition is exactly fulfilled, A 
the exact solution for the evolution of the observed magnetization Mgz (t) is 



M s *(t)=Ti\§ s &(t)\ = M Q f 



[1 — cos(u;ot)] 



where 



with 



u 




K 



if h 
if £>! 



and 



-b 2 
b 2 = 



f 



(2.84) 



(2.85) 



(2.86) 



The natural frequency u^o of the polarization transfer corresponds to the transitions 
between the eigenstates mentioned above. Now, it is clear that the symmetry of the 
system manifests directly in the frequency, where the difference between the two situa- 
tions is represented through the k parameter. So, the relative signs of the heteronuclear 
couplings lead to a characteristic contribution of the homonuclear coupling being three 
times bigger when b\ = b 2 than when b± = —b 2 - The constant 



M = 



(2.87) 



corresponds to the initial magnetization of one / spin. Eq. (2.84) shows that the 
magnetization of S is attenuated by the factor /, and it takes its maximum value Mq 
when d = 0, i.e. when there is no I\ — I 2 interaction. The fact that the homonuclear 
interaction decreases the transferred magnetization was already noticed in ref. [PUL96J. 



We can see that the constant term in eq. (2.84) is proportional to the differences in 



2 Rcmcmber that this initial condition is given in the high temperature approximation. 
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populations between the relevant eigenstates of the system, while the time dependent 
term corresponds to the coherences representing the transitions from {|Mj, Ms)} to 
{\Mj ± 1, Ms =F 1)}- The magnetization in the spins I\ and I2 is given by 

M q (t) = Tr{ifa(t)\ \ f 1 [1 - cos (o; t)] 1 „ 

L \ \ = Mj« (t) = M <l - - f > = M - -M s - (t) 

M q (t)=Tr{iia(t)} J ' U I 2 7 2 / 2 U 

(2.88) 

and the total magnetization by 

M tot (t) = Tr { (> + It + J|) o- (t) } = 2M . (2.89) 

Thus, the total magnetization is given by the initial state and the mean magnetization 
in each site is given by M tot (t) /3 = §M . Because of the symmetry of the system, each 
of the proton transfers forth and back the same polarization to the carbon-13 that is 
half magnitude of the magnetization observed at site S. Figure [273 shows typical curves 



of the Ms- (t) and Mjz (t) magnetization. There, we can see curves for different values 
of the ratio d/b, i.e. different factors /. (the higher the ratio d/b, the lower the value 
of /). The red lines show the difference, with the same values of b and d, between the 
evolution with b\ = b 2 (solid line) and 61 = —62 (dashed line). The mean magnetization 
at each site is show by the blue line. 

2.2.2 A three-spin system coupled to a spin-bath 

In this section, we add to the three-spin system some interaction with other spins using 



an extension of the model proposed by Miiller et al. [MKBE74], see section § 2.1.2 The 
model assumes that the dipolar interactions of the S spin with the I spins are neglected 
except for the coupling to I\ and I 2 . The interaction of these particular spins with the 
bath or the infinite reservoir of I spins is considered in a phenomenological way. All 
kind of spin-lattice relaxations are neglected. The system-environment (SE) interaction 
Hamiltonian can be represented by 



£ SE = Yl [ a W + p ( w + 

[ai° k F° k + |/3 + i+F k ~)] (2.90) 



fc=l,2 



E 

fc=l,2 



with 



^ = u = x,y,z (2.91) 

l>2 

and 

$t = (k ± ^ ) H = Pi (2.92) 



where the subscript I corresponds to the spins within the bath. As in section § 2.1.2 
is an Ising interaction if ft /a = and a XY, isotropic (Heisenberg) or the truncated 
dipolar interaction^] if a/ ft = 0,1,-2 respectively. Following the procedure of section 

3 It corresponds to the fact that we have neglected the non-secular terms with respect to the RF 
field. 
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M/f), dlb= K=3 
M f (i), dlb= K=3 
M s (t),dJb= 1 k=3 
M s {t),dlb= 1 k=1 
M g (t), dlb=\0 Krf 
M f {t),dlb= 1 Krf 




thlb 



Figure 2.3: Temporal evolution of the polarization in a 3-spin system. The black, red and 
blue lines are the 5* (solid) and / (dotted) magnetization with b x = b 2 = b = 3) under 
the ratio d/b = 0, 1 and 10 respectively. The thin red line represents the 5* magnetization 
with d/b = 1 and b\ = —b 2 = b (n = 1) . The blue line is the mean magnetization at 
each site M tot (t) /3 = 2/3M . 



§ 2.1.2, in the semi-classical theory by tracing on the bath variables, F% are treated as 
temporal functions F% (t) representing classical random processes. However, the NMR 
experimental conditions justify a high temperature approximation, and hence the semi- 
classical theory coincides with a quantum treatment |Abr61j . Then, the random SE 
interaction Hamiltonian is written as 

£se (t) = [ aF k (*) K + ¥ { F k (t) It + Fl it) If)] . (2.93) 

fc=l,2 

Here, the interaction of the system with the spins of the bath has been taken into 
account. Any influence of the bath coming from others degrees of freedom (rotations, 
translations, etc.) will manifest through this interaction. These random processes satisfy 



(2.94) 
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where the bar denotes time average, and their correlation functions are 



9 k 



(u,v) 



T) = F£(t)F»* (t + r). 



(2.95) 



The dynamics of the reduced density operator a (t), following the usual treatment to 
second order approximation, is |Abr61j IBlu81| IEBW91] 



dt 



o (*) = ~[Hs, a (t)) - ir {a (t) - a^} . 



(2.96) 



The relaxation superoperator T generated by TCse (t), that accounts for the dissipative 
interactions between the reduced spin system and the bath, drives the density operator 
towards its equilibrium value a^. 

In the following, we assume that the correlation times of the fluctuations are ex- 
tremely short compared with all the relevant transitions rates between eigenstates of 
the Hamiltonian, i.e. frequencies of the order of S/2 and uq. In this extreme narrowing 
regime we obtain 



k u,v 



where 



(u,v) 



dr 
~h 



g { k' V) (i")exp{-i^r} 



is the spectral density and 

€u,v = (a$u,z + (3 {8 u ,x + &u,y)) (oc8 v>z + (3 (8 ViX + 5 v>y )) . 
If we suppose the spatial directions statistically independent, i.e. 

g { r ] (r) = if u ± v, 



(2.97) 



(2.98) 



(2.99) 



the superoperator F can be written as 



where 



fc=l,2 



TX 

1 ki 



T x 2- 



4V 



J k u = ^Jt' u) (o) • 



, (2.100) 



(2.101) 



Here, as in the two-spin case, the axial symmetry of 7is around the z axis leads to the 
impossibility to evaluate separately J£ and J%, so they will appear only as the averaged 
value 

J? = {JZ + JZ) (2-102) 

Taking into account the symmetry of our system b\ = ±b 2 , an extra simplification can 
be done by 



2.2 Three-spin quantum dynamics 



37 



Thus, we obtain 



fc=l,2 



a 



+ P 2 J 



T x 



+ 



jy 



a 



(2.104) 



Although we could absorb the constant a 2 and /3 2 in J7 Z and J xy respectively, we will 



keep it to emphasize the different sources of the anisotropy in eq. (2.104). As we dis- 
cuss in section § 2.1.2 the most usual approximation is to consider J x = J y = J z 
(identical correlations in all the spatial directions) and a = (3 = 1 (isotropic interaction 
Hamiltonian) [MKBE74], however, a better approximation considers a dipolar interac- 



tion Hamiltonian, i.e. a = —2(3 = —2. As in eq. (2.47), we define 

p 2 J xy . 



zz 



a 2 J z 



and 



XY 



(2.105) 



2.2.2.1 Neglecting non-secular terms 



Following the formalism in Abragam and Ernst et al. books |Abr61l IEBW9lj that was 

used in section § 2.1.2.1 we write the superoperator V using the basis of eigenstates of 
the Hamiltonian (2.74). After neglecting the rapidly oscillating non-secular terms with 
respect to the Hamiltonian, i.e., r X Y,r zz <C |6| , a block structure results. The first 
block couples the populations and off-diagonal elements with AM = 0, Zero Quantum 
Transitions (ZQT), of the density matrix. Each of the following blocks couples one order 
AM > 1 of off-diagonal elements of the density matrix among themselves. Because the 



Hamiltonian (2.74) does not have degenerate eigenenergies, all the non-diagonal terms 



coupling the population block with the ZQT block are non-secular and can be neglected. 



As the initial condition (2.82) does not contain coherences with AM > 1, we only need to 



study the evolution of the density operator into a Liouville space restricted to populations 
and ZQT. When there are no degenerate transitions, the secular ZQT block is diagonal. 
However, in our case there are degenerate transitions between eigenstates within the sets 
with M = ±1/2. Thus, some non-diagonal terms in the ZQT block cannot be neglected. 

In the final condition, the S7 2 system reaches the temperature of the I spins reservoir 
as was described in the section 



2.1.2 



(7n 



i + (3 B hn QJ (s z + r 
Mi} 



(2.106) 



It is easily seen that commutes with 7is, not containing coherences with AM > 1. 

By using the present formalism under the considered approximations, we will solve 
eq. (2.39) for the cases relevant to our liquid crystal study. 



Isotropic system-environment interaction rate. Considering 

r Z z = r XY = r (j z = j xy and a = p = i), 

the time evolution of the S magnetization results 

M s * [t) = M [1 - A + e~ R+t - A^e- R - 1 - A c cos (oot) e~ Rct ] , 



(2.107) 
(2.108) 



38 Chapter 2. Many-spin quantum dynamics within the density matrix formalism 



where 



u = UJ , 

R± = x±^/K 



with 




/-I +/ 



(2.109) 
(2.110) 
(2.111) 



(2.112) 
(2.113) 



and 



A, 




± 



[f/(l 



I) 



1] 



(!/-!)+/ 



(2.114) 
(2.115) 



The figure 2.4 shows typical curves of eq. (2.108 ) for different values of the ratio Y / (uqK) . 



As we observed in the 2-spin case, we see that the oscillations are attenuated when the 
ratio is bigger and the net transfer of polarization is faster. The black lines compare 
two different curves with / = 1 (thick line) and / = 0.7 (thin line) for a fixed value 
of Y / (uoh) . We observe that the maximum of the oscillation is smaller as / decreases 
but the final magnetization is the same for both curves. The figure 2J3 shows the 
dependence of the coefficients Ai and the relaxation rates Ri as a function of /. Notice 
that Xi'Xc — 0- Using the initial condition Ms* (0) = 0, it is easy to see that the 
positive constants A + , A_, A c satisfy 1 — A + — A_ — A c = 0. In general A + <C A_ and 
R + > R-, so the first exponential term can be neglected as can be observed in fig. 2.6 
This approximate solution is excellent for / <C 1, but even in the worst case (/ ~ 1) 
differs about 7% from the exact solution (see fig. 2.6[ ). 

The first maximum in the magnetization Ms* (i) is approximately / Mq and the os- 
requency ujq as can be observed comparing with the isolated evolution in 



it 



cillat 


ion 


has 


figs. 


2.4 


and 



2.6 



The oscillations have an amplitude ^e~ Rct that represents the atten- 
uation of the coherences of the SI2 system mounted over non-oscillatory terms. These 
terms take into account the effect of the bath, not only by transferring magnetization 
but also breaking coherences and leading to a quasi-equilibrium. This quasi-equilibrium 
state is given by 



a qe = exp{-H$/(k B T qe )}/Tr{exp{-H s /(k B T qc )}} 



with 



T qe = 3/2(O li7 /^ 0l /)T 
the temperature of the three-spin system. 



(2.116) 
(2.117) 
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Figure 2.4: Typical curves of the magnetization of the S spin in a three-spin system 
coupled to a spin-bath under an isotropic SE interaction. The thick red, black and 
green lines are the polarization evolution with / = 1 and the ratios Y / {(jJqK) = 0.01, 
0.06 and 0.2 respectively. The black thin line represents the magnetization for / = 0.7 
and Y / {uiqK) = 0.06 and the light and dark gray dotted lines are the isolated dynamics 
for / = 1 and 0.7 respectively. 



In the particular case when / = 1, i.e. no I1-/2 interaction, A± = \ T J> 

A c = f/2, R± = \ (9 ± y/Vf) Y/h , R c = 7/4Y/H and u = V2b/h, showing that only 
under this condition the frequency given in ref. [PR 96J is valid. But even under this 
condition, our results show that the equation obtained by Miiller et al. for the SI case 
cannot be directly applied to the SI 2 system. In this last case the attenuation of the 
oscillations and the transfer of polarization to the system is slightly faster than in the 
SI case. 

Anisotropic system-environment interaction rate. Considering 

Tzz + r XY , (2.118) 

we obtain 

M S z (t) = M [1 - A ie - Rlt - A 2 e~ R2t - A 3 e~ R3t - A c cos (ut) e~ Rct ] , (2.119) 



40 Chapter 2. Many-spin quantum dynamics within the density matrix formalism 




Figure 2.5: Coefficients Ai [panel a)] and relaxation rates Ri [panel b)] of the S polar- 



ization expression (2.108) for an isotropic system-environment interaction as a function 



of the relation between homo and heteronuclear spin couplings /. 



where the Ai, 2 = 1,2,3 are functions of / and Tzz/^xy and 



uj = uj , 

R c = (2-f)T XY /h 
1 



Ac 



J- 



(2.120) 
(2.121) 

(2.122) 



The expressions for Ai = Ai(f, Tzz/^xy) and Ri = Ri(f, Tzz/h, T-xy/K) are too long to 
be included here but they are available as supplementary material. 

The transfer of polarization from the bath to the system depends on the non- 



oscillatory terms of eq. (2.119). In the r zz /r X Y > 1 case, at long times (R c t » 1), 
only one of the three exponential terms contributes. In this regime, the transfer is es- 
sentially given by Txy, although there is a slight dependence on Tzz- This differs from 
the SI behavior where the polarization transfer from the bath depends exclusively on 
Txy (see section § 2.1.2). This is a consequence of the fact that in the SI system the 



quasi-equilibrium S z polarization, (1/2)M (the mean magnetization at each site), co- 
incides with the time averaged value of the isolated system. As Txy is associated to the 



flip-flop term in the SE interaction Hamiltonian (2.90), its role transferring polarization 
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Figure 2.6: Polarization evolution of the S magnetization (black lines) of a 3-spin system 
coupled to a spin-bath. The figure shows comparisons between Ms* (t) and the aproxi- 
mate solution (red line) where was neglected one of the pure exponetial term. The pure 
exponetial terms that are contained in the M^z (t) expression are showed by the blue 
and green line. It is evident why the term described by the green line could be neglected. 
The light gray line shows the isolated evolution of the S magnetization to observe that 
the first maximun is essentially /Mo. Panel a) shows the evolution for f — 1 and panel 
b) for / = 0.2 where in the last one the aproximate solution is more better. 



can be easily interpreted. The effect of Tzz is more subtle, it can be associated to a pro- 
cess where the environment observes the system breaking its coherences. This process 
that involves the operator F Z I Z in 7Yse, which is a like number operator, is discussed in 
chapter [3] 



Figure 2.7 shows typical curves of the S polarization (2.119) for different SE inter- 



actions and values of / (the relation between b and d). In the limit T-iz 3> Txy [highly 
anisotropic case, panel d)] it is possible to distinguish two time regimes: one in which 
the system decoheres and reaches a quasi-equilibrium, [7ig,cr qe ] = [SHE98J, character- 
ized by Tzz and other in which polarization transfer from the bath is completed with 
a rate r X y/^- In this situation, it is possible to see that R c is proportional to r Z z, as 
it occurs with the Rt corresponding to the fastest exponential terms. The dependence 
of the non-oscillatory terms on T Z z observed in the SI 2 system can be assigned to the 
fact that the quasi-equilibrium carbon polarization (2/3) Mo (the mean magnetization 
at each site) does not coincide with its time averaged value (//2)M in the isolated 



three-spin system [eq. (2.84) 



Equations (2.108) and (2.119) will be used to fit the experimental data in order to 
extract the relevant parameters of our system. 
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Figure 2.7: Typical curves of the S polarization (2.119) for different SE interactions 



and values of / (the relation between b and d). The dipolar, XY and isotropic SE are 
showed in panel a), b) and c) respectively. In the limit Tzz ^> I\y [highly anisotropic 
case, panel d)] it is possible to distinguish two time regimes: one in which the system 
decoheres and reaches a quasi-equilibrium, [H s , a qe ] = 0, and other in which polarization 
transfer from the bath is completed with a rate I\y / h. 



2.2.3 Many-spin quantum dynamics during Cross-Polarization 
in 8CB 

Nuclear Magnetic Resonance experiments were carried out by people of our group in 4- 



n-octyl-4'-cyanobiphenyl, also called 8CB (see fig. 2.8) obtained from Sigma (Chemical, 
Co) and used without further purification. This system presents the mesophases smectic 
A (SA) and nematic (N), with transition temperatures at 294.5 K (K-SA), 306.5 K (SA- 
N) and 313.5 K (N-I). 

1 H- 13 C cross-polarization measurements as a function of contact time t c were per- 
formed in the smectic and nematic mesophases. In all the cases, the CP sequence was 
performed in static conditions and combined with the sequence SPINAL-64 to perform 
an efficient proton decoupling during acquisition without appreciable heating of the 
sample jFKEOOj . 

In the smectic phase, standard CP experiments were performed at 300 K in a Bruker 
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Figure 2.8: Chemical structure of 4-n-octyl-4'-cyanobiphenyl (8CB) showing the labels 
used in the 13 C-NMR spectra. 



MSL-300. The acquisition time was 92 ms, with 60 ms for decoupling and a recycling 
time of 15 s. The Hartmann-Hahn (H-H) condition [HH62, Sli92] was set with an RF 
field amplitude for carbons corresponding to Qi t s/2ir = 67.7 kHz. During the experiment 
the contact time was varied in the range 2 fis < t c < 5 ms. 

In the nematic phase at 311.5 K two types of CP experiments were performed in a 
Bruker AVANCE DSX-500. The first was a standard CP with protons on-resonance. 
The second was a CP experiment with irradiation for protons in the Lee-Goldburg (LG) 
condition, i.e. the off-resonance for protons was set to have an effective field at the 
magic-angle with the static field H . The acquisition and decoupling times were 74 ms. 
In the standard CP, the H-H condition was set with Qi^/^tt = 60.3 kHz while in the 
Lee-Goldburg Qi f s/2n = 74 kHz. In both sets of experiments the H-H condition was 



optimized for C("y) (see figs. 2.8 and 2.9), and the contact time t c varied up to 2 ms. 

In all these experiments the temperatures were calibrated using the N-I temperature 
transition. Also a 13 C-NMR spectrum of 8CB in the isotropic phase (at 320 K) was taken 
as reference using a single pulse sequence with X H decoupling (direct 13 C polarization). 

2.2.3.1 Comparison between experiments and theoretical results 



The 13 C-NMR spectrum of 8CB in smectic phase can be seen in fig. 2J3 In the inset 
of this figure, the alkyl part of the spectrum at a temperature corresponding to the 
nematic mesophase is displayed. The aromatic part of the spectrum keeps the same fea- 
tures up to the nematic-isotropic transition temperature. The position of C(u), methyl 
group, which in the isotropic phase is at 14.1 ppm has been taken as reference because, 
as experimentally observed, it does not vary at the working temperatures due to its 
high mobility. The aromatic part of the carbon spectrum for 8CB is similar to those of 
other members of the nCB series, so we consider the up-dated assignments reported for 
5CB in previous works [FKEOOj. A detailed temperature dependence of the 13 C chem- 
ical shifts in 8CB as well as their complete assignments have been reported previously 
|CCH+02j . The assignments for the alkyl carbons are supported by the segmental or- 
der parameters obtained from experimental CP frequencies (see below), and deuterium 
NMR experiments |Don971 [CEHL85] . 

In order to analyze the experimental data, it is necessary to correlate the geome- 
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Figure 2.9: 13 C-NMR spectra of 8CB in the smectic mesophase at 300 K. In the inset we 
see the aliphatic part of the spectrum corresponding to 311.5 K (nematics), which is the 
only part that changes at the working temperatures. The assignments for the aliphatic 
part of the spectrum have been done in previous works. Labels refer to fig. 2.8 



try/symmetry of the molecular interactions with the cases presented in the theoretical 
section associated to different values of k, eq. (2.86). For each methylene carbon in 
the aliphatic chain, the geometry of the molecule and the rapid rotations around the 
C-C bonds, which lead to the trans-gauche isomerizations, allows us to take a single 
averaged value for both heteronuclear couplings i.e. b\ ~ b 2 . Then, for carbons C(a) 
to C(r)) in the aliphatic chain the relation of signs of the heteronuclear couplings corre- 
sponds to the case where k = 3. It is also possible to see from simulations done with a 
similar molecule (4-n-pentyl-4'-cyanobiphenyl or 5CB) and from geometrical considera- 
tions that the homonuclear dipolar interaction d between protons belonging to the same 
methylene in general, cannot be neglected [SKSM01J. For each non-quaternary carbon 
in the phenyls rings C(3'), C(2'), C(2) and C(3), we have one directly bonded proton 
interacting with a dipolar coupling b\\ however, a careful analysis indicates that neither 
the 1 H- 1 H interaction d nor the coupling between the 13 C and the nearest non-bonded 
1 H, 6 2 , can be neglected. Considering the rigidity of the phenyl ring and the orientation 
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of each internuclear vector with respect to the external field, we see that b\ ~ —62- Then, 
each non-quaternary aromatic carbon is part of a three-spin system, where both het- 



eronuclear couplings can be considered having an averaged mag nitude b = y/(bl + bj) /2 
and different relative signs, leading to k = 1. 

Firstly, let's consider the experimental results corresponding to the smectic phase. 
In this case, typical oscillations and relaxation of the 13 C polarization vs. contact time 
t c are shown in fig. 2.10 We remark that neither the quaternary carbons nor C(uj) show 
oscillations in the CP experiment. In the last case, the reason is the high mobility of the 
methyl group that averages to zero the effective carbon-proton interaction. In all these 
experiments, the 13 C polarization has essentially reached its asymptotic value at 5 ms. 
The absence of a decay in these curves allows us to neglect spin-lattice relaxation in the 
rotating frame (T" lp ) in the time regime analyzed. We also note that the CP frequencies 
corresponding to methylene groups are higher than those of the aromatic cores. This is 
due to a particularly unfavorable angle (~ 60°) between the internuclear carbon-proton 
vector and the external magnetic field in the case of the phenyl rings. Besides, the 
contribution of the homo nuclear coupling to the frequency is much smaller when k — 1. 



In each dynamical curve shown in fig. 2.10 the C polarization at the first maximum 



is lower than its asymptotic value. As we have seen in the theoretical section, this fact 
indicates that the homonuclear coupling is not negligible and it will allow us to evaluate 
it. 

We have fitted the experimental CP data to the equations derived for the isotropic 



(eq. (2.108)) and anisotropic models (eq. (2.119)) presented in section § 2.2.2.1 In 
the last case, we have distinguished a purely dipolar anisotropy (a = —2/3 = —2 and 
jxy _ jz ^ j e r zz /r XY = 4) from the most general case. For the smectic phase, 



it is seen in fig. 2.10 that the isotropic model separates from the experimental points 



after approximately the time of the second maximum, fitting very poorly the asymptotic 
behavior. The dipolar model constitutes an improvement over the isotropic one, without 
adding extra free parameters. However, a much better fitting is obtained using the 
anisotropic model which follows very closely the first oscillations of the magnetization 
and it is the best suited in the asymptotic regime. 

In the case of the standard CP experiments performed in the nematic phase up 
to 2ms (not shown in the figure), the behavior is already well fitted with the purely 
dipolar model, not showing any appreciable improvement by the use of the completely 
anisotropic model. 

Taking into account the fitting parameters corresponding to each experiment, we 
can obtain the effective dipolar couplings and the relaxation constants for each carbon 
showing oscillations during CP. In tables 2A_ and 2J2 we show the heteronuclear and 
homonuclear dipolar couplings at different temperatures. The values for the heteronu- 
clear couplings are comparables with those obtained for 8CB in ref. |FPG + 86] , where a 
different experimental technique was applied. 

It is important to emphasize here, that the frequency of the oscillations uj is a quite 
independent and robust parameter, leading to values which vary less than 4% using 
different models. This allows one to determine the heteronuclear coupling b with very 
small error. The homonuclear coupling, however, is much more sensitive to the relation 
between the oscillatory and the asymptotic regimes and so, more dependent on the 
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oarbons 


kHz 


kHz 


C(3') 


1.30 ±0.03 


4.9 ±0.3 


C(2') 


1.49 ±0.04 


4.8 ±0.3 


C(2) 


1.58 ±0.06 


4.9 ±0.3 


C(3) 


1.45 ±0.04 


5.3 ±0.4 


c(P) 


4.6 ±0.1 


4.5 ±0.6 


c(t) 


4.8 ±0.1 


5.9 ±0.7 


C(5) 


4.7 ±0.3 


2.9 ±0.3 


C(e) 


4.3 ±0.2 


4.0 ±0.5 


G(0 


3.1 ±0.1 


3.7 ±0.4 


Cfo) 


3.2 ±0.1 


2.2 ±0.3 



Table 2.1: Effective heteronuclear and homonuclear dipolar couplings b and d obtained 
by fitting the data of the standard CP experiment performed at 300K (smectic phase) to 
the anisotropic model. The signal corresponding to C(a) does not appear in the smectic 



spectra. The errors have been assigned taking into account eqs. (2.123) for b, d in terms 
of the fitted parameters. 



Carbons 


Standard CP 
kHz 


LG-CP 
kHz 


Standard CP 
kHz 


LG-CP 
kHz 


C(a) 


4.0 ±0.2 


4.18 ±0.04 


2.7 ±0.5 


2.9 ±0.6 




2.7 ±0.1 


2.95 ±0.03 


3.5 ±0.6 


2.8 ±0.6 


c(t) 


3.0 ±0.2 


3.34 ±0.03 


3.5 ±0.7 


2.5 ±0.5 


C(e) 


2.8 ±0.1 


2.95 ±0.03 


3.1 ±0.5 


2.8 ±0.6 


c(?) 


2.5 ±0.1 


2.29 ±0.03 


2.4 ±0.5 


2.8 ±0.5 


c(v) 


2.4 ±0.1 


2.17 ±0.03 


1.5 ±0.3 


2.3 ±0.4 



Table 2.2: Effective heteronuclear and homonuclear dipolar couplings b and d for the 
aliphatic carbons obtained from the experiments performed at 311. 5K (nematic phase). 
In the standard CP experiment the values are obtained by fitting the data to the purely 
dipolar model (see text). In LG-CP the values of b are obtained directly from eq. 



(2.124) while the values of d requires a combination of the standard CP with the LG-CP 



frequencies (eq. (2.127)). In this mesophase the signal of C(/3) merges to the signal of 



C(5), so the same value of coupling has been assigned for both carbons. 
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Figure 2.10: 13 C polarization as a function of contact time t c for aromatic and aliphatic 
carbons in a standard CP experiment at 300 K (smectic phase). Fittings of the exper- 
imental data to the anisotropic, isotropic and purely dipolar models described in the 
text, are displayed. 



model and the extent of the experimental data. To see clearly this fact, we notice that 
the frequency and the attenuation factor / given in eq. (2.85) are the fitting parameters 



that yield b and d regardless of the model. The dependence between these parameters 
and the couplings are 



1 b 
2nh 



1 u 



./ 



and 



1 d 
2irh 



4 u 

K 2lT 



(2.123) 



As in our cases, / is always bigger than 1/2, the error in / (approx. 2% for the experiment 
in smectic and 5% for the experiment in nematics) affects less to b than d. Now, taking 
into account that the CP frequency to is the best parameter, fitted with an error lower 
than 2%, we can assign the error of the other parameters of interest. Then, the coupling 
b is obtained with a relative error between 3% and 5% while the relative error of d 
reaches in some cases a value of about 20%. 



Fig. 2.11 displays the tendency of the CP frequency obtained for the different exper- 
iments. In particular, for the aliphatic carbons, the expected zig-zag pattern is observed 
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|Don97l ICEHL85] . For the aromatic part of the molecule we have u}(3') J$ uj(2') ps 
to (2) ps a? (3), this is expected because C(3/), being the closest to the cyano group forms 
a bigger C-C-H angle, giving rise to a smaller dipolar interaction [CHEP96]. On the 
other hand, when increasing temperature, i.e. going from smectic to nematic, a fur- 
ther averaging of the dipolar interactions occurs, manifested in the decrease of the CP 



frequencies. For comparison, fig. |2.11| shows the results of the CP experiment at the 
Lee-Goldburg condition, performed in the nematic phase at the same temperature. In 
the LG experiment the CP frequency (ujlg) is only related to the heteronuclear coupling 
|Sli92] because the homonuclear contribution has been quenched, then 



^lg = V 2 (W 2 (sen^ m ) 2 , (2.124) 

where 8 m = 54.7° is the magic angle [Sli92j. This angular factor comes from the pro- 
jection of the RF field into the direction of the effective field for protons irradiated 



off-resonance at the LG condition |Sli92] . Then, considering eq. (2.85), we expect the 
relation 

lj > cjLG/sen^ m . (2.125) 



As shown in fig. 2.11 this relation is not satisfied for aromatic carbons. Consid- 
ering that C(7) is exactly on- resonance and that the H-H condition was optimized for 
that carbon, the disagreement can be attributed to a non negligible mismatch of the 
Hartmann-Hahn condition for carbons in the aromatic part of the spectrum. Taking 
into account a mismatch A for the aromatic carbons, the CP frequency for the LG 
experiment becomes, 



u LG = ^/2(&/ft) z (sen# m ) 2 + A 2 . 
In the standard CP experiment, two frequencies appear, where the observable one is 



u,= U (2.126) 



with 



w± ^ /M ^ + 2(Vfi)2 



The other modulating frequency is too low to be observed in the presence of relaxation. 
Although it is difficult to quantify exactly the magnitude of A, we can see that the 
effect of the mismatch is greater for uj^q than for u where there is a partial compensation 
between the two contributing frequencies uj . This explains why o>lg ~ w for the carbons 
in the aromatic part, in contrast with eq. ( |2.125 ). 



For carbons irradiated on-resonance (aliphatic part), the values of the fitted pa- 
rameters h and d obtained from standard CP experiments can be compared with the 
parameters obtained from the LG-CP. In the later case, the parameter of interest is 



ujlg- Then, b is obtained in a direct way from expression (2.124) and d can be obtained 
comparing o;lg with the corresponding value of u at the same temperature. Using both 
experiments, d is calculated from 
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Figure 2.11: Cross-Polarization frequencies obtained from the LG and standard CP 
experiments for all the non-quaternary carbons in the molecule. Note that the signal 
corresponding to C(cu) does not appear at 300 K; also note that the signal of C(/3) merges 
with that of C(S) at 311.5 K so the same value of frequency has been assigned to both 
carbons. 



The values of d calculated in this way with k = 1 (aliphatic carbons) can be compared 
with those obtained by fitting eq. (2.119) to the standard CP data, i.e. coming from a 
single experiment. 

Table 2J2 displays the values of the homonuclear and heteronuclear couplings for the 
aliphatic carbons obtained from the standard CP experiment, and combining this with 
the LG-CP performed at the same temperature. Fig 2.12 allows for the comparison of the 
13 C- 1 H and 1 H- 1 H couplings obtained from the two methods. As expected, an excellent 
agreement can be observed for the values of b. The novel methodology to estimate 
d values seems to yield good results within an error of around 20%, which could be 
easily improved by taking longer time data. However, there are no many experimental 
estimations of these homonuclear couplings and those obtained directly from the 1 H 
spectra are in good agreement with the values obtained here. 

With regard to the relaxation of the system in the CP experiments, we can observe 
an anisotropic behavior for both phases, quantified by the ratios Tzz/Txy > 1 (see fig. 



2.13). Besides, we note that there is an important difference between the behaviors of 
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Figure 2.12: Effective homonuclear (l/2ir) (d/h) and heteronuclear (1/2tt) (b/K) dipo- 
lar couplings for carbons in the alkyl chain, obtained from the LG and standard CP 
experiments performed at 311.5 K (nematic mesophase). 



the nematic and smectic phases. In the nematic phase, the anisotropy can be explained 
with a purely dipolar model. The average anisotropy factor for all the carbons in the 
molecule is (4 ± 1). In contrast, in the smectic phase the factor Tzz/Txy > 4 reveals a 



highly anisotropic behavior for most of the carbons. This can be appreciated in fig. 2.10 
where the amplitude of the first maximum is higher than the following ones, specially 
for the aliphatic carbons. This fact is not observed in nematics, giving support to 
the purely dipolar anisotropy. Different factors can produce this high anisotropy, one 
could be Hi inhomogeneity, another could be that our main system is actually larger 
than S7 2 . Although both factors would effectively increase the anisotropy, the effect 
should be comparable in both phases. Moreover, under the same Hi inhomogeneity, we 
have not observed such anisotropy ratios in molecular crystals. Alternatively, the large 
anisotropy observed in the smectic phase can be originated on the lack of fast fluctuations 
in this more rigid phase, which would prevent the application of the extreme narrowing 
approximation. If we release this approximation and assume that the spectral densities 
J7"(0) and j7(S/2), although different between them, are approximately constant in an 
interval of width 2c<Jo, our calculations indicate that Tzz is proportional to J7"(0), while 
Txy J7"(S/2). As usually J7"(0) > j7"(£/2), this could explain the higher anisotropy in 
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the smectic phase, where motion is more hindered than in the nematic phase. 
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Figure 2.13: Relaxation factors T zz and T X y obtained by fitting the standard CP data 



to the anisotropic model [eq. (2.119) 



2.3 Summary 

We solved the spin dynamics of many-spin systems. First, we started solving the dy- 
namics of a two-spin system interacting with a fast fluctuating spin-bath. We extended 
ALP07 the solution of MKBE [MKBE74J starting from the system-environment Hamil- 
tonian to obtain the relax ation su peroperator. We solved the spin dynamics for an 



anisotropic SE interaction [ALP07] . where the anisotropy is given by the ratio between 
the Ising and XY contributions. Miiller, et al. used a phenomenological isotropic relax- 
ation for the SE interaction. Here, we emphasize the different sources of the anisotropy 
in the decoherence and relaxation processes. The main difference is that while the XY 
interaction takes the systems to the tot al system equilibri um, the Ising SE interaction 



takes it to an internal quasi-equilibrium CAL + 03[ IALP07 . In the particular case of an 
isotropic SE interaction, we reproduce the MKBE solution. These solutions were ob- 
tained neglecting the non-secular terms of the SE interaction with respect to the system 
Hamiltonian, i.e. under th e condition b 3> Tse- We extended the solution including the 
non-secular terms ALP07 as will be fully discussed in chapter 3 
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In section § 2.2 we further extended the MKBE model to solve th e spin dy namics 



of a three-spin system interacting with a fast fluctuating spin-bath |CAL + 03] . Our 
calculations were performed for three-spin configurations directly applicable to the 8CB 
molecule. This molecule has two possible configurations for the heteronuclear couplings 
that led us to realize that in each space of M = ±1/2 there are only two of the three 
eigenstates involved in the dipolar transitions that give rise to the oscillations. This is a 
consequence of the symmetry of the system, i.e. the flip-flop can only occur between the 
carbon and one combination of the proton states (the symmetric or the antisymmetric) 
depending on the relative signs of the heteronuclear couplings (pi = b 2 or b\ = —b 2 ) 



CAL + 03j . We showed that this affects the oscillation frequency that one measures 



experimentally. 



In section § 2.2.3 we introduced 13 C- 1 H cross-polarization experiments which com- 



plemented with the detailed spin dynamics calculations allowed us to o btain sepa rately 



the homonuclear and heteronuclear dipolar couplings in CH 2 systems CAL + 03 . The 



reliability of the obtained values was tested with a direct determination of the heteronu- 
clear couplings using CP under Lee-Goldburg conditions. Comparing the standard CP 
experiments with the LG-CP one, we concluded that the last one is better to deter- 
mine directly the 13 C- 1 H couplings. However, the standard CP allows one to obtain 
homonuclear couplings in addition to the heteronuclear couplings. This experiment is 
also important to get further information of the system as for example the relaxation 
phenomena. 

From the theoretical analysis of the CH2 dynamics, we recognized two different time 
scales for the decoherence behavior given by Tzz and Txy- Besides, the CP data showed 
that the rate of attenuation of the oscillations is much faster than the rate of polariza- 
tion transfer from the bath. This anisotropy could be explained in the nematic phase 
by assuming a dipolar interaction Hamiltonian between the three-spin system and the 
bath within the extreme narrowing approximation. In the smectic phase, however, the 
anisotropy is much more pronounced and it seems that the short time fluctuations of 
the bath approximation is not appropriate. Consideration of a slow motion regime leads 
to a better agreement with the experime ntal obser vations without resorting to other 



mechanisms which operate in both phases [CAL + 03 



Chapter 3 



Spin dynamics within another 
perspective: The Keldysh formalism 



As we discuss in previous chapters, the characterization and control of spin dynamics 
in open and close spin systems of intermediate size remains a problem of great inter- 
est [MLL03J. However, the quantum interferences of these systems become damped 
by the lack of isolation from the environment and one visualizes this phenomenon as 
decoherence. Indeed, the inclusion of the degrees of freedom of the environment may 
easily become an unsolvable problem and requires approximations not fully quantified. 
This motivates a revival of interest on previous studies in various fields such as Nuclear 
Magnetic Resonance [MKBE74J, quantum transport [DP90J and the quantum-classical 
correspondence problem |PZ99^ IZur03] with a view on their application to emergent 
fields like the quantum computation |WL02t IDRKH03~1 |TSS + Q0j and molecular electron- 
ics [MKRMlESlN03lEnil[Zim03]. 

The most standard framework adopted to describe the system-environment inter- 
action is the one that we used in the previous chapter, the generalized Liouville-von 
Neumann equation |Abr61t IEBW9lj in the fast fluctuation approximation. There, the 
degrees of freedom of the environment are traced out to yield a quantum master equation 



(QME), eq. (2.39). Interactions with the environment given by the Fermi 

Golden Rule (FGR) providing a dissipative mechanism that could induce a non unitary 
dynamics into the system. An overall (conservation) balance condition is obtained by 



imposing a convergence into the thermal equilibrium state (2.49). While sufficient for 
the most traditional applications, this approximation leaves aside important memory ef- 
fects and interferences in the time domain produced by the coherent interaction between 
the system and the bath |TIL03j . 

A less known alternative is provided by the Keldysh formalism |Kel64j in the integral 
representation proposed by Danielewicz |Dan84] . On one side, it uses the well known per- 
turbation to infinite order in selected terms provided by the Feynman diagrams where, 
under certain conditions that go beyond the Fermi Golden Rule, the dynamical feedback 
effects become relevant. On the other, this integral representation has the advantage 
of being able to profit from a Wigner representation for the time-energy domain. This 
last representation is particularly meaningful in the fermionic case since it allows one 
to define energy states and their occupations simultaneously with the physical time 
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|Pas92] . In that case, one can transform the Danielewicz equation into the generalized 
Landauer-Biittiker equation (GLBE) [Pa s9Tl IPas92] to solve the quantum dynamics of 
the system. The Keldysh formalism already inspired original experimental and theo- 
retical developments in coherent spin dynamics involving quantum interferences in the 
time domain. In particular, it was used to develop the notion of polarization waves 
leading to mesoscopic echoes [PLU951 IPUL96] . to establish the influence of chaos on 



time reversal (Loschmidt echoes) [LUP98] |PLU + 00| IJPOlj and to propose a spin pro- 



jection chromatography |DPL04] . A rough account of many-body decoherence enabled 
the interpretation of anomalies in "spin diffusion" experiments as a manifestation of 
the quantum Zeno effect |PU98j . This technique was applied to a case with an exact 
analytical solution [DPL02J and where more standard approximations can be obtained 
to show the potential of our proposal. Our model represents a single fermion that can 
jump between two states while an external fermionic reservoir is coupled to one of them. 
This environment provides decoherence due to a through space Coulomb interaction 
and can feed the system with an extra particle through tunneling processes. While 
the parameters and approximations involved in this model are especially designed to be 
mapped to a problem of spin dynamics (the two-spin system of chapter [2]), it could also 
be adapted to represent a double quantum dot charge qubit [VMB05] . In that case, a 
double dot is operated in the gate voltage regime where there is a single electron which 
can jump between the two coupled states. Only one of these states is coupled to an 
electron reservoir. We introduce fictitious interactions to obtain a common interaction 
rate which leads to a homogeneous non-hermitian effective Hamiltonian. In the specific 
model considered, we analyze how different SE interactions, e.g. tunneling to the leads 
and through space Coulomb interaction, modify the quantum evolution. A particular 
advantage of the fictitious symmetrization is that it leads naturally to a stroboscopic 
representation of the SE processes. This leads to a very efficient numerical algorithm 
where the quantum dynamics is obtained in a sequence of time steps. Finally, we resort 
to the Jordan- Wigner mapping between fermions and spins to apply the procedure to 
a spin system. This allows us to give a first-principle derivation of the self-energies to 
explain, in chapter |4| the puzzling experimental dynamics observed in a spin swapping 



gate |LUP98j (see fig. 1.6). 



3.1 Two level system dynamics 
3.1.1 The system 

Consider a two fermion state system interacting with the following Hamiltonian Tig, 

H s = E c^c Q + E x c\c x - V 01 (cfcj + c^c ) , (3.1) 

with ^(cj the standard fermionic creation (destruction) operators. The Ei are the ener- 
gies of the z-th local states whose spin index is omitted and Vyi the hopping interaction. 
In matrix representation, if we have only one particle in the system, we have a 2 x 2 
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matrix for the Hamiltonian 

|1,0> |0,1) 

f E -Vox \ |1,0) (3.2) 

Hs ~ {-V 01 E x ) |0,1> - 

Here |1, 0) and |0, 1) denote the states with the particle in the level and 1 respectively. 
3.1.2 System evolution 

We are interested in the study of the evolution of an initial local excitation in the system. 
For definiteness, we consider the initial excitation on site 1 which is described by the 
non-equilibrium state |^(0)) = cf |0, 0) = |0, 1) where |0, 0) is the vacuum state. The 
probability to find the particle in the state and 1 is 

P i(*) = |<l,0|*(t)| 2 , (3.3) 

Pn(t) = |(0,l|*(t)| 2 J (3.4) 

where 

|*(*)) = exp{-i# s *}|*(0)). (3.5) 

Thus, we obtain if E = E\ 

P i (t) = o _ \ cos ( w o<) , (3-6) 

Pn(t) = 2 + 2 COs(a;o ^' (3 ' 7) 
where io§ = 2Voi/fr gives the natural oscillation frequency of the transition between the 



states and 1. The fig. 3.1 show the dynamics of these probabilities. 



3.2 A two level system interacting with a particle 
reservoir 

3.2.1 The system 

Let's consider the electron two-state system of the previous section asymmetrically cou- 



pled to an electron reservoir, as shown in fig. 3.2 a), with the total Hamiltonian 



TC = Hs + He + T~Cse- (3.£ 



The system Hamiltonian is the same described in the previous section (3.1) 



Hs = E c^c + E 1 c J [c 1 - Vqi (c Q ~c 1 + cfc ) , (3.9) 
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Figure 3.1: Evolution of the probabilities to find a particle in site (black line) and site 
1 (red line) with the initial conditions at site 1. 



with c^(Cj) the standard fermionic creation (destruction) operators. The hopping inter- 
action Voi gives the natural frequency, uo = 2Voi/ft, of the transition between the states 

and 1. The environment has a similar Hamiltonian, 

oo oo 

H E = E&tt - V ij ffij + d l d d > ( 3 - 10 ) 

i=2 i=2 

where the determines the topology of the interaction network in the environment 
states. The system-environment (SE) interaction is described by 

^SE = [^12 ^ ^l,/3^1,/3^2,a^2,a + ^12 ^ ^l,a^2,a^2,a^l,a ~ ^12 {pl,ap2,a + &2,apl,ct) ) 

(3.11) 

The first two terms on the rhs represent the Coulomb interaction of an electron in site 

1 with an electron in site 2, the first site of the reservoir. U^ 1 '^ is the standard direct 
integral and U^ cla ^ is the small exchange integral which we include for completeness 
The third term is the hopping interaction between sites 1 and 2. 
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Figure 3.2: a) System-environment (SE) representation. Dashed circles and solid circles 
represent the system and the environment states respectively. Dashed lines are hopping 
interactions while wiggels lines are through-space Coulomb interactions, b) Self-enegy 
diagram summing up the different interactions with the environment in a local basis. 
The lines with arrows are exact Green's functions in absence of SE interactions. The 
double dashed line represents the effective SE interaction, c) Retarded Green's function 
at site 0. The interaction with the environment is to infinite order in the self-energy 
given in b) . d) Particle density function at site 0. The double dashed lines represent 
the effective interactions local in time and space summed up to infinity order. 
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3.2.2 System evolution 

3.2.2.1 Quantum dynamics within the Keldysh formalism 

Now, we have a many-body system and we consider it in the thermodynamical many- 
body equilibrium state |^ e q.)- hi the infinity temperature limit, it has all the states 
equally occupied. We are interested in the study of the evolution of an initial local 
excitation in the system. Thus, as we consider previously for definiteness, let's consider 
the initial excitation with a particle on site 1 and a hole in site which is described by 
the non-equilibrium state 

|* n .e.) =C+C |^ cq .). (3.12) 

The evolution in a complete norm preserving solution is described by the particle and 
hole density functions, 

Gfj (t 2 ,ti) = l ~ (* n .e.| C+ (t x ) C % (t 2 ) |^ n . e .) (3.13) 

and 

G% (t 2 M) = ~ (*n.e.| q (t 2 ) C+ (h) |* n . e .) , (3.14) 

that describe spatial and temporal correlations. In these expressions, the creation and 
destruction operators are in the Heisenberg representation. Notice that, in contrast with 
the equilibrium definitions of Gfj(t2,ti), now they have an implicit dependence on the 
initial local excitation. The probability amplitude of finding a particle in site i at time 
t 2 when initially was in site j at time t% is described by the retarded Green's function 
of the whole system 

G R (t 2 ,h) = 9{hM) [G% (f 2 ,*i) ~ Gfj (h,^)} 

= MdtiM* ■ ( 3 - 15 ) 

The reduced density function G < (t,t), where matrix indices are restricted to i,j G 
{0,1}, is equivalent to the single particle 2x2 density matrix and G R (t2,^i) is an 
effective evolution operatoiQ If the system is isolated, the Green's function in its energy 
representation is obtained by a Fourier transform (FT) respect to the time interval 
St = t 2 - ti 

G 0R (e,t) = J G 0R (t + \5t,t- ISt) exp{is5t/h]d5t, (3.16) 

where t — ^ (t 2 + ti) . In a time independent system: 

G 0R (e,t) = G 0R (e) = [el - U s }-\ (3.17) 

After including SE interactions, the Green's function defines the reduced effective Hamil- 
tonian and the self-energy S R (^) [LPD90J, 

Heff.(e) = el - [G R (e)}' 1 = H s + S R (e). (3.18) 



The characters in bold are matrix representation of the respective operator. 
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Here, the exact perturbed dynamics is contained in the nonlinear dependence of the self- 
energy $] R on e. For infinite reservoirs the evolution with H e s. is non-unitary, hence, the 
Green's function has poles at the "eigenenergies" , £„, that have imaginary components 
[DP90] . 

- 2 Im S R (e„) /h = l/r SE = 2Y^/h. (3.19) 

These account for the "decay rates" into collective SE eigenstates in agreement with 
a self-consistent Fermi Golden Rule (FGR) |FP06j . Similarly, ReS R (e i/ ) = Ree^ — e Q v 
represent the "shifts" of the system's eigenenergies e° v . 

The evolution of the density function for the reduced open system is described using 
the Keldysh formalism |Kel64t IDan84j . The density function in the Danielewicz form 
|Dan84] is 



G< (t 2 , h) = KG (t 2 , 0) G< (0, 0) G (0, h) 

rt 2 rti 



+ / / dt k dt l G H (t 2 ,t k )i: < (t k ,t l )G A (t l ,t 1 ). (3.20) 
'o Jo 

The first term is the "coherent" evolution while the second term contains "incoherent 
reinjections" through the self-energy function, 'S < . This compensates any leak from the 
coherent evolution reflected by the imaginary part of £ R (see [Pas92j). The key to solve 



eq. (3.20) is to build up an expression for the particle (hole) injection and retarded 

(3.21) 



self-energies, (ti,t 2 ) and 

s R (h,t 2 ) = 9{t 1 ,t 2 ) [s > (t 2 ,t 1 ) -s < (t 2 ,t 1 )]. 



For this p urpose, we use a perturbative expansio n on 7^se f° r the Coulomb interaction 
DALP07 and the hopping interaction DPA 05 . The first order in the perturbation 
expansion is the standard Hartree-Fock energy correction which does not contribute to 
S < because it is real. We focus on the second order terms, with Feynman diagrams 



sketched in fig. 3.2 b). 

The injection self energy is 



sg (t k ,t l ) = \U la 



2 h 2 Gf 2 {t k ,ti) G% (t h t k ) G§ (t k , t{) 5^ + \V y 



12 



\ 2 Gf 2 (t k ,ti) 5 iX 8 ip 



where U\ 2 



-2U 



(dir.) 
12 



(3.22) 

^(exch.) . g ^ e Coulomb interaction between an electron in 



the system and one in the reservoir. The direct term contributes with a fermion loop 
and an extra spin summation which is represented in the —2 factor |Dan84j . The first 



term in eq. (3.22) corresponds to the direct and exchange self-energy diagrams shown in 



the last line of fig. 3.2 b). The first two diagrams schematize the creation of an electron- 
hole pair in the environment and its later destruction. The last term in eq. (3.22) and 



the last diagram of the same figure is the hopping to site 2 which allows the electron 
to perform a full exploration inside the reservoir. To take into account the different 
time scales for the dynamics of excitations in the system and in the reservoir, we use 
the time-energy variables |Pas92] : the physical time U = \ (t k + 1{) ,and the domain of 
quantum correlations 5ti = t k — U. This last is related to an energy e through a FT 
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|Pas92] . Thus, in equilibrium, 



G<{e, t i ) = i2nN 2 {e) f 2 (Mi) , 
G 2 > 2 (e,t i ) = -i27r N 2 {e) [1 - f 2 (M)] 



(3.23) 
(3.24) 



where N 2 (e) is the local density of states (LDoS) at the surface of the reservoir. Assum- 
ing that the environment stays in the thermodynamical equilibrium and k^T is much 
higher that any energy scale in the bath (high temperature limit), the occupation factor 
is 

f 2 Mi)=f 2 . 
Fourier transforming on e one obtains 



G 22 \t'i 



2 ' 



Sti 



) = i2vr g 2 (Sti) f 2 



and 
where 



G 22 (ti 



5ti 
2 i 



92 (Sk 



k - f ) = -i27r g 2 (Sti) (1 - f 2 ) 
de 



N, e) e 



2nh' 



(3.25) 

(3.26) 
(3.27) 
(3.28) 



Replacing in eq. (3.22) 



(*i + f ,*i - f ) = \U 12 \ 2 h 2 (2tt) 2 [g 2 (5ti)fi 2 [1 - f 2 ] Gfi (^ + f,ti~ f ) SaSy 



±\V l2 \ z i2ng 2 (Sti) 



1-f, 



My. (3.29) 



where the ( ^ 



associates f 2 with E < and (1 — f 2 ) with S > . 
In summary, we are left with the task to evaluate the time dependent self energies 
and the integral in eq. (3.20). We will focus in the parametric regime corresponding to 
the experimental conditions of the spin swapping gate (cross-polarization) described in 



section 5 2.1.2 



3.2.2.2 An environment in the wide band or fast fluctuation regime 

As occurs with the generalized Landauer-Biittiker equations for linear transport, an 
essential ingredient is the possibility to assign a Markovian nature to the environment. 
We are going to see that this appears naturally from the formalism when the dynamics of 
excitations within the environment is faster than the time scales relevant to the system. 
In order to separate the different physical time scales involved in the problem, we start 



changing to the time-energy variables in eq. (3.20). Evaluating in t 2 = t\ = t, the 
integrand becomes 



fdti f d^G R (Mi + f ) £< {ti + fA- f ) G A (ti - f,t) 
Jo J-t 



(3.30) 



The environment unperturbed Green's function g 2 (Sti) decays within the time scale 
h/Vs where Vb is the characteristic interaction inside the reservoir. In the wide band 
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regime (Vb 3> Vqi) H/Vb becomes much shorter than the characteristic evolution time of 
Crfi (ti + ^,t\— given by h/V m . Then, as the main contribution to the integral on 
Sti of eq. (3.20) is around the time scale H/Vb we can replace Gf x (t ; + ^,ti — ¥f) by 
G% Following the same assumption we replace G R (t,ti + ¥f\ by G R {t, t ; ) and 

G A (ti - ^,t) by G A (t h t) . In this fast fluctuation regime, only Eg (ij + ^,t; - ^) 
depends on <5tj leading to 



|t/ 12 |V (2tt 



[<?2 (^i)] 2 d^ 



± |V 12 | 2 i2vr 



f 2 [l-f 2 ]Gf x (Mi) S a S 



1.7 



<? 2 {SU) dSti 



f 2 
l-fa 



Mu, (3-31) 



which is local in space and time. This assumption for the time scales can be seen in fig. 



3.2 b) as a collapse of a pair of black dots, along a vertical line, into a single point. This 



justifies the expansion of fig. 3.2 c) and the use of the ladder approximation containing 
only vertical interaction lines in fig. 3.2 d). 

Assuming Ei = for % = 1, .., oo we obtain for the decay rates 



1 

T"SE 



i (Eft - EfJ 



(En E-q, 



\U 



121 



(27T) 



[^2 mfdsu 



f a [i-f 2 ] + iiv 12l 



2tt 



I 0V + r v ) 



where we have used t H/Vb to define 



h\u 



121 



2vr 2 



[<? 2 ((Jti)] 2 d5ti 



f 2 [1 - f 2 ] 



the Coulomb decay rate, and 



1^12 1 n 



92 {SU) dSU 



(3.32) 



(3.33) 



(3.34) 



the hopping decay rate. If one assumes that the environment ( 3.10) can be represented 
by a linear chain with all the hoppings equal to Vb the LDoS is DPA05 : 



N 2 {e) = l/(nV B )Jl 



2Vr 



(3.35) 



Thus, the Green's function 



92 {Sk) 



2ttVu 



SU 



(3.36) 
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is proportional to the first order Bessel function and decays within a characteristic time 
h/Vs- Assuming that f2 = 1/2 and the integration limits in the T's expressions are taken 
to infinity because t ~ fr/Voi 3> h/V E (wide band approximation), one obtains 

2 



r^xl^l 2 ^ (3-37) 



and 



^ 2 | 2 ^-. (3.38) 



It is important to remark that the wide band limit can be relaxed because the FGR holds 
when time t is in the range [FP06J t B < t < t R ~ ah/T SE ln(VB/r SE ), where a depends 
on the van Hove singularities of the spectral density Ji{e) = / N 2 (e — e')N 2 (e')de' . Here, 
£b = ^1^2(0) — h/Vs is the survival time of an electron-hole excitation at the surface site 
and £r characterizes the transition to a power law decay arising from memory effects. 
Hence, as long as Tse, Voi <C Vb, the FGR is valid for times much longer than H/Tse- 
Since the interaction is local in time, the reduced density results: 

G < (t,t) = /i 2 G R (t,0)G < (0,0)G A (0,t) + f dt i G R (t,t i )S < (t i )G A (t i ,t), (3.39) 

Jo 

which is complemented with 

S< = ( 2Y U UG< 1 (t u t-) + 2T v h (if 2 ) ) • (3 - 40) 

Here, the propagators G R (t, 0) and G A (0, t) that enter in both terms are obtained from 
the effective Hamiltonian of the reduced system, 

H-. - ( \ It", ) • (3.4D 

where Tse is energy independent. This results in an equation of the form of the GLBE. 
However, the Hamiltonian is asymmetric in the SE interaction complicating the form 
of the associated propagator. The apparent complexity to solve this equation contrasts 
with the homogeneous case where the evolution of the GLBE was obtained |Pas91] 
through a Laplace transform. Our strategy will be to induce such form of symmetry. 

3.2.2.3 A fictitious homogeneous decay 



The main difficulty with the eq. (3.39) is that it involves multiple exponentials. In order 
to create propagators with an homogeneous decay, i.e. a single exponential factor, we in- 
troduce fictitious interactions, S R C , with the environment. The symmetric Hamiltonian 
becomes 

H sym . = H c ff. + S R C 

-Voi \ , ( -i^r SE 



(3.42) 



-Voi -ir SE J V qr SE 



-i^TsE —V12 
— V12 — i^TsE 



3.2 A two level system interacting with a particle reservoir 



63 



Here Sg c includes the fictitious interactions which, in the present case, produce a leak 
of probability in site at a rate Tse/^ while in site 1 inject probability at the same 
rate. Both states of H sym . interact with the environment independently with the same 
characteristic rate I^e/^. Note that this rate is half the real value. The propagators of 
eq. (3.20) have now a simple dependence on t as 



G R (t,0) = G OR (t,0)e~^ t/h , (3.43) 
where 



«(^0) = G?f(t,0) = ^cos(^) (3.44) 



and 



G^(t,0) = O^0r = ^seu(^) (3.45) 
are the isolated system propagators. The reduced density evolution is now, 

G< (t, t) = h 2 G 0R (t, 0) G< (0, 0) G 0A (0, t) e-* /2TSE 

+ / dt.G * (t, ti) S< ym . (ti) G 0A (t u t) e-W***, (3.46) 
J o 

which is similar to the GLBE |Pas91l lPas92] . It is easy to see that the introduction of 
negative/positive imaginary parts in the diagonal energies of the effective Hamiltonian 
produces a decay/growth rates of the elements of the density function which, being 
fictitious, must be compensated by a fictitious injection self-energy 

£fic« = -him + SLj,) Gg (ti, ti) • (3.47) 

In our case, this results in an injection that includes the compensation effects for the 
symmetrized interaction, 

s< m . (to = s< (t ; ) + s< . (to 

o o \,(r m hG&(ti,t{) o 

2T u hG< 1 (ti, tO + 2T v h (if 2 ) J + \ -T SE hG< {t u ti) 

(3.48) 

Here, the second term is proportional to the local density functions Gf t (ti,tj) injecting 
and extracting density on sites and 1 respectively to restore the real occupation. We 
can rewrite the last expression to separate the processes that involve density relaxation 
(through injection and escape processes) and pure decoherence (through local energy 
fluctuations): 

S< m . (tO = £< (ti) + £< (ti) 



ir 



SE 



(f 2 -jG£(ti,tO) J V f^(ti,ti) 



(3.49) 
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Here 

jG< fo, ti) = j J G< (s, ti ) £L = f 1 ( ti ) (3.50) 

and 

jG<(t i ,t i ) = f Q (t i ), (3.51) 

while 

Pv = r y /r SE (3.52) 

is the weight of the tunneling rate relative to the total SE interaction rate. As the initial 
state has the site 2 occupied we have that 

-G<.(0,0) =5 ll 5 lr (3.53) 



Introducing eq. (3.49) into eq. (3.46) and using 

1 



— = lr SE (3.54) 

T"SE 



we get two coupled equations for G^ and Gf x 
*G<(t,t) = \hG°f(t,0)\\- t /^) + 



+ j |«3ff(t,t 1 )re-^^^[fG^(t ll t 1 )]. (3.55) 

fG<(t,t) = |/iG^(t,0)|V* /(2rsB) + 

j \hG<(f (t, ti )\ 2 e-^™ 2pv _^A_ [k _ k g < {tiA 



+ J I^^Ore-^^^IfG^^O]- (3-56) 

In each equation, the first term is the probability that a particle initially at site 1 be 
found in site (or 1) at time t having survived the interactions with the environment 
with a probability e~ t ^ 2l ~ SE \ The second term describes the process of injection/escape 
of particles enabled by the hopping from/towards the reservoir, where the last of such 
processes occurred in the time range (U,ti + dti) with a probability 2py^:- The injec- 
tion/escape is produced on site 1 and fill/empty the site to level it to the occupation 
factor f 2 . The third and fourth terms take into account the last process of measurement 
at time £; due to the SE interaction with a probability This confirms our inter- 

pretation that in eq. (3.49) the dissipation processes are in (t) while (i) involves 
pure decoherence. It is clear that by iterating this formula, one gets a series in the form 
represented in fig. 3.2 d). 
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3.2.2.4 The dynamics of a swapping gate 



The solution of the coupled eqs. (3.55) and (3.56) involves a Laplace transform. We 



consider a parameter range compatible with the spin problem where t% < 1 while we 
allow the tunneling relative weight pv in the range [0,1]. In a compact notation, the 
density function results: 



Poi (*) = f G< (t, t) = l- a e- Rot - ai cos [{u + i V )t + (p ] e~ Rlt . 



(3.57) 



Here, the decay rates Rq, R\ and r], and the oscillation frequency u are real numbers 
associated with poles of the Laplace transform. The amplitude ao is also real while, when 
u = 0, the amplitude a% and the initial phase ipo acquire an imaginary component that 
warrants a real density. These observables have expressions in terms of adimensional 
functions of the fundamental parameters in the model. Denoting 



X = UqTse 



(3.58) 



and remembering that 



pv = Ty/r 



SE, 



(3.59) 



we define 



and 



<t>{pv,x) 



1/2 2 l/i \2 

- I x -p v --(1-Pv) 



(3.60) 



X (pv, x) = {4(1- p v ) (9x 2 - 2 (1 - p v f + I8p 2 v ) 

+ 12 [3 (Ax 6 - ((1 - pvf + I2p v ) 



x 



+4p 2 v (5 (1 - pvf + 3p 2 v ) x 2 - 4p 2 v ((1 - pvf - P v )' 



(3.61) 



The observable "frequency" , 



VS fl 6{p v ,x) 

U + l7 l = TT I ftX \PV, X) + 6— r 

2x \6 x(Pv,aO 



(3.62) 



is purely real or imaginary, i.e. ur) = 0. Also, 



Rq 

Ri 



4>(pv,x) 1 / x , , 1 /, A 1 

c - ^X[Pv,x) +p v + - (1 -Pv) — : 

X{Pv,x) 6 3 J T SE 



(3.63) 
(3.64) 
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and 



a 



a-2 



a 3 



12(u 2 



2R\ 



jft) + (jR() _ Rl) 



2 ' 



1 (2^^ - uj 2 ) (R -Ri) + 2 (uj 2 - 7] 2 ) R 



2(lo + irj) (u 2 
1 uj 2 + 2R 2 - ARq^ 



r, 2 ) + (R - R x ) 



2 2 2 

a 1 = a 2 + a 3 , 



tan (d)n 



) 2 ' 

«2 

a 3 ' 



(3.65) 

(3.66) 

(3.67) 
(3.68) 



The oscillation frequency in eq. (3.62) has a critical point x r at a finite value o f x 
showing, as we called, a quantum dynamical phase transition ADLP06[ IDALP07 for 



which oo and rj in eq. (3.57) exchange their roles as being zero and having a finite value 



respectively. A full discussion of this issue for a spin system will be presented in chapter|4| 
Here, the dynamical behavior changes from a swapping phase to an overdamped phase. 
This last regime can be associated with the Quantum Zeno effect [MS77] where frequent 
projective measurements prevent the quantum evolution. Here, this is a dynamical effect 
[PU981 IPN94] produced by interactions with the environment that freeze the system 
oscillati on. 

Fig. 3.3 shows typical curves of f Gqo m the swapping phase. The different colors 
correspond to different SE interactions rates, pv = 0, 0.5 and 1, which are Coulomb 
(Ty = 0), isotropic (Ty = IV) and pure tunneling (IV = 0) interactions rates. The 
hopping interaction does not conserve the net energy in the system inducing a dissipation 
which is manifested through the non conservation of the number of particles in the 
system. This is the case of py ^ where the final state of the system has the occupation 
probability of the sites equilibrated with the bath occupation (f2). In fig. 3.3, this is 
manifested as the asymptotic normalized density (occupation probability) of 1. However, 
if py = 0, tunneling is forbidden and the system goes to an internal quasi-equilibrium 
as we described in the eq. (2.63), i.e. the local excitation is spread inside the system. 



In this case the asymptotic occupation probability of site is 1/2. 



3.3 Stroboscopic process 



Eq. ( 3.39[ ) has two main difficulties for a numerical implementation: The first is the 



evaluation of the system non-unitary propagators under inhomogeneous perturbations. 
The second is to keep track of all previous states of the system to enable the integration 
over previous times. We will show that the decay homogenization enables the implemen- 



tation of an efficient numerical algorithm. First of all, we identify in expression (3.46) 
that e~ t ^ 2rsE ^ = s(t) is the system's survival probability to the environment interrup- 
tion, i.e., the probability that the system remains coherent, and dtj/(2rsE) = Q (t{) dti is 
the "interruption" probability in a differential time around t{. The interaction of the en- 
vironment is discretized in intervals r str . where it acts instantaneously. This stroboscopic 
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Figure 3.3: Occupation probability, Poi (t) =^Gq Q (t) /h, to find a particle at site when 
initially was at site 1. Each line corresponds to different kinds, pv, of SE interactions. 
The plots correspond to x = VqiTse/^ = 10 belonging to the swapping phase and f 2 = 1. 



interaction leads to: 



s(t) = (l-p) n ^ 

oo 

<7(*) = X><J(*- 



rar, 



str.; 



(3.69) 
(3.70) 



m=l 



where 

n it) = int (t/r str .) . (3.71) 

Here, the stroboscopic interruptions may occur at the discrete times mr str . with a prob- 
ability p. At time t there were n (t) possible interruptions. In the joint limit r str . — > 
and p — > such that 

p/r str . = l/(2r SE ), (3.72) 

we can recover the continuous expression. To do this, notice that if n (t) = n, one can 
write eq. (3.69) as 



(" T atr.) 
S (t) — (1 — P) T str- 



Tstr. 

2^ 



(nT str .)/T st r. 



(3.73) 
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If t = nr^r then 



*(*)= ! 



7"str. 



2r< 



SE 



t/r str . 



(3.74) 



By taking the limit r str . — > the variable t becomes continuous yielding 

t/r s tr. 



s (t) = lim 1 



r str. 



2r SE ) 
exp [-*/ (2r SE )] , 



T s tr.— *0 



(3.75) 



which is the continuous expression for s (t) 



Then, by substituting p = r str J(2t se ) in eq. (3.70) we have 

q{t) = Yj r str .5(t - mr str .). 

2t S e , 
m=l 



(3.76) 



In the limit r str — > 0, t m — mr str , becomes a continuous variable and we can convert the 

:adh 



sum into an integral, leading to continuous expression of the GLBE (3.46): 

dt m 1 



2t S e Jo 



T s t T 8(t — t m ) 



Tstr. 2r S E 



(3.77) 



Coming back to the discrete version, we introduce eqs. (3.69) and (3.70) into the 
reduced density expression (3.46) to obtain 

G< (t, t) = h 2 G 0R (t, 0) G< (0, 0) G 0A (0, t) (1 - p) n{t) 

+ / dt i rsE^5(t i -t m )G OR (t,t i )S< m .(tOG OA (t i ,t)p(l-p) n(t -* i) . (3.78) 
Jo i 



m=l 



After integration, we obtain 
G< (t, t) = h 2 G 0R (t, 0) G< (0, 0) G 0A (0, t) (1 - p) n 

n 

+ h 2 Y, G° R (*, tm) SGf ni . (t m , t m ) G 0A (t, t m ) p (1 - p) n - m , (3.79) 



m=l 



where n = n(t) , t m = mr str and 

<5G< nj . (t,t) 



M M - 27 " SE y< 

^2 sym 



^sym. (0 • 



(3.80) 



G OR (t,0) = -^exp[-iH s t/fr] 



In this picture, the evolution between interruptions is governed by the system's propa- 
gators 

(3.81) 
(3.82) 



h 



and 



G 0A (0,t) = G UK (t,0) T . 



OR ( + n ^t 
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The spin-bath stroboscopically interrupts the system evolution producing the decay of 
the coherent beam. This decay is compensated through the reinjection of probabil- 
ity (or eventually of coherences) expressed in the instantaneous interruption function, 
(t,t), which also contains actual injection/decay from/to the bath. 
The first term in the rhs of eq. ( |3.79[ ) is the coherent system evolution weighted 



by its survival probability (1 — p) n . This is the upper branch in fig. 3.4 The second 
term is the incoherent evolution involving all the decoherent branches. The m-th term 
in the sum represents the evolution that had its last interruption at mr stL and since 
then survived coherently until nr str . Each of these terms is represented in fig. |3.4|by all 



the branches with an interrupted state (gray dot, red online) at the hierarchy level m 
after which they survive without further interruptions until fir str .. This representation 
has an immediate resemblance to that introduced by Pascazio and Namiki to justify the 
dynamical Zeno effect |PN94j . 



As mentioned above, the solutions of eqs. (3.79) and (3.46) are both computation- 



ally demanding since they involve the storage of all the previous steps and reiterated 
summations. Thus, taking advantage of the self-similarity of the hierarchy levels in the 



interaction with the environment, we rearrange expression (3.79) into a form optimized 
for numerical computation, 



h G < (t n+ i, t n+ i) — G° (i n +ij t n ) G < [t n , t n ) G° [t n , t n+ \) (1 — p) 



K 2 



+ G 0R (t n+1 ,t n ) (t n , t n ) G UA (t n , t n+1 )p. (3.83) 



OA 



in.] 



This equation provides a new computational procedure that only requires the storage 
of the density function at a single previous step. Besides, it avoids random averages 
required in models that include decoherence through stochastic or kicked-like perturba- 
tions |TFL+03[lDCM92] . 



3.3.1 A nice physical interpretation: The environment as a 
measurement apparatus 

We introduce our computational procedure operationally for an Coulomb interaction 
form (Vi2/t/i2 = 0) of 7^se- The initial state of the isolated two-level system evolves 
with Tis- At r str ., the particle reservoir interacts instantaneously with the "system" 
interrupting it with a probability p. The actual physical time for the SE interaction is 
then obtained as tse = T str /p. Considering that the dynamical time scale of the bath 
(tb — fr/Vs) is much faster than that of the system (fast fluctuation approximation), the 
dynamics of site 2 produces an energy fluctuation on site 1 that destroys the coherence 
of the two-level "system". This represents the "measurement" process that collapses 
the "system" state. At time r str ., the "system" evolution splits into three alternatives: 
with probability 1 — p the state survives the interruption and continues its undisturbed 
evolution, while with probability p the system is effectively interrupted and its evolution 
starts again from each of the two eigenstates of cfc 1 . These three possible states at r str 
evolve freely until the system is monitored again at time 2r str . and a new branching of 



alternatives is produced as represented in the scheme of fig. 3.5 a). 
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Figure 3.4: Quantum branching sequence for the stroboscopic evolution. Red dots 
represent states with interrupted (incoherent) evolution while the black dots are coherent 
with their predecessor. The horizontal continuous arrows represent the isolated evolution 
and the vertical dashed lines are the instantaneous interruptions. Notice the self-similar 
structure. 
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Figure 3.5: Quantum branching sequence for the swapping dynamics. Panel a) stands 
for a Coulomb system-environment interaction and b) a pure hopping one. Single states 
represent states with interrupted evolution (incoherent) while pairs of states are coherent 
superpositions. Notice the self-similar structure. 



A similar reasoning holds whe n V\% ^ 0. The sequence for isotropic interaction rates 
(IV = IV) is shown in fig. 3.5 b). The hopping part of 7^se can inject a particle. 
When an interruption occurs, the bath "measures" at site 2 and, if found empty, it 
injects a particle. The inverse process is not considered because we are assuming that 
the occupation factor in the reservoir is over 1/2 and we are normalizing the occupation 
probability to this value. In the figure, this can be interpreted as a "pruning" of some 
incoherent branches increasing the global coherence. In the next section, we will map 
the fermion system to a spin system to describe spin dynamics. Thus, this "pruning" 
explains why the rate of decoherence is greater when the Ising (Coulomb) part of 7Yse 
dominates over the XY (hopping) part. This occurs with the dipolar interaction and, 
in less degree, with the isotropic one. This contrasts with a pure XY interaction where 
the survival of spin coherences is manifested by a "spin wave" behavior. [MBSHf2Z] The 
empty site of the bath is refilled and de-correlates in a time much shorter than the 
swapping between sites and 1 (fast fluctuation approximation). Consequently, the 
injection can only occur from the bath toward the system. 
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3.4 Application to spin systems through the Jordan- 
Wigner transformation 



We apply this procedure to the spin system described in section § 2.1.2 providing a first 



principle derivation of the phenomenological equations employed there. We consider 
a system with M = 2 spins 1/2 coupled to a spin environment with the following 
Hamiltonian 



H = Hs + He + Hse, 



(3.84) 



where the system Hamiltonian, Hs, is described by eq. (2.30) 



h s = hn z (s z + !*\ + 16 (s+ii + s-if 
= hQ z (s z + if) + b (s x r$ + s y if) . 



(3.85) 



The environment Hamiltonian is described by 



i>% 



(3.86) 



This equation, in contrast with eq. (2.32), does not contain the Ising term. Thus, we 



have here an XY interaction inside the spin-bath instead of the dipolar interaction of 



section § 2.1.2 The other difference is that the SE interaction described by 



H 



SE 



a n im + \b l2 it / 2 - + /f/ 2 4 



(3.87) 



has only one connection between the environment (I2) and the system (Ji) while the given 



by eq. (2.31) has many. The purpose of the last simplification is only for a simplified 



description, but it could be generalized. Remember that, this spin-spin interaction is 
Ising if 612/012 = 0, and XY, isotropic (Heisenberg) or the truncated dipolar (secular) 
if a\ilb\i = 0, 1, —2, respectively. 

In the previous sections of this chapter we have observed the similarity of the results 
with the obtained ones in section 8 12.1.21 Thus, to use the results obtained within 



Keldysh formalism we map the spin system into a fermionic system using the Jordan- 
Wigner Transformation (JWT) [LSM61J. 



/ = c+exp { 



(3- 
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The previous Hamiltonians become^] 

H s = hn z (c+c + c+c l - i) + \b (c+c x + c+c ) , (3.89) 
H E = E «?z (^q. - |1) + E % {cfc 3 + cp t ) , (3.90) 

i>2 *> 2 

j>i 

Use = a 12 (cfc 1 - |l) (c^c 2 - |i) + |fo 12 (c^c 2 + cjcj , (3.91) 

where the index represents the spin S. Here, the system interacts with the environment 
through the site 2 (the surface site of the bath). In the last Hamiltonians, the terms 
proportional to the identity do not contribute to the dynamics because they only change 
the total energy by a constant number. This Hamiltonian, as we described in section § 



2.1.2, is a standard cross-polarization experiment (swapping gate) in NMR [MKBE74J. 
In this experiment, the site S is a 13 C and the site I\ a 1 H while the environment is a 
1 H spin bath. The typical experimental Hartmann-Hahn condition [HH62J equals the 
values of the effective energies at the 13 C and the X H sites to optimize the polarization 
transfer. The SE interaction has terms linear in the number operators cfc l and c\c 2l 
that only change the energy of the sites 1 and 2 respectively. Thus, the Hartmann- 
Hahn implementation, compensates the change of energy produced by the environment 

through these linear terms. Finally, we have Hamiltonians equivalent to those in eqs. 

b v — — h -ii- 7v( dir -) — n 

2' — 2 ' 12 — u 12 



(|3.9 3.10 3.111) where the site energies are equal, and V i = —I, = — U- ] 



and 0- 

The spin dynamics of the system is described by the spin correlation function |DPL04 
IDan06j : 

<*„.|/f(0)/f(0)|*„.)' 
which gives the local polarization at time t on the i-th spin with an initial local exci- 
tation on the 1-th spin at time t = 0. Here, |^/ e q.) is the thermodynamical many-body 
equilibrium state and 

I*(t) = eiwt/R/* e -i«t/R ( 3 93) 

are the spin operators in the Heisenberg representation. After the JWT, the initial local 
excitation on site 1 is described by the non-equilibrium state 

|*n.e.) = C+ |* eq .) . (3.94) 

In the experimental high temperature regime, k&T much larger than any energy scale 
of the system, the spin correlation function becomes [Dan06j 

Pa(t) = fG<(t,t)-l. (3.95) 

Notice that G^(t,t) implicitly depends on the initial local excitation at site 1. Here, 
Gfi (t, t) is the reduced density function of sites and 1 and can be split into the contri- 
butions G u {N) fa, ti) from each subspace with iV particles (or equivalently iV spins up) 



2 Note that the Jordan- Wigner transformation maps a linear many-body XY spin Hamiltonian into a 
system of non-interacting fermions. This leads us to solve a one-body problem, reducing the dimension 
of the Hubert space from 2 N to N states that represent local excitations LSM61J. 
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in the following way, 



M 



/M-1X 
I N-l J 



(3.96) 



N=l 



and analogous for the hole density function, 
described by 

^(o,o) = i(S 



The initial condition in this picture is 

■■ + M S * S u) > ( 3 - 97 ) 
where the first term is the equilibrium density (identical occupation for all the sites) and 



ij 



the second term is the non- equili 



we have an expression like (3.20 



system the —1 term of eq. (3.95 



Drium contribution where only site 1 is excited. Thus 



for each iV-th subspace |DPA0 5 . For this two-spin 



is canceled out by the background evolution, i.e. the 



evolution of the first term of eq. (3.97) plus the evolution of the second term of eq. (3.20 ) 
for the N = 2 subspace. As a consequence, the observable dynamics only depends on 
the initial local excitation at site 1, 



Gj w (0,0) = iMy, 



and evolves in the 1-th particle subspace, 

P«i(f) = ?G* w (M)- 



(3.98) 



(3.99) 



Finally, the solution of the polarization Poi(t), with T X y ^ and r zz <-> IV, is the 



same that was obtained in eq. (3.57) 



3.4.1 



Keldysh formalism versus the generalized quantum mas- 
ter equation 



Typical solutions of the quantum master equation for a spin swapping |CAL + 03[ IALP07] , 
described in chapter |2j were obtained following that of Miiller et al. [MKBE74] . They 
considered an isotropic interaction with the spin environment, represented by a phe- 
nomenological relaxation rate R = Txy — r zz — 1/ (4tse) • Within the fast fluctuation 
approximation and neglecting non-secular terms, this leads to 



exp [—Rt] — - cos(u;t) exp 



-Rt 



(3.100) 



that is obtained from eq. (2.58) in the isotropic case. This expression is used in most 



of the experimental fittings [PR96~l ILUP98] . Our eq. (3.57) reproduces this result with 
l/r^ = R ~ 1/ (4r SE ) by considering an isotropic relation between r zz and Txy, i-e. 
Pxy = 1/2 under the condition 1/ (4tse) b/h. Note that our microscopic derivation, 
for an XY chain as environment, does not imply that T zz = Txy comes from an isotropic 



SE interaction. This conclusion arrives from the expressions (3.37) and (3.38). 



Returning to the comparison between expressions (3.100) and (3.57) at short times 



t tse, one can see that the MKBE swapping probability growths exponentially with a 
rate 1/ (4tse)- m contrast, our solution manifests t hat the p olarizatio n gro wths quadrat- 



ically on time, {\b/ft) t 2 , as we anticipated in § 2.1.2.2| In eq. (3.57) this is made 
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possible by the phase 0o in the cosine. In the parametric region, 6tse <C h, where 
MKBE is not valid, our model enables the manifestation of the quantum Zeno effect 
|MS77t IPU98t |FP02| . This means that the bath interrupts the system through measure- 
ments too frequently, freezing its evolution. Here, this is a dynamical effect |PU98t lPN94j 
produced by interactions with the environment that freezes the system oscillation. At 
longer times, t 3> tse, one gets 

1 - P i (t) ~ (1 + 2 (b/h)\i E ) exp [-(b/h) 2 T SE t] , (3.101) 

and the quantum Zeno effect is manifested in the reduction of the decay rate oc 
ip/h) 2 tse as tse gets smaller than h/b. This surprising dependence deserves some inter- 
pretation. First, we notice that a strong interaction with the bath makes the X H spin to 
fluctuate, according to the Fermi golden rule, at a rate l/r SE . The effect on the 13 C is 
again estimated in a fast fluctuation approximation as 

l/r^(b/h) 2 r SE oc {b/hf [(alJtf + blJh 2 )^]' 1 . (3.102) 

This "nesting" of two Fermi golden rule rates is forma lly obtained from a continuous 



fraction evaluation of the self-energies jLPD9Q[ IDPA05 involving an infinite order per- 
turbation theory. Another relevant result is that the frequency depends not only on b, 
but also on t^e- A remarkable difference between the quantum master equation and 
our formulation concerns the final state. In the quantum master equation must be 
hinted beforehand, while here it is reached dynamically from the interaction with the 
spin-bath. Here, the reduced density, whose trace gives the system polarization, can 
fluctuate until it reaches its equilibrium value. 



We mentioned in section § 2.1.2.2 that by including the non-secular terms of the 



SE interaction in the generalized QME, we obtain the expression (3.57) for the polar- 
ization evolution. However, here we obtained the solution from a microscopic model 
using an XY linear chain as the spin environment, without forcing the final state. The 
Keldysh formalism gives us another perspective to discuss about the physics of quantum 
evolutions and the interpretation of the approximations made. 

It is important to remark that within the fast fluctuation approximation, the tem- 



poral and parametric behavior of expression (3.57) does not depend on the environment 
Hamiltonian. This only appears as an intensity in the values of Tzz and Txy and the 
relation between them (the anisotropy). 



3.5 Memory effects of the spin-bath 

The 13 C polarization, Poi(t), in the Keldysh formalism arises from the coherent evolution 
of the initial particle density, for which the environment is a "sink" , and an incoherent 
contribution where the bath acts as a particle "source" . This can be compared with the 
complementary framework. Instead of dealing with a "particle" problem we consider it 
as a "hole" problem [fig. |3.6|b) and c) respectively]. On these grounds, at t — 0, all the 



sites are occupied except for the "hole" excitation at the 0-th site. See fig. 3.6 c) where 



the black color stands for the hole excitation. At later times this excitation evolves in 
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Figure 3.6: a) Schematic representation of the spin system at time t — 0. The shaded 
region stands for the thermodynamic equilibrium state at high temperature and establish 
a background probability level. The black filling represents the excess of probability over 
the equilibrium state which is responsible for the observed dynamics. In b) the same 
system as in a) after the JWT, that is, under the particle point of view. Note that in 
this situation the background contribution is removed and the dynamics is described 
by the excess of probability AP. In c) we represent the complementary problem of the 
case b). Here the black filling stands for the hole that represent the excitation. In this 
representation it is easier to calculate the memory effects in the bath. 



the system and also propagates through the reservoir. The "environment" does not have 
holes to inject back into the "system" but those evolved coherently from the initial hole 
(i.e. S < = 0). Here the environment is a perfect "sink". Thus all the dynamics would 
be coherent, in the sense previously explained. If we add the result obtained in this case 
with that of eq. (3.57) we obtain a one for all times consequence of the particle-hole 
symmetry. This is a particularly good test of the consistency of the formalism because 
in each result the "environment" is set in a different framework. It also shows that the 
background polarization does not contribute to the dynamics as we obtained in (3.99). 

This "hole" picture can help us to get a very interesting insight on the dynamics in 
a case where the memory on the environment becomes relevant. Consider, for example, 
the case Voi = V12 = = V and E$ = E\. The finite version of this effective Hamil- 
tonian applies to the actual experiments reported in ref. |MBSH + 97j . In this case, the 
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simplifying approximations of the fast fluctuations regime are not justified. However, 
the exact dynamics of the system can be analytically obtained if one considers an infi- 



nite chain. This enables the use of eq. (3.35) to evaluate E , eq. (3.21). Then we can 



calculate the propagator through eq. (3.18) to solve eq. (3.20). The integration gives 
the first Bessel function, hence: 



Poi(t) = 1 



h_ 
tv 



j^tv/h) 



(3.103) 



A first observation is that the frequency above is roughly increased by a factor of two 
as compared with that in eq. (3.6). Since the maxima of Pqi are zeroth of the Bessel 



function it is clear that the frequency increases slightly with time. These are memory 
effects of the environment that are dependent on the interplay between the spectral 
density of the bath and that of the system. 

We notice that the memory effect can also appear in other condition for the bath. 
For example, if the proton nuclei have random polarizations and the density excitation 

., .. representing the 1 H sites filled up 



is at site 0, i.e. in fig. 



3.6 



| for n - 

to the shaded region; and the 13 C site with an occupation \ + AP. In this case the 
excitation propagates over a background level (shaded region) that does not contribute 
to the dynamics. The schematic view of this initial condition is equivalent to that of fig. 
c) where now the black filling represents a particle excitation. The solution of the 



3.6 



polarization is the first Bessel function, 



Poo(t) 



h 

tv 



j^tv/h) 



(3.104) 



as we obtain in (3.36). Apart from the finite size mesoscopic effect, this is precisely the 
situation observed in ref. [MBSH + 97j . although without enough resolution for a quanti- 
tative comparison. The effect of a progressive modification of the swapping frequency is 
often observed in many experimental situations such as CP experiments. Depending on 
the particular system, the swapping frequency can accelerate or slow down. Reported 
examples are fig. 5 on ref. |LUP98j and fig. 4 on ref. |LCP+04j . We can observe the 



same effect in the figure 2.10 for a three-spin system. This simple example solved so far 
shows that environmental correlations have fundamental importance in the dynamics 
and deserve further attention. 



3.6 Summary 



We have solved the Schrodinger equation within the Keldysh formalism for an open 
system within the w ide band regime (fast fluctua tion approximation) i n the environment 



DPA05 DALP07 



DALP07 ADLP07a that involves the 



We have shown a method 
transformation of the density function expressed in the Danielewicz integral form into 
a Generalized Landauer Biittiker Equation. This was possible by resorting to Wigner 
time-energy variables to perform the fast fluctuation approximation for the environment 
which leads to interactions local in time. This results in an injection of quantum waves 
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without definite phase relation with the initial state. The model proposed allowed us 
to consider the effect of the environment over the system via the decay of the initial 
state followed with an incoherent injection. Further on, we effectively symmetrized 
the syst em-environ ment interactions transforming them into a spatially homogeneous 



process 



ADLP07a 



This has an uniform system-environment interaction rate leading 
to a simple non-hermitian propagator. The original multi-exponential decay processes 
are recovered by an injection density function. Moreover, through discretization of the 
GLBE, we built a stroboscopic process which is the basis for an optimal numerical 



algorithm where the quantum dynamics is calculated in discrete time steps |ADLP06 



ADLP07a 



We applied |ADLP06[ lADLP07a] these techniques to a two-spin system coupled 



to a spin-bath improving the result obtained through the application of the secular 
approximation [MKBE74J in the standard density matrix calculation given in section § 
2.1.2 One improvement is the manifestation of the quantum Zeno effect that leads to 



novel interpretation of previous experiments |LUP98j which will be discussed in chapter 
4 The arising of the Zeno effect showed in the decoherence time, oc ib/ft) t~se, can 
Toe interpreted as a "nested" Fermi golden rule rate emphasizing the non-perturbative 
nature of the result 



ADLP06 



We observed that the solution of the Keldysh formalism 
is also obtained within the generalized QME if we include the non-secular terms of the 
SE interaction 



ALP07 



However, within the Keldysh formalism we derived it from a 
microscopic model for the entire system (system plus environment) and the final state 
must not be hinted beforehand. The Keldysh formalism gives us another perspective 
to discuss about the physics of the quantum evolution and the interpretation of the 
approximations made. Moreover, of particular interest, the Keldysh formalism allows 
us to include temporal correlations within the spin-bath in a model which has exact 
solution. On one side, it enabled us to show a novel res ult: memory effects can produce 
a progressive change of the swapping frequency |PPA05 . On the other side, these results 
will serve to test approximate methods developed to deal with complex correlations. 

In general, our analytical results based in the spin-particle mapping allow a deeper 
understanding of the polarization dynamics. They may constitute a starting point for 
the study of other problems. 



Chapter 4 



Environmentally induced Quantum 
Dynamical Phase Transition 



Experiments on quantum information processing involve atoms in optical traps |MKT + 00j . 
superconducting circuits |VAC + 02] and nuclear spins [BE Y+031 IPSM+03] among oth- 
ers. As we discussed in previous chapters, the system to be manipulated interacts 
with an environment |MKT+00l IGFMB03L IZur031 ISB04j that perturbs it, smoothly 
degrading its quantum dynamics with a decoherence rate, l/r^, proportional to the 
system-environment (SE) interaction H/tse- Strikingly, there are conditions where the 
decoherence rate can become perturbation independent |UPL98l lPLU + 00] . This phe- 
nomenon is interpreted |JP014 IJSB014 ICPJ04] as the onset of a Lyapunov phase, where 
1/-70 = min [l/r SE , A] is controlled by the system's own complexity A. Describing such 
a transition, requires expressing the observables (outputs) in terms of the controlled 
parameters and interactions (inputs) beyond the perturbation theory. We are going to 



show that this is also the case of the two-spin system treated in sections § 2.1.2 and 



3.4[ a simple swapping gate that is an essential building block for quantum information 
processing. While the swapping operation was recently addressed in the field of NMR in 
liquids [MBE98, LHBF99J with a focus on quantum computation, in the introduction we 
showed that the pioneer experiments were performed in solid state NMR by Miiller et al. 
[MKBE74J. They obtained a swapping frequency u determined by a two-spin dipolar 
interaction b, and a decoherence rate l/r^ = R that, in their model, was fixed by in- 
teractions with the environment 1/ (4tse)- This dynamical description was obtained by 



solving a generalized Liouville-von Neumann equation as in section § \2.1.2 As we antic- 



ipate previously, more recent experiments which span the internal interaction strength 



|LUP98j (see fig. |L6|) hinted that there is a critical value of this interaction when a dras- 
tic change in the behavior of the swapping frequency and relaxation rates occurs. Since 
this is not predicted by the standard approximations in the quantum master equation 
[MKBE74J, this motivated us to deepen into the physics of the phenomenon leading to 
develop the new theoretical result of chapter [3] 

In the first part of this chapter, we present a set of 13 C- 1 H cross-polarization NMR 
data, swept over a wide range of a control parameter (the ratio between internal inter- 
actions and SE interaction strengths). These results clearly show that the transition 
between the two expected dynamical regimes for the 13 C polarization, an oscillating 
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regime and an over-damped regime, is not a smooth cross-over. Indeed, it has the 
characteristics of critical phenomena where a divergence of the oscillation period at a 
given critical strength of the control parameter is indicative of the nonanalyticity of this 
observable [HL84, SacOlj . The data are interpreted by solving the swapping dynamics 
between two coupled spins (qubits) interacting with a spin bath. With this purpose we 



use the microscopic model proposed in section § |3.4| to describe the cross-polarization 
(swapping operation) using the Keldysh formalism. Within this picture, the overdamped 
regime arises because of the quantum Zeno effect [MS77t IPU9 8, FP02], i.e. the environ- 
ment "measures" the system so frequently that prevents its evolution. Such quantum 
freeze can arise as a pure dynamical process governed by strictly unitary evolutions 
|PN94l IPU98j .The analytical solution confirms that there is a critical value of the con- 
trol parameter where a bifurcation occurs. This is associated with the switch among 
dynamical regimes: the swapping phase and the Zeno phase. Consequently, we call this 
phenomenon a Quantum Dynamical Phase Transition. 

A major challenge in quantum control is the decoupling between the system and the 
environment. Many techniques |TED + 05[ |PJT + Q5l |MTA + Q6] are developed to avoid 
the loss of information induced by decoherence. We will show here how the criticality 
described above is useful to "isolate" a spin-pair. For this purpose, we will extend the 
model to a three interacting spin system where only one is coupled to the environment. 
We show that beyond a critical interaction, the two spins not directly coupled to the 
environment oscillate with their bare frequency and relax more slowly. In a two-spin 
system there is always a critical point that depends on the anisotropy relation of the 
SE interaction quantified as the ratio between the Ising and XY terms. However, in 
the three-spin system, the decoherence rate has a smooth cross-over from proportional 
to inversely proportional to the SE interaction. This cross-over approaches a critical 
transition as the anisotropy of the SE interaction goes from a purely XY to an Ising 
form. 



4.1 Experimental evidence 

The cross-polarization experiments exploit the fact that in polycrystalline ferrocene 



Fe(C5H 5 ) 2 (see fig. 4.1), one can select a pair of interacting spins, i.e. a 13 C and its 
directly bonded 1 H, arising on a molecule with a particular orientation. This is because 
the cyclopentadienyl rings perform fast thermal rotations (~ ps) around the five-fold 
symmetry axis, leading to a time averaged 13 C- 1 H interaction. The new dipolar con- 
stant depends [SH92] only on the angle 9 between the molecular axis and the external 
magnetic field Hq and the angles between the internuclear vectors and the rotating axis, 
which in this case are 90°. Thus, the effective coupling constant is 

1 /io 7H7cfi 2 (3cos 2 fl- 1) 

2 47rr| c 2 1 ' ' 

where 7's are the gyromagnetic factors and r HC the internuclear distance. Notice that 
biff) cancels out at the magic angle 9 m ~ 54.74°. As the chemical shift anisotropy of 13 C is 
also averaged by the rotation and also depends on 9 as (3 cos 2 9 — 1) , it is straightforward 
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Figure 4.1: Crystalline structure of ferrocene, Fe(C5H 5 ) 2 in its (room temperature) 
monociclic form with space group P2i/a. Two unit cells are included in the graph. 



to assign each frequency in the 13 C spectrum to a dipolar coupling b. Thus, all possible 
b values are present in a single polycrystalline spectrum. The swapping induced by b is 
turned on during the "contact time" t c , when the sample is irradiated with two radio 
frequencies fulfilling the Hartmann-Hahn condition [HH62, Sli92j. At t = 0, there is no 
polarization at 13 C while the 1 H system is polarized. The polarization is transferred forth 
and back in the 13 C- 1 H pairs while the other protons inject polarization into these pairs. 
We show the raw experimental data of 13 C polarization as a function of the contact time 



and b{6) in fig. 4.2 a). In fig. 4.3 a), the polarizations have been normalized to their 
respective values at the maximum contact time (3 ms) for each 9 when it saturates. It 
can be appreciated in the figure that the oscillation frequency is roughly proportional 
to |6| , showing that this is the dominant interaction in the dynamics. This is consistent 
with the fact that the next 13 C- 1 H coupling strength with a non-directly bonded proton 
is roughly 6/8 and, as all the intramolecular interactions, also scales with the angular 
factor (3 cos 2 9- 1). 

A noticeable feature in these experimental data is the presence of a "canyon", in 
the region \b\ < 2 kHz, where oscillations (swapping) disappear. The white hyperbolic 
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Figure 4.2: Raw experimental data of 13 C polarization as a function of the contact time 
and spin-spin coupling b(9) for a spin swapping dynamics in a 13 C- 1 H system. 



4.1 Experimental evidence 



83 




Figure 4.3: Spin swapping dynamics in 13 C- 1 H. a): Experimental 13 C polarization in 
Fe(C5H 5 ) 2 as a function of the contact time t c and spin-spin coupling b(Q). b): Numerical 
simulations of the 13 C polarization obtained from eq. (4.12) for different values of b, 
a dipolar system-environment interaction (| 012/^12 1 — 2) and a constant value for tse 
(tse = 0.275 ms) obtained by fitting the experimental data in the regime where the 
MKBE expression is valid. Projection plots in the b — t plane show a canyon where the 
oscillation period diverges indicating a Quantum Dynamical Phase Transition. 
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stripes in the contour plot at the bottom evidence a swapping period 2ti /uj that diverges 
for a non-zero critical interaction. This divergence is the signature of a critical behavior. 

The standard procedure to characterize the cross-polarization experiment in fer- 
rocene and similar compounds was described in section § 2.1 .2| that is derived from the 
MKBE model [MKBE74]. There the 13 C polarization exchanges with that of its directly 
bonded 1 H, which, in turn, interacts isotropically with other protons that constitute the 



environment. Their solution is (3.100) 



1 exp 



t 

T 4> 



cos (cut) exp 



3 t 
'2t> 



(4.2) 



where the decoherence rate becomes determined by the rate of interaction with the 
environment \T^/h =1/ (4tse) — ► l/ r <£ = R, while the swapping frequency is given 
by the two-spin dipolar interaction, b/h — > uj. A dependence of the inputs b and tse 
on 9 should manifest in the observables u and r^. However, working on a polycrystal, 
each tse(6 i ) value involves a cone of orientations of neighboring molecules and a rough 
description with single average value for the SE interaction rate is suitable. 

We have performed non-linear least square fittings of the experimental points to the 
equation P (t) for the whole 13 C spectra of ferrocene in steps of ~ 80 Hz and contact 
times ranging from 2 /is to 3 ms. The and uj parameters obtained from these fits are 
shown as dots in fig. 4.4 The proportionality of the frequency with b for orientations 



that are far from the magic angle is verified. In this region a weak variation of 
around 2.2 kHz reflects the fact that 1/ (2tse) does not depend on 6. A drawback of this 
simple characterization is that it tends to overestimate the width of the canyon because 



of limitations of the fitting procedure when eq. (4.2) is used around the magic angle 



In spite of the MKBE theoretical prediction, one observes that the frequency be- 
comes zero abruptly and the relaxation rate suddenly drops with a quadratic behavior 
when 6 C — 2 kHz. The minimum of the parabola occurs at the magic angle, when 6 = 0. 
Then, all the polarization reaching the 13 C at this orientation originates from protons 
outside the molecule. Then, the rate of 0.5 kHz obtained at this minimum constitutes 
an experimental estimation of this mechanism. This has to be compared with the al- 
most constant value of l/(4r SE ) = — 2.2 kHz observed outside the magic angle 
neighborhood. This justifies neglecting the J-coupling and the direct relaxation of the 
13 C polarization through the dipolar interaction with protons outside the molecule. In 
the following we describe our stroboscopic model that accounts for the "anomalous" 
experimental behavior. 



4.2 Theoretical description 
4.2.1 The system 

Let us consider M coupled 1/2 spins with a Hamiltonian: 



H = H Z + J2 <HjI*Ij + K [Itlj + lilt) / 2 
i<j 1 v 7 



(4.3) 
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Figure 4.4: Decoherence rate l/r^ and frequency uo in the spin swapping of a 13 C- 
1 H system. Data points are obtained by fitting cross polarization experiments to the 
expression P m ^(t). The zero plateau in the frequency and the parabolic behavior of 
I/t^ in the region 6tse <C h are indicative of an over-damped Zeno phase. Solid lines 
are the prediction of our model assuming a constant tse = 0.275 ms. 



where 

M 
i=l 

is the Zeeman energy, with a mean Larmor frequency Qz- As we mention in previ- 
ous chapters, the second term is the spin-spin interaction: b^/ciy = is Ising, and 
dij/bij = 0, 1, —2 gives an XY, an isotropic (Heisenberg) or the truncated dipolar (sec- 
ular), respectively. 

In order to describe the experimental system, we use the model described in sections 
Hence, let us take the first N = 2 spins, I = S (a 13 C) and I\ (its 



2.1.2 and S 3.4 



directly bonded H), as the "system" where the swapping ||,T) ^ IT; I) occurs under 
the action of 6qi- The other M — N spins (all the other 1 H), with M — > oo, are the spin- 
bath or "environment" . This limit enables the application of the Fermi Golden Rule or 
a more sophisticated procedure to obtain a meanlife tse for the system levels. We will 



not need much detail for the parameters of the spin-bath in eq. (4.3) except for stating 
that it is characterized by an energy scale d-Q which leads to a very short correlation time 
tb — h/d-Q. Thus, this spin system can be mapped into a fermion particle system using 
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the Jordan- Wigner transformation |LSM61] as we described in section § 3.4 Under the 
experimental conditions, SQi = 0, aoi = and 6qi — b, the system Hamiltonian becomes 



(3.89) 



H s = hn z (c+c + cfc 1 - 1) + \b (c+Cx + c+c ) . (4.5) 

The Jordan- Wigner transformation maps a linear many-body XY spin Hamiltonian 
into a system of non-interacting fermions. Spins li with 2 < i < M — 1 are interacting 
among them. To simplify the presentation, and without loss of generality, we consider a 
single connection between the system and the spin-bath ay = b\j = 0, j = 3 . . . oo and 

a j 



J 0j 



0, j = 2 . . . oo to use the SE interaction described by eq. (3.91 ) 
7~Lse = ai2 (cfc 1 — \ 



!)( 



c 2 C 2 



(4.6) 



In the experimental initial condition, all spins are polarized with the exception of S 
[MKBE74J (see section § 2.1.2). Thus in the high temperature limit {Mlz/k-QT e s < 1), 
the reduced density operator is 



<7(0) = ?G< (0)= (i + s/^/Trji} 



which under the Jordan- Wigner transformation becomes 



Tr{l} Tr{l} 



(4.7) 



(4- 



Since the first term does not contribute to the dynamics (§ |3.4 ), we retain only the 
second term and normalize it to the occupation factor. This means that site 1 is empty 
while site 2 and sites at the particle reservoir are "full". This describes the tiny excess 
Ap above the mean occupation 1/2. To find the dynamics of the reduced density matrix 
of the "system" 

a (t) = -G< (t) , (4.9) 
i 

we will take advantage of the particle representation and use the integral form |Kel64[ 



Dan84] of the Keldysh formalism, eq. (3.20), instead of the standard Liouville-von Neu- 
mann differential equation, eq. (2.39), as we described in section § 3.4. There, any 



perturbation term is accounted to infinite order ensuring the proper density normaliza- 
tion. The interaction with the bat h is local and, because of the fast dynamics in the bath, 



it can be taken as instantaneous DALP07 ADLP07a . Hence, the evolution is further 



simplified into an integral form of the G eneralized Landauer B iittiker Equation (GLBE) 



[Pas91[ IPas92j for the particle density DALP07[ lADLP07a . There, the environment 



plays the role of a local measurement apparatus § 3.3.1 



4.2.2 Analytical solution 



After all these assumption, we have the spin system used in section § |3.4| where we apply 
the spin-fermion mapping to use the Keldysh formalism. In that section, we obtain that 
the observable dynamics is consequence of a initial local excitation on site 1, 

GS (1) (0,0) = iMii, (4.10) 
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and evolve in the 1-st particle subspace, 

P*i(*) = ?GS w (M). (4-11) 
After that, the solution of the polarization Poi(t) (experimental 13 C polarization) is the 



same that was obtained in (3.57) 



P i(t) = l- a e 



-R t 



a\ cos [{uj + irj) t + 0o] e" 



-Rxt 



(4.12) 



where V m = -f, 



■|, Vij = — ^y, Utf U ' J = a 12 and ' — and we substitute in the 
present solution Txy ^ IV, Tzz <-> IV and pxy ^ Pv- In spite of appearance, the 
last equation has a single fitting parameter. This is because the real functions u, r], R 
and Ri as well as a , ai and 4> depend exclusively on b, r SE and pxy- Besides, 6 and 
Pxy are determined from crystallography and the anisotropy of the magnetic interaction 
(Pxy = 1/5 for dipolar) respectively. The phase transition is ensured by the condition 
u>r] = 0. The aniso tropy ratio IV y/IVe — > Pxy accounts for the observed competition 
(see section § 2.1.2 and § 2.2.2.1) between the Ising and XY terms of 7Yse- The Ising 
interaction drives the "system" to the internal quasi-equil ibrium (2.63 ). In contrast, the 
XY term allows the thermal equilibration with the bath |CAL + 03~f ! 



4.2.3 Comparison with the experiments 

In order to see how well our model reproduces the experimental behavior we plot the 13 C 
polarization with real istic parameters. Since the system is dominated by the dipolar SE 
interaction (CAL + 03j . as we described in chapter [j] we take |ai2/£»i2 1 = 2. We introduce 
b with its angular dependence according to eq. (4.1) and we select a constant value 
for tse = 0.275 ms representative of the b ^> H/tse regime. Since here, we are only 
interested in the qualitative aspects of the critical behavior of the dynamics, there is 
no need to introduce tse(#) as a fitting parameter. These evolutions are normalized at 



the maximum contact time (3 ms) experimentally acquired. They are shown in fig. 4.3 
b) where the qualitative agreement with the experimental observation of a canyon is 
evident. Notice that the experimental canyon is less deep than the theoretical one. This 
is due to intermolecular 13 C- 1 H couplings neglected in the model. We will show that 



the analytical expression of eq. (4.12) allows one to determine the edges of the canyon 



which are the critical points of what we will call a Quantum Dynamical Phase Transition 
(QDPT). 



4.3 Quantum Dynamical Phase Transition 

Our quantum observable (the local spin polarization) is a binary random variable. The 



dynamics of its ensemble average (swapping probability), as described by eq. (4.12), 
depends parametrically on the "noisy" fluctuations of the environment through tse- 
Thus, following Horsthemke and Lefever |HL84j . one can identify the precise value for tse 
where a qualitative change in the functional form of this probability occurs as the critical 
point of a phase transition. This is evidenced by the functional change (nonanalyticity) 
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of the dependence of the observables (e.g. the swapping frequency u>) on the control 
parameter br^/h. Since the control parameter switches among dynamical regimes we 
call this phenomenon a Quantum Dynamical Phase Transition. 

It should be remarked that the effect of other spins on the two spin system introduces 
non-commuting perturbing operators (symmetry breaking perturbations) which produce 
non-linear dependences of the observables. While this could account for cross-over among 
the limiting dynamical regimes, it does not ensure a phase transition. A true phase 
transition needs a non-analyticity in these functions which is only enabled by taking the 
thermodynamic limit of an infinite number of spins |Sac01] . In our formalism, this is 
incorporated through the imaginary part of the energy, H/tse, evaluated from the Fermi 



Golden Rule (3.19) 



When the SE interaction rates are anisotropic (T^z 7^ Txy), there is a functional 
dependence of u on tse and b yielding a critical value for their product, 6tse/^ = ^ XY , 
where the dynamical regime switches. One identifies two parametric regimes: 1- The 
swapping phase, which is a form of sub-damped dynamics, when br SE /k pxY > h (r/ = 
in eq. (4.12)). 2- A Zeno phase, with an over-damped dynamics for brs E /k pxY < h as a 
consequence of the strong coupling with the environment (zero frequency, i.e. us = 0, in 



eq. (4.12)). In the neighborhood of the critical point the swapping frequency takes the 
form: 



to 



PXY 







2 

SE 



bTsF,/k PxY > h 

frrSEAp XY < ft 



(4.13) 



and r] becomes 



PXY 



bTs E /k PXY > h 



fc p XY / T SE - ( h / h f br SE /k PxY < h 



(4.14) 



where 



h 2 

>XY 



12 



[(Pxy - I) 2 X (Pxy)] ~ 3 + C (pxy) + 19p XY + 



(PXY ~ 1) 3 C (PXY 

[X (Pxy)] 5 



(4.15) 



and 



X (Pxy) = -5291p| Y - 1084p| Y + 546p| Y - Ap XY + 1 + 2AV3 P xy \j (28p 2 XY - 2p XY + l) 3 , 

(4.16) 

C (Pxy) = -215pxv _ 2 Pxy + 1, 



1 (fx 2/3 + 36/ 2 ) (-fsA 2/3 + 36/ 2 / 3 + /i/ 4 ) 



PXY 



fi 5/3 u 



fx = 36^ _ 8 + 24p XY + 48p| Y - 36k 2 p v XY - 64p XY + 12 f 4 , 



h = 1/3^ XY - l/3p XY - 1/9 (1 - pxy) 2 , 

„2 1C „4 , nnJ mJ , i o„2 ;„2 o 1.2 „ , ;2 



h = -/4+PXY/4-6^ XY -16pxY + 20pxY-10pxY + 13PxY^ XY -2^»XY + fc 



Pxy' 



(4.17) 

(4.18) 

(4.19) 
(4.20) 
(4.21) 
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h = (l2k 2 pxY + 96p XY k 2 pxY - 120p XY k 2 pxy + m P 2 XY k 2 pxy - 39k 2 pxY p 2 XY 

+Qkl XY p X y - Sk 2 pxY - A8p XY + A8p 3 XY - 12p 2 XY ) 1 . 



In equation (4.18) the functions fx, f 2 , f% and / 4 only depend of pxy- Thus, the pa- 
rameters a PxY and fc PxY only depend on p XY which is determined by the origin of the 
interaction Hamiltonian. For typical interaction Hamiltonians the values of these param- 
eters (pxy, k PxY , a PxY ) are: Ising (0, |, l), dipolar (|, 0.3564, 0.8755), isotropic 
and XY (1,1,1). Fig. 4.5 shows the parameters a PxY and function of Pxy- 



1 n 1 

2' U ' V2 




Figure 4.5: a) Critical value of the swapping frequency as a function of p XY (anisotropy 
of the SE interaction), b) The parameter a PxY of the swapping frequency as a function 
of pxy- 



It is important to note that when the SE interaction rates are isotropic (r X Y = T zz ) 
the critical point is zero, thus, the oscillation frequency is finite for all values of the SE 
interaction rate Tse/^- While the oscillation could be very attenuated it is ensured a 
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swapping with a frequency proportional to b. The swapping period is 



T 



T, 



3/2 



0c 



V2i 



Oc 



-1/2 



where 



PXY 



Tn 



2irh 



is the isolated two spin period and its critical value 

2vrr SE 



Oc 



(4.22) 



(4.23) 



(4.24) 



determines the region where the period T diverges as is observed in fig. 4.3 The 
estimated value of tse = 0.275 ms and dipolar SE interactions yield a critical value for 
the ^C^H coupling of b c /h = 2vr /T 0c = 1.3 kHz . 

The complete phase diagram that accounts for the anisotropy of the SE interactions 



is shown in fig. 4.6 There, the frequency dependence on pxy and br^/h is displayed. 



At the critical line, kp xY as a function of Pxy, the frequency becomes zero setting the 



limits between both dynamical phases. 

The two dynamical phases can now be identified in the NMR experiments which up 
to date defied explanation. The experimental setup provides a full scan of the parameter 
6tse through the phase transition that is manifested when the frequency goes suddenly 



to zero (fig. 4.4 b) and the relaxation rate [fig. 4.4 a)] changes its behavior decaying 
abruptly. The fact that tends to zero when b <C Ji/tse confirms the Zeno phase 
predicted by our model. In this regime, l/r^ is quadratic on b as prescribed by eq. 
P1021 ), 

1/T0 oc (b/h) 2 r SE . 



(4.25) 



4.3 quantitative, we fit the 



To make the comparison betw een the two panels of fig. 

MKBE 

predicted dynamics of fig. 4.3 b) with P (t) , following the same procedure used to 
fit the experimental data. The solid line in fig. 4.4 show the fitting parameters l/r^ and 



uj in excellent agreement with the experimental ones. 



We point out that eq. (4.2) is used to fit both the experiments and the theoretical 



prediction of eq. (4.12 ) because it constitutes a simple, thought imperfect, way to extract 
the "outputs" (oscillation frequency uj and a decoherence time r^). While the systematic 
errors shift the actual critical value of the control parameter, b/h, from 1.3 kHz to 2 kHz, 
eq. (4.2) yields a simplified way to "observe" the transition. 



4.4 Signatures of a QDPT in a three-spin system 
coupled to a spin-bath 

We consider a three-spin system with the following Hamiltonian 

n s = m z (f, + i* + /*) + \b (/+ i + /z 1 / + ) + \b (i+r x + / -/+) , (4.26) 
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Figure 4.6: Quantum dynamical phase diagram for the spin swapping operation. The fig- 
ure shows the frequency dependence on system-environment (SE) interaction anisotropy 
Pxy and the ratio among the internal and the SE interaction br^/h. The projection 
over the br^/h — pxY plane determines the phase diagram where the transition between 
the swapping phase into the Zeno phase (u = 0) is manifested. Values of pxy for typical 
SE interaction Hamiltonian are indicated in the contour plot. 
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where we add one spin to the system. However, in this case there is not a interac- 
tion between the spin i_i and the I\ as there was in section § 2.2 The environment 



Hamiltonian and the SE interaction Hamiltonians remain as in the previous section, 



*>2 4 > 2 V 

j>i 



^SE = 12 ^2 + 2^12 ( Iih + h ^2 



(4.27) 
(4.28) 



Also, the environment is coupled only to one spin of the system. We solve the general- 



ized quantum master equation (2.39), using the same procedure as in section § 2.1.2.2 



without neglecting non-secular terms of the relaxation superoperator. Considering now 
the initial condition 

g(0) - Tr { i } ' ^ 

and the spin-bath polarized, we obtain for the magnetization of site —1 

Af/. 1 (t) = Tr{/V(t)} (4.30) 

= M (1 - a e~ Rot - aie" Rl * + a 2 sen (u 2 t + fa) e~ R2t + a 3 sen (uj 3 t + fa) e"^ 3 *) . 

(4.31) 



The coefficients a i; R iy u)i and fa are real and they are functions of b, l/r SE and pxy- If 
Pxy 7^ the final state has all the spins polarized because a net transfer of magnetization 
from the spin-bath is possible. However, for an Ising SE interaction, pxy = 0, we obtain 
that Rq = and 1 — do = 1/3 (the asymptotic polarization) because the final state is the 
quasi-equilibrium of the 3-spin system as we described in section § 2.2.2.1, eq. (2.116). 
In fig. 4.7 we show the frequencies u>2 and u 3 and the different relaxation rates as a 
function of ibr^/h)~ l when the SE interaction is Ising (pxy = 0). Two changes, each 
resembling the critical behavior shown by two spin systems are observed (see fig. 4.4). 



The same phenomenon occurs in figure 4.8 a) where the coefficients are shown. The 



polarization evolution of an isolated 3-spin system is 



M Iti (t) 



(3 + 4 cos (u%t) + cos (wgt)) 



where 



and 



V2b 



(4.32) 



(4.33) 



(4.34) 



are the natural frequencies. When {br^/K) 
ai -> 1/3 - 3/8 = -1/24, a 2 -> 1/2 and a 3 
spin dynamics. The dependence of u 3 as a function of 
the swapping frequency of fig. WA 



« 1, we observe that u 2 — >• uj%, u 3 



UJ 



3' 



1/8 as expected for an isolated 3- 
t se/^) _1 is similar to that of 
However, instead of becoming zero when the SE 
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Figure 4.7: a) Frequencies involved in the temporal evolution of the polarization in the 3- 
spin system as a function of (6r SE //i) _1 . Dashed lines represent the isolated system. Dot 
line correspond to two spins decoupled from the environment, b) Different relaxation 
rates of the polarization. 
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Figure 4.8: a) Coefficients (weights) of the different terms of eq. ( 4.30| ). At the critical 
region there is a switch between the 2-spin and the 3-spin regime, b) and c) Temporal 
evolutions of the polarization in the 2-spin and 3-spin regimes respectively for different 
Tse- In b) b/h = 2n x 1 kHz and Tse = 1.43 ms for the thick line and t$e = 10 ms for the 
thin line. In c) b/h = 2n x 1kHz and t$e = 0.1ms for the thick line and tse = 0.01ms 
for the thin line. 
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interaction increase, it suddenly stabilizes at ojq = b/h~, the bare 2-spin Rabi frequency. 
At the same point u>2, R2 and R3 also have a sudden change. While R2 and R3 initially 
grew as (&tse/^) _1 , there R2 increases the growing speed while R3 begins to decay as 
in the Zeno phase (see fig. 4.4). Moreover, looking at the behavior of the coefficients 
a 2 and a 3 , we observe a form of switch between them, a 2 suddenly goes down and a 3 
goes up. These coefficients are the weight of the different frequency contributions in 
the time evolution. Both changes together, on the decoherence rates and on the weight 
coefficients of the different oscillatory terms, beyond the critical interaction (region) leads 
the system to oscillate with the bare Rabi frequency of the two spins decoupled from the 
environment. If we continue increasing the control parameter ipT^/K)" 1 , this effect is 
enhanced by the next transition. After the second transition, R\ begins to decrease as in 
the Zeno phase behavior of fig. 4A As the term of eq. 4.30 that relax with R\ leads the 
system to the 3-spin quasi-equlibrium, when it goes down this final state is approached 
more slowly. This effect tries to avoid the interaction between the 2-spins not coupled 
to the environment and the I2 spin. After the second transition, the coefficient a% goes 
to zero abruptly leading to a more pronounced "isolation" of the 2-spins. After that, 
we can characterize two dynamical regimes: One which is characterized by the 3-spin 
dynamics when (6tse/^) 1 < 1 and the second one, when (&tse/^) 1 > 1, have the 
2-spin behavior. 

Fig. 4.8 b) and c) show the temporal evolution of the magnetization of eq. (4.30) on 
the 3-spins and the 2-spin regimes respectively. While in fig. 4.8 b) the two frequency 
contributions are evident, in fig. 4.8 c) only the bare Rabi frequency is manifested. 

In fig. 



In 

4.8 b), we show 



each graph we show two curves with different SE interactions 
that increasing the SE interaction the decoherence rates increase. However, in the 2-spin 
regime [fig. 4.8 c)] when it is increased the decoherence rate decrease leading to a better 
"isolation". It is important to take into account that while the relaxation rates goes 
to zero smoothly the swapping frequency acquire the bare value near the critical point. 
Another fact to remark is that this effect is more pronounced when the anisotropy of 
the SE interaction is close to a pure Ising SE interaction while an increase in the XY 
nature leads a further smoothing of the transition. The reason is that, when pxy 7^ 0, 
there is a net transfer of magnetization to the system which is redistributed between 
the three spins, this redistribution begins to be slower at the second transition when 
i?3 goes down. In contrast, for a purely Ising interaction, there is no net polarization 
transfer and a purely decoherent process at site 3 freezes its dynamics but its fast energy 
fluctuations prevent the interaction with the other spins. 



4.5 Summary 



We found experimental evidence that environmental interactions can drive a swapping 
gate th rough a Q uantum Dynamical Phase Transition towards an over-damped or Zeno 



phase ADLP06 . The NMR implementation of a spin swapping in a 13 C- 1 H system 
enables the identification and characterization of this phase transition as a function of 
the ratio br^/h bet ween the internal and SE interaction. In chapter [3l we developed 



a microscopic model |DALP07[ lADLP07aj that describes both phases and the critical 
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region with great det ail, showing that it depends only on the nature of the interaction 



ADLP06[ lADLP07a] . In particular, the phase transition does not occur if the SE inter- 
action gives isotropic interaction rates, Tzz = Txy- The phase transition is manifested 
not only in the observable swapping frequency but also in the decoherence rate 
While a perturbative estimation through the standard Fermi golden rule would tend to 
identify this rate with the SE interaction, i.e. l/r^ = l/r SE ~ i.d 23 /h) 2 tb, as it occurs 
well inside the swapping phase, both rates differ substantially as the system enters in 
the Zeno phase (br SE < k PxY h). Here, the decoherence rate switches to the behavior 
oc {b/h) 2 T SE . In the Zeno phase, the system's free evolution decays very fast with 
a rate r^. In spite of this, one can see that the initial state as a whole has a slow 
decay (its dynamics becomes almost frozen) because it is continuously fed by the envi- 
ronment. Since the r SE has become the correlation time for the spin directly coupled 
to the environment, 1/r, provided by our calculation can be interpreted as a "nested" 



Fermi golden rule rate emphasizing the non-perturbative nature of the result [ADLP06 



Based on the wealth of this simple swapping dynamics, we can foresee applications that 
range from tailoring the environments for a reduction of their decoherence on a given 
process to using the observed critical transition in frequency and decoherence rate as 
a tracer of the environment's nature. This led us to extend the model to a 3-spin sys- 
tem and to show that beyond a critical region the two spins become almost decoupled 
from the environment oscillating with the bare Rabi frequency and relaxing more slowly 



ALP07] . While in the two spin swapping gate the dynamical transition is critical, in 
the 3-spin system the criticality is smoothed out. However, enough abruptness remains 
to give the possibility to use it to "isolate" a two-spin system with a finite system- 
environment interaction. Thus, these applications open new opportunities for both, the 
field of quantum i nformation processing and the general physics and chemistry of open 
quantum systems 



DALP07 



Chapter 5 

Polarization transfer enhancement 
by selective pruning in NMR 



Inspired in the stroboscopic model discussed in section § |3.3[ we propose a new NMR 
pulse sequence to improve the transfer of polarization through a specific pathway in a 
system of many interacting spins. The sequence effectively prunes branches of spins, 
where no polarization is required, during the polarization transfer procedure. In this 
way, a local polarization excitation is transferred more efficiently from a chosen "source 
spin" to a "target spin" . Simulations of the spin dynamics with real values of chemical 
shifts and J-couplings in the 13 C backbone of leucine for a couple of sources and targets 
are performed. Enhancements of a factor of up to 300% in the target polarization 
are observed without sacrificing the relevant times. Possible applications and potential 
fundamental contributions to engineered decoherence are discussed. 



5.1 The pruner sequence 

One standard procedure to transfer polarization within molecules in the liquid state, 
called TOCSY (Total Correlation Spectroscopy) |BE83| IBD85] is based on the isotropic 
interaction (or J-interaction). In general, one has a mixing Hamiltonian which makes 
possible the spin-spin interaction. The mixing Hamiltonian in the rotating frame for an 
isotropic interaction is 

where is the scalar coupling and J" [u = x, y, z) are the spin operators. As this 
Hamiltonian connects all the possible pairs of spins, an initial local excitation would 
spread over all of them. The idea behind our new sequence is to effectively "disconnect" 
those spins which do not belong to the "selected transfer pathway" . In order to achieve 
this goal, one makes use of the differences in chemical shifts in the spin system. 

"The pruner" sequence stroboscopically interrupts the mixing evolution performed 
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under 7i m ° x with free evolutions, that in the rotating frame are given by 



i i^j 



3 ' 



(5.2) 



where Af2j are the off-resonance shifts as defined in 2.7 and A^ = h(ACli — AClj) ^> Jij. 
Thus, there is a time interval, At mix , with the mixing evolution under the Hamiltonian 
^mix followed by a free evolution time step At free . This sequence of evolutions are 



repeated successively as schematized in fig. 5.1 The adjustment of the free evolution 



The pruner sequence 



At 



mix 




Figure 5.1: Schematic representation of the pruner sequence. 



time, Atf r ee, must be a least common multiple of the inverses of the chemical shift 
differences 



j 



between the spins connected in the pathway, i.e. 

At free = I'KUijUj |Ay| 



(5.3) 



(5.4) 



where natural number. Thus, all the spins in the pathway are in phase because 

their relative phase is proportional to 2tt, while the others are dephased in each free 
evolution. 



5.2 Numerical simulation on the L- leucine molecule 

Simulations of a local excitation dynamics under "the pruner" sequence were performed 
for different topologies and initial conditions. The topologies, chemical shifts and J- 
couplings to exemplify the operation of the sequence were extracted from amino acid 
moleculeg^j Here, we present the simulations that were done in the L-leucine molecule 
(see fig. |5.2[ ). We find that the pruning is more effective for 

At mix « h/Jf; th (5.5) 



■"■Particular values for leucine in D2O at pH 6.89 were taken from the National Institute of Advanced 
Industrial Science and Technology (AIST) of Japan database (SDBS- 13 C NMRSDBS No. 1142CDS-00- 
770). 
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Figure 5.2: Molecule of L-leucine showing the numbering of the carbons used in the 
spin dynamics calculation. The carboxylate carbon is not considered in the simulations 
because as it resonates at quite different frequency is easily decoupled from the others. 
The 13 C chemical shifts in a 23 T magnet at pH 6.89 are Att 1 /2n = 5.52 kHz, Aft 2 = 
6.325 kHz, Afi 3 = 10.2275 kHz, Afi 4 = 13.6525 kHz and Afi 5 = 5.7825 kHz, and J/K = 
0.04 kHz for all the 13 C- 13 C nearest neighbor couplings. 



where J? at are the J-couplings in the pathway. However, polarization transfer becomes 
much slower than in a regular isotropic mixing sequence. Thus, after numerical explo- 
rations, we found that for a practical experimental situation, a good compromise for an 
enhancement of the signal is found if the time interval for evolution under 7i m ° x is 

At mix ~ ^ (h/Jf; th ) • (5.6) 

This allows one to observe a polarization transfer that is larger than the one obtained 
by the isotropic mixing. In general, for arbitrary chemical shift differences Ay, it is 
difficult to obtain a short time Atf ree . Then, one must find an approximate solution 
to avoid excessive long times. Therefore, the particular enhancement factor depends 
on the number of spins necessary to reach the target: the more spin steps the lower 
efficiency of the overall sequence. In pathways with three spins, the enhancement factor 
in the polarization transfer is approximately three. In the case of a completely enriched 



carbon-13 backbone of an L-leucine residue (fig. 5.3) with the initial excitation at site 
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■ The Pruner sequence 
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Figure 5.3: Local polarization evolution under the pruner sequence at different 13 C sites 
in L-leucine. The initial local excitation at site 1 (the source) is transfered through the 
selected pathway (1 — 2 — 3), to site 3. The free evolution time is Aif ree = 1.2744ms = 



6 x h/Ais « H/Ai2 « 5 x h/A 23 and At iso = 3.3 ms 



lh 
8 J 



. A gain of approx. 270% 



relative to the first maximum of the isotropic evolution is observed at site 3. 



1, one can get a polarization at site 3 which is 270% larger than that obtained with 
7Y^° X alone. The enhancement factor is calculated by comparing the first maximum in 
each evolution. The black lines in fig. 5^ show the polarization evolution on each site 
of the molecule under an isotropic mixing. The red lines show the evolution under the 
pruner sequence. It is remarkable the effectiveness of the pruning manifested in the non 
observable polarization at sites not belonging to the selected pathway, i.e. sites 4 and 5. 
The gray lines show for comparison the polarization evolution under an isotropic mixing 
for a three-spin system composed exclusively by carbons 1, 2 and 3. This is the optimal 
result that one can expect. In the same molecule, by choosing the source at site 2 and 
site 4 as the target (fig. 5.4) one obtains a gain of approximately 152%. There, we can 
observe that the pruning is not working so good as in the previous example. This seems 
to be a consequence of the fact that carbons 1 and 5 do not differ much in chemical shift 
from carbon 2, requiring more cycles of the sequence to dephase appreciably from the 
selected carbons. 
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■ The Pruner sequence 
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Figure 5.4: Numerical evolution of the local polarization under the pruner sequence 
at different 13 C sites in an L-leucine molecule. The 13 C chemical shifts and nearest 



neighbor J-couplings are given in fig. 5.2 The initial local excitation at site 2 (the 



source) is transfered through the selected pathway (2 — 3 — 4), to site 4 (the target). 
The free evolution time is At free = 0.27292 ms = 2 x H/A 24 ~ h/A.23 ~ 2 x h/A 3A and 
At iso = 2.3ms ~ Yjj. A gain of 152% relative to the first maximum of the isotropic 
evolution is observed at site 4. 



5.2.1 Alternative mixing Hamiltonians 



As we mention in section § |5.1[ the pruner sequence stroboscopically interrupts the spin- 
spin interaction to "disconnect" those spins which do not belong to the "selected transfer 
pathway". Here, we will show what happens if we change the mixing Hamiltonian. We 
consider the following mixing Hamiltonians 

^ = E J v ( 2j W '"'I ^l) > ( 5 - 8 ) 



where the first is an XY (planar) Hamiltonian and the second one is a truncated dipolar. 
We mention in previous chapters that the XY Hamiltonian is less "diffusive", i.e., the 
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intensities of the coherences are less attenuated. Moreover, an XY dynamics can be 
solved exactly in some systems (FBE98} IFR99} [DPL02J as we have done in section § 
3.5| to test new theoretical methods. This Hamiltonian is experimentally achievable in 
liquid |MBSH+97j and solid-state jDMFOOl ICCCR06] NMR while the truncated dipolar 
Hamiltonian is the natural interaction in solid-state NMR |Abr61l IEBW911 ISE92] . We 
apply these mixing Hamiltonians to the L-leucine molecule to show the main differences 
in the pruner sequence performance. Note that with the purpose of comparing the 
effect of the different anisotropies in the mixing Hamiltonian, we have kept the chemical 
shift and coupling values of the isotropic interaction in the liquid-state. A real solid- 
state would require knowledge of the 13 C- 13 C dipolar couplings and orientation 
with respect to the external magnetic field. Figures |5.5 and 5J3 show the polarization 
dynamics for a source at site 1 for an XY and a dipolar Hamiltonian respectively. 
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Figure 5.5: Local polarization evolution with the pruner sequence at different 13 C sites 
in an L-leucine molecule under an XY mixing Hamiltonian. The initial local excitation 
at site f (the source) is transfered through the selected pathway (1 — 2 — 3), to site 3. 
The free evolution time is Atf ree = 1.2744 ms = 6 x H/A 13 m H/A 12 ~ 5 x h/A 2 3 and 
At iso = 3.1ms ~ |4. A gain of approx. 420% relative to the first maximum of the 
isotropic evolution is observed at site 3. 



The pruner sequence (red lines) is optimized to transfer polarization to the site 3 (the 
target) through site 2. It is noticeable the performance of our sequence under the XY 



5.2 Numerical simulation on the L-leucine molecule 
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Figure 5.6: Local polarization evolution with the pruner sequence at different C sites 
in an L-leucine molecule under a truncated dipolar mixing Hamiltonian. The initial local 
excitation at site 1 (the source) is transfered through the selected pathway (1 — 2 — 3), to 
site 3. The free evolution time is Atf ree = 1.2744 ms = 6 x h/A ri « ft/A 12 ~ 5 x ft/A 23 
and At iso = 3.1ms ~ |j. A gain of approx. 234% relative to the first maximum of the 
isotropic evolution is observed at site 3. 



Hamiltonian. It almost reproduces the isolated three-spin system dynamics (gray lines). 
The Ising term of the mixing Hamiltonian originates the decrease of the performance in 
the other Hamiltonians because it produces an effective energy difference between the 
sites of the pathway. This detuning does not allow for very high polarization transfers. 
This can be seen in the polarization dynamics of site 1, where for the XY mixing (fig. 



5.5) the polarization transfer is almost complete, whereas in the dipolar (fig. 5.6) and 



isotropic (fig. 5.3) cases, only around half of the polarization is transferred. In spite 



of this, if the goal is the efficiency of the pruning and not the amount of polarization 
transfer, the pruner sequence works quite well even for isotropic or dipolar interactions. 



5.2.2 Step by step pruning of the branches 



We have shown how to "isolate" a selected group of spins. Now, we will show the 
suitability of this technique to choose a first group of spins during a time interval and then 
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to select another group to transfer the polarization from the previous one. In order to do 
this, firstly, we optimize the pruner sequence to transfer the polarization from site 1 to 
site 2. Then, once the polarization arrived at site 2, we optimize the sequence to transfer 
the magnetization to site 3. Fig. |5.7| shows this step by step pruning for different mixing 
Hamiltonians. We observe only the interaction between site 1 (blue line) and site 2 (red 




Figure 5.7: Local polarization evolution at different 13 C sites in an L-leucine molecule 
under a step by step version of the pruner sequence for XY, isotropic and dipolar mixing 
Hamiltonians. The scattered lines show the polarization evolution with this pruner 
sequence, whereas the solid lines show the evolution under the mixing Hamiltonian for 
the five spins (black) and a three-spin chain (gray). The initial local excitation at site 

1 (the source) is transfered to site 2 in the first step of the sequence. Within this step 
the free evolution time is Atf ree = 1.2422ms = h/A\2 and Ati SO = 0.6ms <C j. Within 
the second part of the sequence, the free evolution time At free = 0.2562 ms = fr/A 2 3 
optimizes the polarization transfer between sites 2 and 3 once the polarization arrived 
at site 2. 

line) during the the first step of the sequence where At£f*° = 2ttA/ |A 12 | . It is very 
similar to the isolated two-spin polarization evolution until these spins are disconnected 
in the second step of the sequence by keeping only the interaction between sites 2 and 
3. Within this step, we use At^f ^ = 2irh/ |A 2 3| to optimize the transfer between sites 

2 and 3. Thus, in the figures, one can observe the two-spin dynamics between the sites 
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2 and 3 while the others remain frozen. The figures show essentially a pure two-spin 
polarization evolution. This resembles the qu antum Zeno phase predicted in chapter|4]for 
the swapping operation ADLP06[ IDALP07 which arises when the system-environment 
interaction is much stronger than the system intera ction. There, the Is ing system- 



environment interaction, in the stroboscopic version [ADLP06[ lADLP07aj . acted as a 



measurement process freezing the quantum oscillations. Moreover, the manifestation of 
this phase in a thre e-spin sy stem leads to "isolate" the two spins not directly connected 



to the environment |ALP07) . Here, the pruner sequence has stroboscopic free evolutions 
that act as an environment producing an effective Ising system-environment interaction 
over the "non-selected" spins, isolating the selected ones. This resembles the applications 
of engineered reservoir techniques for fundamental studies of decoherence |MKT + 00[ 
IPazOlj . 



5.3 Summary 



We developed a new NMR pulse sequence, "the pruner", to improve the transfer o f 



polarization through a specific pathway in a system of many interacting spins APL + 07 



The sequence effectively prunes branches of spins, where no polarization is required, 
during the polarization transfer procedure. We obtained a remarkable enhancement, 
higher than 200%, of the polarization transfer with respect to the standard methods. 
Moreover, by changing the mixing Hamiltonian to an XY interaction, we obtain a gain 
of about 400%. Therefore, it is a very practical tool for NMR applications where a 
signal gain is of great importance. However, from a fundamental point of yiew, it is 
the startin g point of an enginee red application of the stroboscopic model ADLP06 
ADLP07a described in se ction § |3.3| which se ems to manifest a spin dynamics within 
the quantum Zeno phase ADLP0"6j ~ DALP07| . Moreover, it can help to go deeper in 



the understanding of decoherence processes and consequently of the environmentally 



induced quantum dynamical phase transition [ADLP06 



DALP07 
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Chapter 6 



Entanglement as a tool for ensemble 
spin dynamics calculations 



The quantum time evolution of a qubit cluster of intermediate size is of great inter- 
est because its potential applications to quantum information processing (QIP) [Kan98, 
BDOO], an d structur al characterization through Nuclear Magnetic Resonance |PLU95[ 



IMBSH+971 ICAL+03l IRSB+05] . Experimental realizations and control of a pure-state dy- 



namics is still one of the major challenges in nowadays quantum physics |QCR04| . There- 
fore, one generally has to deal with ensemble evolutions, as done in previous chapters for 
the study of NMR spin dynamics and in related works jPLU951 IMBSH+971 IRSB+05] . 
There, the dynamics is from an initial mixed-state that, as discussed in the Introduction, 
gives rise to the development of the ensemble quantum computation |VSC04l SSB05J. 
Hence, it is necessary to find a way to calculate or simulate the ensemble evolution effi- 
ciently Theoretical solutions of the ensemble dynamics arise from the integration of the 
Liouville von-Neumann equation (see chapter [2]) or alternative methods like the Keldysh 
formalism, disc ussed in chapter [3] which, for some special cases, yield exact analytical 



results |DFL02t IDPA05 . However, in the most general situations, one has to resort to 
numerical solutions. There are different numerical methods to solve the time dependent 
evolution, such as exact diagonalization, the Suzuki- Trotter Product-Formula (STPF) 
and the Chevyshev Polynomial among others [RM04J. These methods are limited by the 
exponential growth of the Hilbert space dimension with the system size |RM04] . The 
STPF allows one to treat systems of greater dimension than the exact diagonalization 
method which is presently limited to around 20 qubits. The use of the STPF to obtain 
the dynamics of a pure initial state avoids the storage of the Hamiltonian matrix to be 
diagonalized. The computation of an ensemble dynamics with an initial mixed-state, 
requires the calculation of the individual evolutions of each of its components. These 
are of the order of the Hilbert space dimension demanding long computational times. 
Here, in order to overcome this limitation, we take profit of the quantum parallelism 
[SKL02J to calculate the ensemble dynamics through a pure entangled state evolution us- 
ing the STPF. Apparently, the underlying physical mechanism that makes possible these 
efficient simulations is the rapid intrinsic decoherence of highly correlated many-qubit 
systems [KS041 IKS061 ICCCR061 ISPL07] . 

We consider a system of M spins 1/2 to calculate the M-qubits dynamics. We 
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compare the time evolution of a local polarization (experimentally observed by NMR) 
[ZME92bl ILUP98] , between an initial ensemble of local excitations and an initial pure 
entangled state. This entangled state is built in as a superposition of each of the ensemble 
components where the complex coefficients of the linear combination are the thermal 
weights with random phases. We show that the contribution of the coherences of the 
pure initial state to the dynamics can be neglected due to a self averaging property 
arising on the destructive interference of the random phases. To show the relevance of 
the number of independent phases, which are of the order of the Hilbert space dimension, 
we also calculate the same evolution with an initial pure product state where the number 
of independent phases are of the order of the system size. We compare the entangled 
state and ensemble calculations in two contrasting topological configurations: one with 
relevant mesoscopic effects and other where the disorder makes these effects negligible. 



6.1 Ensemble vs. pure entangled state evolution 



Let us start considering a system of M spins 1/2. As we are interested in the evolution 
of a local excitation, we take the ensemble of all states representing the many-spin 
states in the product base along the quantization axis (z), with the n-th site polarized 
(spin up). Each of these states has a thermal weight W{. The local spin dynamics through 
this ensemble can be obtained by calculating a generalization of the auto-correlation 
function P°™(t) |PLU95] 



PZ s (t) 



2-M-ln-M— 1 

i=i f=i 



(* / |e- ifit/,l |* i ) 



(6.1) 



The quantity 



-mt 



represents the probability of finding the n-th site polarized 

in the state \^f) at time t provided that the same site was polarized at t = in the state 
The sum over all possible initial and final states and |^/), gives the amount 
of the z component of the local polarization at time t on the n-th site. As the value of 



the polarization runs from —1 to 1, the eq. (6.1) is properly renormalized. It represents 
the standard thermal mixture calculation obtained with the density matrix [MBSH + 97] 
(see section 5 [2l 
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(6.3) 
(6.4) 

(6.5) 



6.1 Ensemble vs. pure entangled state evolution 



109 



Within the Keldysh formalism, using equations (3.92) and (3.95), we obtain 



jens,K 



(*) 



<*„.|^(0)/-(0)|*eq.> 



(6.6) 



The expression (6.1) involves 2 different dynamics for each of the initial states, see 



fig. 6.1 a). This number is directly related to the dimension of the Hilbert space, and 
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Figure 6.1: Quantum evolution schemes of an ensemble [panel a)] and an entangled 
pure-state [panel b)]. Each represents a simple tensor product state. 



thus the number of required evolutions increases exponentially with M, leading to large 
computational times. Our goal is to extract essentially the same information hidden in 
this ensemble calculation but at a lower cost. To do this we exploit the parallel behavior 
of quantum systems |SKL02j . This is achieved by replacing the ensemble by only one 
pure-state, see fig. 6.1 b), obtained as a linear combination of the components of the 
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ensemble, that is, 



/=1 



E 

/=1 



2 1 

~ 2 

2 



i=l 



(6.7) 



(6i 



Here {a} denotes the set of all the ojj = y/wie m , with y?j an arbitrary phase, involved 



in the initial pure state. Note that the substantial difference between eq. (6.1) and eq 



(6.7) is that the sum in the former is outside the square modulus while in the latter is 



inside. Rewriting eq. (6.7), it becomes 



2 M-1 

/,»=! 



2 1 

~ 2 



e- int/h i^j\e mt/n \^f) 



= 0) + 2j] a^^/le-^l^^-le^l*/), (6.9) 

where we can see that the cross terms make the difference between P^ s (t) and Pnn\t). It 
is evident that the choice of the phases (pi has a relevant role in the equivalence between 
eq. (6.7) and eq. (6.1). For M — ► oo and choosing the (fi randomly dis trib uted, where 



we observe that the cross terms in eq. (6.9) average to 



^purey is an entangled state 1 
zero and 

P™(t) = Pt } (t). (6.10) 

This self averaging property of a randomly correlated pure-state, assisted by the large 
dimension of the Hilbert space, suggests an "intrinsic decoherence" characteristic of 
very complex systems. As we are considering the whole system, one does not expect 
to observe decoherence as defined in chapter [Tj Nevertheless, the observable involves a 
single spin of the system and consequently we cannot distinguish between a mixed-state 
and a pure entangled one, which effectively reflects a loss of information. However, for 
finite size systems, this self averaging property depends on their particular characteristics 
as shown below. When the self averaging property is not satisfied, an extra average over 
initial states is mandatory to force the equality between P™ s (t) and Prm{t). Then, 



1 N a 

w 



(6.11) 



1 If the ifi are randomly chosen, the pure state 
state (non- factorable). 



pure / 



has a high probability to be an entangled 
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where N a denotes the number of different realizations of the sets {a}. One has to perform 
N a evolutions in the last expression, while for an ensemble dynamics one needs 2 . 
Thus, the calculation time is reduced by a factor of N a /2 M ~ l where, for an appropriate 
choice of the {a} coefficients, N a <C 2 M_1 . 



6.2 Application to spin-systems with different cou- 
pling networks 



6.2.1 The systems 



In typical situations of high-field solid-state NMR |Abr611 IEBW9T] . the spin interaction 
Hamiltonian Ti can be expressed by eq. (4.3). One can effectively eliminate the Zeeman 



contribution by working on-resonance in the rotating frame |Sli92j . Then, we can focus 
in the spin-spin interaction, 



H 



M 

£ 

Kj 



.J Z J Z _L lh.. ( f 



where, %/ajj = represents an Ising-like coupling, a^/bij = an XY Hamiltonian, 
ciij/bij = 1 the isotropic one, and aij/bij = —2 a dipolar (secular) Hamiltonian truncated 
with respect to a Zeeman field along the z axis. 

In order to show the potential of the proposal summarized in eq. (6.11), we will 



apply it to two different spin systems which have well differentiated kinds of dynamics 
a) A ladder of spins interacting through an XY Hamiltonian, as shown in fig 



6.3 



a). There, a,- 



0, b 



5j+M/2,i+M/2+l 



br, and &„■ 



-M/2 



by. Here, the dynamics 



6.3, due to the ordered 



presents long lived recurrences, showed in the black line of fig. 
topology |PLU95llMBSH+97] . 

b) A star system in which all the spins are interacting with each other through a 
dipolar coupling, a^/b^ = —2, with a Gaussian random distribution with zero mean 
and a 2 variance for the intensities of the couplings [see fig. 6.3 b)]. In this case, the 
local polarization decays with a rate proportional to the square root of the local second 



moment of the Hamiltonian and no recurrences are observed. The black line of fig. 6.4 
shows the local polarization of this system. 

In order to compare eqs. (6.1) and (6.11), and the dependence of eq. (6.11) on 
the choice of the phases tpi, we calculate the evolution for two types of initial states 
in the infinite temperature limit, i.e. Wi = 1/2 . This limit corresponds to the 
NMR experimental condition [Abr61]. The particular initial conditions are states with 
a local excitation on the n-th site over a background level which is determined by the 
zero magnetization of the others M — 1 spins. These states are used in several NMR 
experiments |ZME92bl IPLU951 lMBSH+97j . 

The pure entangled state has the phases (pi randomly chosen and \l/pure^ becomes 
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Figure 6.2: Panel a) shows a ladder of spins interacting through an XY Hamiltonian 
where the horizontal couplings are given by b x and the vertical by b y . Panel b) contains 
a star system in which all the spins interact with each other through a dipolar coupling, 
dij/bij = —2, with gaussian random distribution with zero mean and a 2 variance. 



The correlation function, eq. (6.11), calculated with this state will be called P^' Na (t). 

The second case is a product (unentangled) state, built with the n-th spin up and 
all the others in a linear combination of spins up and down, with equal probability and 
randomly correlated. Assuming n — 1, we have 



where 



M 



= ir>i<8>n 

m=2 

^(ll) + IT) 



m ' 



(6.13) 



-l<Pm 



tfi. Here, the correlation function (6.11) will be identified as P^ d ' Na (t). 



Note that this state is a particular case of (6.12) with a special correlation on the phases 



6.2.2 Quantum evolution 



The local polarization, P 1 e { ls (t), obtained with eq. (6.1), for the ladder system composed 



of 14 spins is shown in fig. 6.3 with black line. Red and green lines correspond to 
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Figure 6.3: Polarization evolution of the local spin dynamics in a 14-spin ladder system 
of an ensemble compared with an entangled pure-state dynamics. The ratio between 
the x and y coupling in the ladder is given by b y /b x = 1/10. The black line shows 
the ensemble local polarization, P 1 e ° s (t), and the red and green lines correspond to the 
evolution obtained with pv™ ' a anc l p^> <* ^ respectively. The upper panel shows 
the dynamics of eq. (6.11) for N a = 1. The lower panel shows the dynamics for 
ppro , a ^ an( j P^ 7 a (t), where N a is the lower value for which each curve reproduces 
the ensemble dynamics. 
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the temporal evolution of eq. (Oil), Pfl° d,Na (t) and P™ t,Na (t) respectively. The upper 



panel shows the dynamics of eq. ( |6.11 ) for N a = 1. The agreement between P™ s (t) 
and P 1 e f t,1 (t) is excellent, while Pf[° ' (t) has a dynamics far away from that of the 
ensemble. The difference between the dynamics of the two initial pure states is due to 
the different number of independent random phases (uncorrelated phases) of each state. 
In the entangled state, there are 2 M_1 independent random phases that make possible 



the cancellation of the second term in the rhs of eq. (6.9). However, the number of 
independent phases for the product state is M — 1. This implies that there are multiple 
correlations between the phases in the cross terms impeding their self cancellation for 
low values of M. 



The lower panel in fig. 6.3 shows the dynamics for P 1 p 1 rod,JVa (t) and P™ t,Na (t), where 
N a is the lower value for which each curve reproduces the ensemble dynamics. Note 
that the relation between iVp ro d and N ent comes from the equivalence in the number of 
independent phases 

N^^-N^M-l). (6.14) 



This equation for the particular case of fig. 6.3 a 



is 



8192 = N ent 2 M ~ L ~ iV prod (M - 1) = 7969. (6.15) 

The statistical theory of the density matrix [Blu8~T] is based in the random correlation 
nature of a real system to describe it as an ensemble state in the thermodynamic limit. 
Here, we observe that even for small numbers like 14, the equivalence between a randomly 
correlated pure state with an ensemble state remains as a consequence of the dimension 
of the Hilbert space which grows exponentially with M. 



In fig. OA it is shown the same calculations observed in fig. 6J3 but for a star system. 
This system shows a spin "diffusion" behavior for the polarization due to the topological 
complexity yielding the lack of recurrences for long times. For N a = 1, the upper panel 
shows that the agreement between the ensemble dynamics and both pv™ d ' l (t) 
(t) curves is good only for values of the polarization higher than 0.2. For 



of fig. [6.4 
and P, 

times where P^ s (t) is close to zero, the cross terms in eq. (6.9) become relevant and 
both pv™ d,1 (t) and P^' l (t) are different from the ensemble curve. Note that in contrast 
with the ladder case, in this case, the evolutions -P^ od,1 (^) an d P^' 1 ^) are similar. 
For higher values of N a one obtains a better agreement as in the ladder system. It is 
important to note that the difference between P^ s (t), P%™ d ' Na {t) and P^' Na (t) would 
not be appreciable in a real experiment with a typical signal to noise relation. 

The ensemble calculation of eq. (6.1) needs to project the evolution of every 2 M ~ 1 



initial states, 2 13 = 8192 for this 14 -spin system, into the same number of possible final 
states. Instead of this, if one starts with an entangled state or with a random product 
state, the number of initial states gets significantly reduced. Even for the worst case in 
which one needs to do N a ~ 630 averages to obtain a good agreement with the ensemble 



curve (low panel of fig. 6.3), this number represents a small fraction, lower than 8%, of 
the 2 13 initial states of the ensemble. For the case where only one pure entangled state 
is enough to mimic the ensemble dynamics, the number of evolutions is reduced to a 
number around the 0.01% of the 2 13 of the ensemble. This shows the potentiality of the 
method. 
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Figure 6.4: Polarization evolution of the local spin dynamics in a star system of 14 spins. 
The ensemble dynamics compared with an entangled pure-state evolution is showed. The 
dipolar interaction a^jh^ = —2 is given by a random normal distribution with zero mean 
and a variance. The black line shows the ensemble local polarization, P^ s (t), and the 
red and green lines correspond to the evolution obtained with Pv™ d ' Na (t) and P^ ,Na (t) 
respectively. Panel a) shows the dynamics of eq. (6.11) for iV a = 1 and panel b) the 
dynamics for P^ A)Na (t) and P^' Na (t), where N a is the lower value for which each curve 
reproduces the ensemble dynamics. 
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6.3 Summary 

In summary, in order to overcome the limitations of the numerical calculations of an 
ensemble spin dyn amics when one increases the number of spins, we develop a novel 
numerical method |ADLP07bj . It exploits the delicate property of the quantum su- 
perpositions and quantum parallelisms |SKL02j to reproduce the ensemble dynamics 
through a pure entangled state evolution. The method is useful to use the Suzuki- 
Trotter product-formula in an ensemble evolution. We showed that the contribution of 
the coherences of the pure initial state to the dynamics can be neglected due to a self 
averaging property arising on the destructive interferences of the randomly correlated 
pure-state. The underlying physical mechanism that makes possible these efficient sim- 
ulations seems to be related to the observation that higher order coherences in highly 
correlated states decay faster than those of lower order |KS04t IKS061 ICCCR06| ISPL07] . 
This suggests that the self averaging property assisted by the large dimension of the 
Hilbert space could be involved in what is called "intrinsic decoherence". The concept 
developed here can be used in two-ways: on one side, it allows for very efficient dynam- 
ical calculations of common experimental situations where big ensembles are involved. 
On the other side, it sheds light on how to experimentally prepare entangled states for 
specific purposes. 



Chapter 7 

Conclusion and final remarks 



In this thesis we had a dive into quantum dynamics focused in the decoherence phe- 
nomenon. Usually, one wants to manipulate a particular system but inevitably other 
degrees of freedom interact with it changing its dynamics. In order to continue using the 
potentialities of the system of interest, this led us to study how these external degrees of 
freedom disturb the system dynamics While this work was done within the NMR field, 
the workhorse of quantum mechanics, our results could be applicable to many fields 
because they involve fundamental concepts of quantum mechanics. 

In the beginning, we studied the cross-polarization technique for molecular character- 
ization purposes. Through the well known methodology of the generalized Liouville-von 
Neumann quantum master equation, typically used in the NMR field, we incorporated 
the degrees of freedom of the environment. This mechanism was described in chap- 
ter [2] where we solved the spin dynamics of many-spin systems. Interested in the 8CB 
characterization, we in terpreted cross-polarization experiments over this liquid crysta 



2.2 



in several mesophases [CAL + 03j . The 8CB molecule, as we described in section 
could be represented by a three-spin system coupled with a spin-environment. Thus, we 
began with a simple system to understand the cross-polarization dynamics. Based on 
the Miiller, et al. model [MKBE74J, we solved the two-spin dynamics interacting with a 
fast fluctuating spin-bath. We reobtained their solution and then we extended the model 
to a three-spin system interacting with a spin-bath. As this molecule has two possible 
configurations for the heteronuclear coupling, this led us to note that in each space of 
M = ±1/2 there are only two of the three eigenstates that are involved in the dipolar 
transitions that give rise to the oscillations. This was explained as a consequence of the 
symmetry of the system, i.e. the flip-flop can occur only between the carbon and one 
(the symmetric or the antisymmetric) combination of the proton st ates depe nding on the 
relative signs of the heteronuclear couplings {pi = b 2 or bi = —b 2 ) |CAL + 03j . Hence, the 
frequency of the oscillation depends of the different configurations. The experimental 
data were well fitted to the analytical polarization expression concerning the oscillation 
frequency. However, the relaxation process was not well described using a direct exten- 
sion of the MKBE model with an isotropic system-environment interaction. This was 
manifested by the fact that the experimental cross-polarization data of the 8CB molecule 
showed that the rate of attenuation of the oscillations is much faster than the rate of 
polarization transfer from the bath. Consequently, we extended the MKBE model to ob- 
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tain a relaxation superoperator that takes in to account this phe nomenon arising from an 



anisotropic system-environment interaction CAL + 03[ IALP07 . The anisotropy is given 



by the ratio between the interaction rates of the Ising and XY term. We emphasized the 
different roles of the anisotropy of the system-environment interaction on decoherence 
and relaxation processes. The main difference is that while the XY interaction takes the 
system to the total system equilibr ium, the Ising sys tem-environment interaction takes 



it to an internal quasi-equilibrium |CAL + 03[ lALP07j . The introduced anisotropy could 



be explained in the nematic phase by assuming a dipolar system-environment interaction 
Hamiltonian within the extreme narrowing approximation. Nevertheless, in the smectic 
phase the anisotropy is much more pronounced. Hence, we extended the model outside 
the fast fluctuation approximation to take into account slower motions. In this way, a 
better agreement with the experimental observa tions was obtained without resorting to 



other mechanisms which operate in both phases CAL + 03 . This detailed spin dynamics 



calculations allowed us to obtain separately the homonuclear and hetero nuclear dipolar 



couplings in CH 2 systems which constitute the 8CB molecule [GALOPS] . We tested the 



reliability of the results with an experimental direct determination of the heteronuclear 
couplings using cross-polarization under Lee-Goldburg conditions. 

Within a fundamental point of view, we observed that the solutions based on the 
MKBE model did not describe the quantum quadratic short time behavior nor the fea- 
tures displayed in fig. 1.6 within the b <C ^se region. In order to overcome these 



limitations, a further improvement was done including non-secular terms of the system- 
environment interaction to extend the solution. This solution describes well the experi- 
mental short time behavior and led us, together with the spin dynamics analysis within 
the Keldysh formalism (chapter [3|, to enrich our perspectives. 

Inspired in the generalized Landauer-Buttiker equation within the Ke ldysh form al- 



ism, we obtained that the initially phenomenological stroboscopic model [ADLP06 



in 



the continuous form, led to the non-secular solution within the generalized quantum 
master equation. This induced us to obtain analytical results of the spin dynamics 
from microscopic derivation and c haracterization of the system-environment interac- 



tion [DPA051 |DALP07[ lADLP07aj . To do that, we solved the Schrodinger equation 
within the Keldysh formalism for an open system under the wide band regime in the 
environment (fast fluctuation approximation). Within this formalism, assisted by the 
Jordan- Wigner transformation that maps a spin system into a fermion one, an exact 
solution for an XY linear chain is obtaina ble. This allowed us no t onl y to include mem- 

but it enabled us 
To 



3.5 



3.4 



ory effects within the spin-bath [DPA05], discussed in section 
to test approximation methods to solve the more complex systems of section 
describe the spin-bath under the fast fluctuation approximation that leads to interac- 
tions local in time, we resorted to the Wigner time-energy variables. This allowed us 
to transform the density function expressed in the Danielewicz integral form into a gen- 
eralized Landauer-Butt iker equation. By ap plying this technique to a two-spin system 



coupled to a spin-bath |ADLP06[ lADLP07a] , we improved the results obtained through 
the secular approximation within the standard density matrix formalism. Further on, 
we effectively symmetrized the system-env ironment in teractions virtually transforming 
them into a spatially homogeneous process ADLP07aj. This involved a uniform system- 
environment interaction rate that leads to a simple non-hermitian propagator while 
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the original multi-exponential decay processes were recovered by an injection density 
function. This method led us to interpret the phenomenon of decoherence and the 
system-environ ment intera ction within a special view: The environment as a measure- 
ment apparatus 



ADLP06 



There, the decoherence takes another perspective: The loss 
of information of the system entity, i.e. how the isolated system dynamics is degraded in 
a characteristic time given by the fictitious homogeneous decoherence time. The homog- 
enization procedure enabled the microscopic derivation of the stroboscopic model for the 



system-environment interaction |ADLP07a) . While, the stroboscopic process may not 



be the best description of reality, it provid es an optimal numeric al algorithm to calculate 



quantum dynamics in discrete time steps [ADLP06[ lADLP07a 



While the dynamics obtained through the Keldysh formalism are reproduced by the 
non-secular solution within the generalized quantum master equation, they derive from 
a microscopic model of the entire system (system plus environment) and the final state 
must not be hinted beforehand. The Keldysh formalism gives us another perspective to 
discuss about the physics of the quantum time evolution and the interpretations of the 
approximations made. For example, the arising of the quantum Zeno effect shown in the 
decoherence time, l/r^ oc {b/h) 2 T SE , could be interpreted as a "nested" Fermi golden 

rule rate emphasizing the non-perturbative nature of the result ADLP06 . 

The manifestation of the quantum Zeno effect led us to a novel interpretation of 
previous experiments |LUP98j shown in fig. 1.6 In chapter El we found experimental 



evidence that environmental interactions can drive a swapping gate (two-spin system) 
through a Quantum Dynamical Phase Transition towards an over-damped or Zeno phase 
ADLP06"] . The NMR spin swapping experiments in a 13 C- 1 H system enable the iden- 
tification and characterization of this phase transition as a function of the ratio br^/h 
betwee n the internal and sy stem-environment interaction. The developed microscopic 
model DALP07[ IADLP07a for the swapping operation describes both phases and the 
critical region wi t h great det ail, showing that it depends only on the nature of the inter- 



action 



ADLP06. ADLP07a 



In particular, it shows that the phase transition does not 
occur if the system-environment interaction gives isotropic interaction rates, T/z — Txy- 
Within the standard approximations typically used to solve the generalized quantum 
master equation, one tends to think that r zz = T X y corresponds to an isotropic system- 
environment interaction. However, a careful microscopic derivation within the Keldysh 
formalism showed us that this is not necessarily correct. It is also important to men- 
tion that the spin-bath occupation factor modifies the Tzz rate. For the description of 
the quantum dynamical phase transition, it is crucial to distinguish between inputs and 
output parameters. One can visualize applications that range from tailoring the envi- 
ronments for a reduction of their decoherence on a given process to using the observed 
critical transition in frequency and decoherence rate as a tracer of the environment's 
nature. For example, we extended the model to a 3-spin system to show that beyond a 
critical region two spins become almost decoupled f rom the environment oscillating with 



the bare Rabi frequency and relaxing more slowly [ALP07 



Inspired in the stroboscopic model and the arising of the environmentally induced 
quantum Zeno phase, we developed a new NMR pulse sequence. "The pruner" strobo- 
scopically interrupts the system evolution to improve the transfer of polarization through 
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a specific pathway in a system of many interacting spins |ADL + 07] . The sequence effec- 
tively prunes branches of spins, where no polarization is required, during the polarization 
transfer procedure. We obtained a remarkable enhancement of the polarization transfer 
with respect to the standard methods. Therefore, it is a very practical tool for NMR 
applications where a signal gain is of great importance. However, from a fundamental 
point o f view, it is the start ing point of an engine ered application of the stroboscopic 
model ADLP06. ADLP07a described in section 8 3.3 



dynamics within the quantum Zeno phase |ADLP06[ IDALP07 



which se ems to manifest a spin 
Moreover, it can help 



to go deeper in the understanding of decoherence processes and consequently of the 



environmentally induced quantum dynamical phase transition [ADLP06[ IDALP07 



Finally, in order to study the spin dynamics of la rger system s, we developed a novel 
numerical method that exploits quantum parallelism |ADLP07b . This provides a tool to 
overcome the limitations of standard numerical calculations of ensemble spin dynamics 
for high number of spins. Hence, this numerical method constitutes the starting point 
to extend the study of the quantum dynamical phase transition to larger systems and 
how it is involved in the irreversibility phenomena. Moreover, the observation that the 
contribution of many coherences of an entangled randomly correlated state to the dy- 
namics can be neglected due to a self averaging property, opens a fundamental question: 
Is this self averaging property involved in what is called "intrinsic decoherence"! 
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